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ABSTRACT 

This  report  provides  rcdor  designers  with  background  information  on  how  the  idenlized 
free-spoce  radar  theory  must  be  modified  to  toke  account  of  reflections  from  the  eorth's 
surfoce.  It  is  primarily  concerned  with  problems  that  arise  in  designing  systems  thot  are 
oirbome  and  contoins  discussions  of  the  following  topics:  effects  of  reflections  from  the 
earth  on  the  signal  received  from  an  elevated  target,  implications  of  these  effects  for 
detection  and  parameter  estimation,  effects  of  surface  roughness,  techniques  for  the  re- 
(kiction  of  clutter,  scattering  from  the  occon,  and  spherical -earth  geometry. 
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INFLUENCE  OF  THE  EARTH*S  SURFACE  ON  RADARt 


INTRODUCTION 

The  purpose  of  this  report  is  to  provide  radar  designers  with  background  material  on  how 
the  idealized  free-space  radar  theory  must  be  modified  to  take  account  of  the  earth's  surface. 
This  surface  constitutes  an  important  component  oi  the  propagation  medium  and  plays  a  signifi* 
cant  role  in  determining  the  extent  to  which  a  given  design  will  prove  successful. 

For  the  radar  designer,  the  presence  of  the  earth's  surface  has  two  important  implications. 
First,  due  to  the  convexity  of  the  surlace,  for  any  given  height  of  the  radar  antenna,  a  portion 
of  the  space  above  the  earth  will  lie  in  shade tv.  Second,  due  to  the  reflective  properties  of  the 
surface  materials,  a  certain  amount  of  the  energy  incident  on  the  surface  will  be  reflected  and 
scattered  back  into  the  space  above  it.  Both  effects  must  be  evaluated  in  order  to  predict  a 
radar's  performance. 

Assuming  that  the  earth  exists  in  an  atmosphere  with  a  constant  index  of  refraction  (e.g., 
free  space)  so  that  electromagnetic  rays  can  be  drawn  as  straight  lines,  one  can  represent  the 
problem  imposed  by  the  earth's  convexity  as  shown  in  Fig.  1.  When  a  target  is  above  the  horizon 


line,  it  will  be  illuminated  directly.  When  it  is  below  the  horizon  line,  it  will  be  illuminated 
only  to  the  extent  that  diffraction  takes  place.  Thus,  if  the  function  of  the  radar  requires  good 
coverage  near  the  surface  of  the  earth,  one  must  either  {i)  try  to  raise  the  height  of  the  radar 
antenna,  thus  extending  the  direct  illumination  region,  or  (2)  try  to  take  advantage  of  diffraction 
Each  of  these  techniques  for  extending  low  altitude  coverage  has  certain  advantages  and  dis¬ 
advantages,  and  each  imposes  strong  constraints  on  the  choice  of  radar  parameters.  For  ex¬ 
ample.  to  make  use  of  diffraction,  one  must  go  to  very  low  frequencies.  On  the  other  hand, 
elevating  the  antenna  an  appreciable  amount  limits  the  size  of  the  antenna  that  can  be  employed. 


tfhU  report  wai  originolty  intended  for  inclMion  fn  o  book,  it  it  inwed  at  o  Technical  Report  becoute  the  pub' 
licotion  of  this  book  hot  been  unovoidobly  deloyed.  The  informotlon  on  which  this  report  it  bated  wot  obtained 
prior  to  1962  and  no  attempt  hot  been  mo^  to  update  It. 
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Fig. 2.  Forword-  and  bockscatterad  energy. 
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U  uno  takes  into  account  the  variability  of  the  index  of  refracion,  then  a  third  possinui' \  ;u  iset: 
(3)  one  can  atlvnf.pt  to  entend  the  iov,'  altitude  coverage  by  n^.aking  use  of  atmospheru  berni.ng 
and  S'  at;ering  'rh;.s  again  implies  strong  limitations  on  the  iho.i  e  of  radar  par am'  ters  iiia.s- 
much  O'J  c?ie  must  choose  parametr-rs  for-  which  th<-re  is  suffit  I'-nt  ticnding  aiiri  sratteting  In 
general.  wliH-h  of  these  techniques  will  hr  most  useful  will  it*  pend  on  the  specifu  tumtion  for 
which  the  radar  is  being  design«*d  In  this  report,  it  will  be  assumed  th.ii  tne  ranar's  fuii.  tion 
IS  to  de*«'ct  and  track  aweraft  at  the  greatest  possible  distances,  and  that  ilicse  aiii  raft  may  be 
flying  at  low  altitudes  as  well  as  high  ones  b'or  this  furKiion,  ilie  solution  whuli  has  been  found 
most  sut  ccasful  is  that  of  elevating  the  antenna,  and  it  is  this  typi-  of  larlar  toward  whn  h  the 
present  discussion  will  be  ontntt  d  More  spec  ifu  ally,  th<-  diseiission  will  tie  limited  to  <  on- 
sidert.tioiio  which  aris.-  in  designing  an  airborne  early-warning  and  lontrol  radar  ,  and  Ui  j  heiiom- 
ena  which  occur  inside  the  direct -illumination  region,  Onl;.  Iho-sr-  problems  will  be  (oiisuieted 
which  involve  energy  transmitted  by  the  radar  itself  (as  opposed  to  ener  gy  t  ransmiitr-d  by  laitsi.-ie 
sources),  and  spe./ial  r  onsideratioa  will  tie  given  to  ‘he  casr-  in  whir  h  the  radar  anti  tar  get  are 
operating  over  wat«T  a.s  orposed  to  land.  It  w.ll  also  be  assumml  that  atmospheric  conditions 
are  such  that  the  bending  of  the  electromagnetic  rays  due  *o  variations  in  the  index  of  refrartion 
can  be  rectified  by  usi^g  an  equivalent  earth's  radius. 

Wjthir:  the  direct  illumination  region,  the  most  important  phenomenon  to  be  consi.k  rerl  in 
connection  with  thr-  earth's  surface  is  that  of  scattering.'^  For  many  mirth  surlaci-  materials  and 
many  frequencies,  only  a  small  portion  of  the  incident  energy  will  be  absorbed  and  most  of  it 
will  be  scattered  back  into  thr  space  above  the  earth.  It  is  the  existent  e  of  this  siatir  red  rmergy 
that  has  caused  ‘.he  direct -illumination  region  to  be  referred  to  as  the  "interlerem  e  region"  In 
general,  at  any  given  point  on  the  earth’s  surface,  this  scafered  enr-rgy  will  propagatr-  in  all 
direetiotis.  The  only  directions  of  importance  to  the  radar.  hov;ever.  are  the  direction  toward 
the  target  {the  "forward-scattered  energy")  and  the  direction  back  to  the  radar  antenna  (the  "back- 
scattcred  energy"  or  "clutter").  A  schenriatic  illustration  of  the  foi’W'ard-  and  backscattered  en¬ 
ergy  is  shown  in  F'ig.  Z.  As  far  as  obtaming  infe-rmation  on  the  aircraft  target  is  concerned, 
whereas  the  forward-scattered  energy  is  in  many  ways  an  advantage  for  the  radar,  the  back- 
scatterel  energy,  being  a  form  of  interfei  ence,  is  inherently  a  disadvantage.  The  relative 
amounts  of  enei  gy  scattered  backward  and  forward  will  depend  upon  the  detailed  structure  of  the 
region  over  which  the  '•adar  is  operating,  the  parameters  of  the  radar  system,  and  the  geometry. 

•Assuming  that  the  surface  is  perfectly  smooth,  one  finds  that  the  effect  of  the  surface  on  the 
target  signal  is  to  create  four  possible  round-trip  paths  oyer  which  the  energy  can  travel  from 
radar  to  target  and  back  to  the  radar.  These  four  paths  are  illustrated  in  Fig.  3.  .Making  use 
of  images  (and  assuming  fo'-  simplicity  that  the  earth  is  llat),  one  can  Interpret  these  four  com¬ 
ponent  target  signals  and  the  smooth-earth  clutter  .signal  as  shown  in  Fig  4.  -According  to  this 
interpretation,  the  siUiation  is  similar  '.o  one  in  which  there  are  two  targets  and  two  sources  in 
free  space  with  the  targets  and  sources  constituti.ng  slaved  pairs  and  with  no  scattering  between 
targets.  .As  the  surface  becomes  rougher,  the  images  tend  to  spread  out. 

F'or  certain  purposes,  the  roles  of  the  forward-  and  backscattered  energy  in  modifying  the 
received  signal  can  best  be  visualized  in  *erms  of  the  filter  concept  i.etting  X{u))  and  A'(w)  be 
the  complex  spectra  of  the  transmitted  and  jvc*  ived  signals,  and  assuming  that  the  system  con¬ 
sist  mg  of  the  radar  antenna  and  su'-rounding  space  can  be  regarded  as  a  fixed  linear  filter  with 
a  tran.sfer  furtet  jo-!  S{w),  one  can  relate  A'fui)  lo  Xiw)  by  the  equation  Y(u;)  S(ui)  X(w)  (see  F'ig.  5). 

t Throughout  this  report,  the  word  "tcatiaring"  wU>  be  'tted  to  ref«r  to  all  ancrgDr  which  is  not  obsarbod. 
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Letfing  FUv)  bft  the  space  filter  under  the  assumption  that  there  is  no  ear‘h  present,  and  C(w) 
the  space  filter  under  the  assumption  that  tl  ere  is  no  target  present,  one  can  write  S(w)  ■ 

Fiv)  E(w)  +  C{o)),  where  Eiu;)  is  the  earth-mcdification  factor  for  the  target  signal.  In  this 
language,  the  introduction  of  the  earth  into  the  radar  problem  is  equivalent  to  the  filter  trans¬ 
formation  F{w)  •'  F(w)  E{ti))  +  C(w),  and  the  task  of  describing  the  influence  of  the  earth  on  the 
recei’'8d  signal  Yiw)  is  equivalent  to  the  task  of  describing  the  functions  E(w)  and  C(uj;. 

The  discussion  in  this  report  will  be  divided  into  four  main  sections;  I.  Influence  of  a 
Smeeth  Earth;  II.  General  Remarks  on  Influence  of  a  Rough  Earth;  III.  Scaitr  ring  From  ihc 
Ocean:  and  IV,,  Spherical-Earth  Formulas,  The  main  task  in  Sec.  I  will  b''  to  show  how  the 
free-spvice  radar  equation  must  be  modified  to  take  ic’o  account  the  presenc,  of  the  earth's  sur¬ 
face  under  the  idealized  assumption  of  spec’Iur  reflection,  and  to  point  out  some  of  the  imph- 
CEtions  of  this  modification  for  various  radar  functions  such  as  de»ection  and  height -finding. 
Inasmuch  as  the  interference  caused  by  the  clutier  signal  in  the  snnooth -earth  case,  has  no  im¬ 
portance  for  the  airborne  early-warning  problem  (the  source  of  the  clutter  signal  being  confined 
to  the  specular  point  directly  beneath  the  radar),  t,he  discussion  in  Sec.  I  will  be  devoted  ex¬ 
clusively  to  the  target  signal.  Ir  Sec.  II,  the  assumption  that  the  eari.h  is  smooth  will  be  elim¬ 
inated,  and  the  effects  of  surface  roughness  on  the  radar-design  problem  will  be  considcrec 
without  reference  to  the  particular  scattering  characteristics  of  the  rough  surface.  In  Sec.  Ill, 
the  rough  surface  will  be  specialized  to  that  of  the  ocean  ano  a  detailed  description  wil'  be  given 
of  the  scattering  results  for  this  surface.  In  both  Secs.  11  and  III,  attention  will  be  focused  on 
the  problem  of  backsoattering,  the  modifications  of  the  forv/ard-scattered  energy  being  consid¬ 
ered  only  vary  briefly  Finally,  Sec.  IV  will  contain  a  list  of  formulas  relating  the  various 
geometric  variables  whose  values  are  needed  in  order  to  p*  eaici  a  radar's  performance  once 
the  physical  paratriPters  of  the  situation  have  been  adequately  determined. 

I.  INFLUENCE  OF  A  SMOOTH  EARTH 

As  already  stated,  if  it  is  assumed  that  the  earth's  surface  is  a  smooth,  homogeneous 
sphere,  the  total  returned  signal  from  an  elevated  target  will  consist  of  the  sum  of  four  com¬ 
ponent  signals  {DD,  II,  DI,  ID),  each  of  which  corresponds  to  a  particular  round-trip  path  from 
radar  to  target  and  back  to  the  radar  again.  This  signal  will  depend  upon  (a)  the  transmitted 
signal,  (b)  the  radar  antenna  pattern,  (c)  the  target  cross  -section  pattern,  (d)  the  reflective 
properties  of  the  surface,  and  (e)  the  relative  locations,  orientations,  and  motions  of  the  antenna, 
target,  and  earth. 

A.  Vertical -Plane  Ccordinates 

If  the  orientations  of  the  antenna  and  target  are  held  constant  with  respect  to  some  fixed, 
arbitrary  reference  system,  and  the  locations  of  the  antenna  and  target  ai  e  constrained  to  a 
fixed  vertical  plane,  all  the  geometric  variables  entering  into  the  equation  for  the  received  sig¬ 
nal  will  be  determined  by  the  locations  of  the  antenna  and  target  within  this  vertical  plane. 
Ordinarily,  in  order  to  specify  these  locationd,  one  would  need  to  use  f  ir  cviordinates  —  two 
for  the  antenna  and  two  for  the  target.  However,  since  it  has  been  assumed  that  the  earth's 
surface  is  a  smooth,  homogeneous  sphere,  the  n^turned  signal  will  be  independent  of  the  par¬ 
ticular  portion  of  the  surface  over  which  the  radar  is  operating  and  will  thus  be  determined  by 
only  three  coerdinates.  Four  vertical -plane  coordinate  systems  which  are  of  frequent  use 
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(see  Fig.  6)  are  (h,  H,  d),  (h,  H.  R),  (n,  y,  R),  and  (h.  A.  R),  where  h  =  radar  height,  H  =  target 
height,  d  -  ground  distance  between  radar  and  target.  R  =  length  of  direct  path,  y  ~  angle  be¬ 
tween  direct  path  and  horizontal  at  radar,  and  A  =  pathlcngth  difference  between  direct  and  in- 
uireoi  paihs.  Whereas  the  first  of  these  systems  is  the  most  natural  one  from  the  vantage  point 
of  the  earth’s  surface,  the  third  and  fourth  are  the  most  natural  from  the  vantage  point  of  the 
radar;  all  the  variables  h,  y,  A,  and  R  can  be  measured  directly  in  terms  of  the  signal  re¬ 
turned  from  the  target. 


Fig.  6.  Vertical -picme  coordinates. 


Assuming  for  purposes  of  illustration  that  the  earth  is  flat,  one  has 


R  +  (H 


A  =  (R^  +  4hH)^/^  -  R 


y  -  sin"^  ( ^-) 


ii 


(A/2)^  h^  -  {A/2)^ 


Assuming  h/d  and  H/d  are  small,  one  obtains  the  approximations 

2hn 


R  ^  d 


H  -  h 
r  =  — - 


A  - 


H  = 


R 


A/2 


[h^-(A/2)V''^ 


(1) 


(2) 


Spherical-earth  relationships  for  these  vt.riables,  as  well  as  other  variables,  can  be  found  in 
Sec. rv. 


B.  Complex  Spherical -Earth  Reflection  Coefficients 

In  order  to  give  a  quantitative  description  of  the  component  signals  along  the  paths  II.  DI. 
and  ID  (see  Fig.  3).  it  is  necessary  to  compute  the  effect  of  ihe  reflection  process  on  the  am¬ 
plitude  and  phase  of  the  incident  rays.  In  general,  the  electromagnetic  properties  of  the  air 
and  earth  that  need  to  be  evaluated  for  this  computation  are  ihe  permittivity  c,  the  conductivity 
cr.  and  the  permeability  In  the  following  discussion,  it  will  be  assumed  that  p.  c,  and  a 
for  the  air  are  the  same  as  for  free  space,  and  that  p  for  ‘.he  earth  is  the  same  as  for  free 
space. 


1  Plane  Reflection  Coefficients 

Assume  that  the  incident  wave  is  a  plane  wave  of  wavelength  y  and  that  the  earth  is  flat. 
For  any  vector  z,  let 


6 


h(  z  )  -  component  of  z  normal  to  plane  of  incidence  (the  plane 
normal  to  the  reflecting  surface  which  contains  the  in¬ 
cident  ray) 


v{  z  )  component  of  z  in  plane  of  incidence 
h  (v(?'I  =  component  of  v(?)  parallel  to  reflecting  surface 
V  {v{z'I  -  component  of  v(z)  normal  to  reflecting  surface. 


Also,  let 


&  =  complex  dielectric  constant 

=  — — ihOku  where  c  and  a  are  evaluated  for  the  earth  sur- 
^o 

face  material.  Cq  is  the  value  of  f  for  free  space,  and  X  is 
given  in  meters  and  a  in  mhos/meter 

n  =  unit  propagation  vector 

4  =  grazing  angle 

K  -  electric  field  vector 

H  ■  magnetic  field  vector 

-  reflection  coefficient  for  horizontal  polarization 

-  reflection  coefficient  for  vertical  polarization. 

Using  the  subscripts  1  and  2  to  denote  the  incident  and  reflected  waves,  one  can  illustrate  the 
cases  of  horizontal  and  vertical  polarizations  as  shown  in  Fig.  7(a-b).  Horizontal  polarization 
is  defined  by  the  condition  h{Kj)  -  and  vertical  polarization  by  vdl^)  -  K^.  Making  use  of  the 
above  definitions  and  assumptions,  one  can  show  that  the  reflected  wave  will  also  be  a  plane 
wave  of  wavelength  X  and  that 


11(^2)  =  h{n2)  = 

vfK,)  ^  -Ty  f  n2  5<  [iTj  X  vfFj))]  v(n2)  = 

^  .sintf-(£-cos"ti);/^.^<^2>.^f^<»2>l  . 

^  sinrf  +  {&  -  cos^  4)^^^  h(Ej)  v 

&sini  h(Hj)  v  fvfE^)] 


-r^  fn2X  In,  X  v(Hj)l) 


-h  (v(H2)) 
h  [v(H,)l 

-Ti  (v(E2)) 

h  [v(E,)| 


(3) 


Note  that  if  t'  =  =  '*■>  is  small,  then  vv  =  v  so  that  v(E,)  =  F  v(E.)  and  v(H,)  =  F.  v(H,).  If 

iz  zviznl 

!6i  »  1  (which  is  usually  the  case  except  for  very  dry  ground),  the  formulas  for  Fj  and  F^ 
reduce  to 


sinv>  - 

li  + 


sin 


^  sin  j  -  1 


(4) 


Note  also  that 

lim  F.  =  lim  F  =  -1 

lim  F,  -  - 1  and  lim  F  -  1 

c  h  f. 


when  ('<"*> 


when  !i  >  0 


(5) 


The  variations  of  and  with  frequency  and  grazing  angle  for  the  case  of  sea  water  are  shown 

in  Fig.  8(a-b).  For  a  detailed  discussion  of  these  plane  reflection  coefficients,  the  reader  is  re- 
2 

ferred  to  Stratton. 


2.  Divergence  Factor 

In  the  above  equations  for  the  reflection  coefficients,  it  is  assumed  that  a  plane  wave  is  ' 
cident  on  a  plane  surface.  A  more  realistic  assumption  for  propagation  over  the  eartn  is  that  a 
spherical  wave  is  incident  on  a  spherical  surface.  A  second-order  "geometrical-optical”  approxi¬ 
mation  has  been  worked  out  for  this  case  by  van  der  Pol  and  Bremmcr.^  They  liave  shown  that 
the  reflection  coefficient  p  for  this  case  can  be  approximated  by 

p  =  Jr  (6) 

where  J,  the  divergence,  is  a  purely  geometric  factor  which  describes  the  extra  divergence  of 
a  beam  of  rays  due  to  the  reflection  from  a  spherical  surface.  More  specifically, 

3)=  lim  (7) 

n-*o  ^ 


where  Q'  and  Q  are  the  corresponding  curved  and  flat  earth  cross  sections  of  the  reflected 
beam  (see  Fig.  9),  and  il  is  the  solid  angle  subtended  by  Q.  Letting  •  denote  the  equivalent 
earth's  radius,  one  notes  that 

0  <  J  <  1  for  all  h.  H.  <>.  . 

J-^las.-**  or  h-*0  or  H-*0 

S  -*  0  as  0  .  (81 

Unless  4i  is  quite  small,  for  most  practical  purposes,  J  can  be  approximated  by  unity.  Ex¬ 
plicit  formulas  for  computing  J  can  be  found  in  Sec.  IV. 


Fig.  9.  Dlvergsnc*  foctor. 
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C.  Description  of  Received  Signal  in  Frequency  Domain 

Assume  for  simplicity  that  (1)  the  transmitting  and  receiving  antennas  are  identical;  (2)  the 
polarization  of  ihe  antenna  is  pure-horizontal  or  pure-vertical;  (3)  the  antenna  and  target  are 
motionless;  (4)  the  hypothesis  of  the  reciprocity  theorem  is  satisfied.^  Let 

t  =  time 

w  =  ang>ilar  frequency 

c  =  velocity  of  light 

t  =  direct  path  D  or  indirect  path  I 

g(f )  ~  complex  amplitude  gain  of  antenna  evaluated  along  t 

a(t^,  t^)  =  complex  scattering  length  of  target  evaluated  for  incident  wave 
along  it  and  reflected  wave  along  (the  target  cross-section 
a  is  related  to  3  by  a  =  4*  j  s  |  ^) 

X  =  complex  coefficient  of  exp(iu>t)  on  transmission 

Y  =  complex  coefficient  of  exp(iut)  on  reception 

Yp  =  free-space  component  of  Y  (i.e.,  component  of  Y 
resulting  from  path  OO) 

F  =  Yp/X  =  free-space  factor 

E  =  Y/Yp  =  earth-modification  factor. 

The  quantities  g,  s,  X,  Y,  Yp.  F,  E.  and  p  are  all  functions  of  the  frequency  w,  and  all  but  X 
are  also  functions  of  the  geometry.  In  order  to  simplify  the  notation,  these  dependencies  will  be 
made  explicit  only  where  necessary.  Assume  also  that  (5)  s(D,  D)  s(I,  I)  =  s^(D,  I).  [This  assump¬ 
tion  will  be  sure  to  be  satisfied  provided  the  angle  between  D  and  I  at  the  target  is  small  with 
respect  to  the  angular  variations  in  the  target  scattering  length,  for  then,  approximately,  s(I,  I)  = 
s(D,D)  =  s(D.I).l 

Applying  assumptions  (1)  through  (5)  (and  assuming  a  refer*.. .ce  impedance  level  of  one  ohm 
so  that  the  power  is  given  by  the  square  of  the  amplitude  of  the  voltage),  one  can  show  that 

Y  =  XFE  =  YpE  (9) 

where 


F  =  exp{- 


E  =  [1  +  K  exp(- 


and  Z  and  K  arc  defined  by 


Z  =  I  g^(D)  s(D.  D) 


The  term  K  exp(-iwA/c)  occurring  in  E  is  nothing  more  than  the  ratio  of  the  complex  signals 
corresponding  to  the  paths  I  and  D.  The  DD  componeni  Yp  of  Y  is  obtained  from  the  term  1 
in  E.  the  II  component  from  the  term  [K  exp<-iwA/c))^.  and  the  DI  and  ID  components  (which 
are  equal)  from  the  term  2K  exp(-iwA/c).  When  E  =  1  (i.e.,  K  =  0),  Y  reduces  to  the  free- 
space  equation  Y  =  Yp.  The  precise  role  of  each  of  the  assumptions  (1)  through  (5)  in  obtaining 
the  above  equation  will  become  apparent  in  Sec.  I-I  where  this  equation  will  be  generalized. 

Writing  any  complex  number  z  as  |zi  exp{i(f>  )  and  defining  q  by  q  =  -  wA/c,  one  finds 

that  the  quantities  lYf  and  <0^  are  given  by 


tFor  o  ditcustion  of  fh*  raciptoclty  thoortm,  tf,  for  •xompie,  Schclkunoff  and 
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UU) 


wh^ro 


i..i2  •..iZ 

1T|  =  :A1  Il'l  IE.| 


f|2  . 


lK|^  ^  (1  +  1K|^  +  2  |K|  cosq)^ 


and 


</>y  =  ■>■ 


(11) 


where 


o„  =  2tan-^^4r7^  ) 

£<  1  +  !Ki  cosq 


Assume  now  that  the  physical  properties  of  the  antenna,  target,  and  earth,  and  the  orienta¬ 
tions  of  the  antenna  and  target,  are  held  fixed,  and  that  the  locations  of  the  antenna  and  target 
are  constrained  to  a  fixed  vertical  plane.  With  this  assumption,  all  the  quantities  entering  into 
the  equation  for  Y  can  be  regarded  as  functions  of  h,  R,  A,  and  u).  If  h,  R,  and  A  are  held 
fixed  and  w  is  allowed  to  vary,  the  above  equations  constitute-  a  description  of  the  frequency 
spectrum  of  the  received  signal  for  the  given  h,  R,  and  A.  liquation  (9)  gives  the  complex 
spectrum  Y(a;),  liq.  (10)  gives  the  energy  spectrum  iY(a,')|^,  and  Rq.dl)  gives  the  phase  spec¬ 
trum  ^  ts  held  fixed  and  h,  R,  and  A  arc  allowed  to  vary,  the  above  equations  con¬ 

stitute  a  description  of  how  the  received  signal  dept*nds  upon  the  positions  of  the  antenna  and 
target  within  the  vertical  plane  for  the  given  a:. 

Assume  in  addition  that  attention  is  confined  to  a  region  of  (h.  R,  A,  u.-)  space  in  which 
!  K(h,  R,  A,  «)  I  varies  slowly  in  comparison  to  cos  (q(h,  R.A,  a>)|.  Within  this  region,  the  earth- 
modification  factor  K  will  be  locally  periodic  in  q  of  period  2t.  Graphs  of  |K(q)|  and  e>jj(q) 
for  q<  ir  and  various  fixed  values  of  jKl  are  shown  in  Kig.  lO(a-b),  It  is  easy  to  show  that 

max  |E(q)|^  =  (1  +  |k1)"^  occurs  at  q  =  Znx 


min  |K(q)|^  •-  (1  -  1k!)”^  occurs  at  q  =  (Zn  ■*  1)  ir 
ave  !K(q)i^  =  1  -^  |  Kl"*  -t  4  | 


|E(q)|  =  1  and  V’p(q)  -  0  when  ]K1  -  0 
iE(q)|^  -  16  cos^  (q/3)  and  <?p(q)  -  q  when  |Ki 


^  1 


I'Kqll'^  —  jKl"*  and  V’,.-(q)  —  Zq  when  1K|  -*  « 


(12) 


If  attention  is  further  confined  to  a  region  of  (h.  R,  A.  w)  space  in  which  <Pp(h,  R,  A,  w)  is  slowly 
varying  in  comparison  to  uiA/c,  then  q.  and  therefore  E,  will  be  locally  symmetric  in  ui  and 
A.  This  implies  that  u  and  A  play  dual  roles  in  the  modification  of  the  free-space  signal  im¬ 
posed  by  the  presence  of  the  earth.  (They  do  not,  of  course,  play  dual  roles  in  the  total  re¬ 
ceived  signal  Y  since  they  do  not  occur  symmetrically  in  Y^,.)  Letting  A  remain  fixed,  E  is 
periodic  in  a)  of  period  Zirc/A  with  a  phase  determined  by  j>j^.(A).  Letting  u  remain  fixed,  E 
is  periodic  in  A  of  period  Zr-c/w  with  a  phase  determined  by  This  duality  between  w 

and  A  will  be  referred  to  again  below. 
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n 

(b) 


Fig.  10.  Earth-modificotion  foe  tor  E;  (o)  |E(q}|^,  (b)fg(q). 
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One  important  case  in  which  K  is  approximately  constant  over  a  sizeable  region  of  (h,  R,  A,  u) 
space,  and  the  equation  for  E  becomes  particularly  simple,  is  that  in  which  (11  the  polarization 
is  horizontal,  (2)  the  reflecting  surface  consists  of  sea  water,  and  (3)  the  antenna  and  target  both 
have  broad  vertical  patterns.  In  this  case,  for  a  fairly  large  region  of  (h,  R,  A.  w)  space,  one 
has,  approximately,  p  =  -1,  R/(R  +  A)  =  1,  g{I)/g(D)  =  1.  and  s(D,  I)/s(D,  D)  =  1.  Thus  K  =  -1 
and 

E  =  (l-exp(-if^)l2 


|E|^  =  16  sin^(^) 


=  - 


wA 

c 


(13) 


D.  Description  of  Received  Signal  in  Time  Domain 

Let  x(t),  y(t),  and  yj.(t)  be  the  time-domain  counterparts  of  X(w),  Y{u),  and  Yp(u;);  x(t), 
y(t).  and  yj.(t)  be  the  complex  representations  of  x(t>,  y(t),  and  yj.(t);  x(t),  y(t),  and  yj.(t)  be  the 
complex  modulations  of  x(t),  y(t).  and  yj.(t)  around  Thus, 


y(t) 


Y(a»)  exp(iut)  dw 


y(w)  exp(iwt)  doj 


y  (t)  =  ^  ^  y(w)  exp|i(«  -Wq)!!  dw 


y(t)  =  y(t)  expdw^t) 

y(t)  -  2  Re  ly(t)|  -  2  |y(t)|  cos  (w^t  +  <(>y(t)I 


(14) 


quency  span  of  X(w).  Writing  these  constants  as  Z/w^  and  K,  and  making  use  of  Eqs.  (9)  and 


and  similarly  for  x(t)  and  yj.(t).  Assume  nc  v  that  Z(u>)/w  and  K(u»)  are  constant  across  the  fre¬ 
quency  span  of  X( 

(14),  one  obtains 

y(t)  = 


^  !x(t  -  ^)  +  x(t  -  ^  +  2K  5c(t  -  ^  -  ^)1 

C  V  V.  C  C 


w  R 
o 


z  /  ^“^0*^X1-,,  2R,  ,  „2  2R  2A,  /  “*‘^o^\ 

y(t)  =  — ^  exp^-  — -  — )  +  x(t  -  —  -  -^)  exp  (- 


w  R 
o 


(15) 


The  equivalent  free-spare  equations  (i.e..  the  equations  for  yj.(t)  and  yj.(t)I  can  be  obtained  from 
these  equations  by  setting  K  -  0.  If  A  is  sufficiently  small  to  allow  one  to  ignore  its  effect  on 
the  modulation  function  x(t)  (i.e.,  exp(  -iwA/c)  is  constant  across  the  frequency  span  of  X(w)|, 
then  the  equation  for  y(t)  can  be  simplified  to 


y(t)  -  jc(t  -  ^)  F(w^)  E(u;^) 


(16) 
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The  complex  time-domain  eaulvalents  of  the  eouations  Y(w>  =  Y^{w)  E/  '  and  Y..(a;)  = 

r  r 

X(w)  F(w)  of  Sec.  I-C  are  given  by  the  convolution  products 

y(t)  =  yj{t)  .  e(t)  y^it)  =  x(t)  •  f<t)  (17) 

where  f(l)  and  e{t)  are  the  complex  representations  of  the  inverse  transforms  of  F(w)  and  E(w). 
Insofar  as  their  effects  on  y(t)  and  y^it)  are  concerned,  f(t)  and  e(t)  can  be  approximated  by 

f(t)  =  — ^  «(t  - 

e(t)  =  «(t)  +  «(t  -  ^)  +  2K  6(t  -ft  (18) 

where  6{t)  is  the  Dirac  function. 

Letting  A-(t)  and  A-(t)  denote  the  complex  autocorrelation  functions 
x  y 


A-(t)  =  ^  x’*'(t)  x(t  +  t)  dt  =  ^  ^  |X(w}|"  exp(iajT)  dw 
Ay(T)  =  ^  y*(t)  y(t  +  t)  dt  =  ^  |  Y(ti))|^  expdwx)  dw 


one  can  show  that 


(19) 


+  |K|^  +  4!k|^)  A-(t)  +  A-(t  -  •— )  +  A-(t  + 

y  X  X  c  X  c 

o 

+  2K(1  +  |Kl^)  A_(t- |)  +  2K^(1  +  lK|^)  A3j(t+  |)|  .  (20) 

The  autocorrelation  function  Ay(T)  is  the  complex  time-domain  representation  of  the  output  of  a 
filter  which  is  matched  to  the  earth-modified  input  signal.  The  autocorrelation  functions  A^fr) 
and  A^(t)  (defined  in  a  manner  similar  to  A_(T)j  are  related  to  A_(t)  by  the  formulas 

y  y  y 

A_(t)  A~ir)  expdw^T)  A^It)  =  2  Re  [A_(t))  (21) 


E.  Relations  Between  Free-Space  Signal  and  Earth-Modification  Factor 

Let  p  denote  the  bandwidth  of  Ypfw)  and  T  the  duration  of  y^ft).  To  the  extent  that  the  tar¬ 
get  is  a  simple  point  target  (i.e.,  the  impulse  response  of  the  target  is  an  impulse),  P  and  T 
will,  in  most  practical  systems,  be  determined  primarily  by  the  transmitted  signal  X.  Assum¬ 
ing  that  K(w)  is  slowly  varying  across  P,  one  sees  that  E(w)  will  be  locally  periodic  in  w  of 
period  2irc/^  and  the  number  of  cycles  in  P  will  be  pi^/Zitc.  In  the  event  that  Ap/2irc  «  1.  E(w) 
will  be  approximately  constant  across  Yp(u.’)  and  the  effec>.  of  the  earth  will  be  merely  to  multiply 
Yp(w)  by  a  constant.  If  p^/Zrc  »  1,  then  E(w)  will  go  through  many  cycles  within  p  and  the 
effect  of  the  earth  will  be  to  change  the  shape  of  Yp(w).  If  Yp(w)  is  such  that  pT  =s  2ir  (e.g.,  yj.(t) 
is  a  simple  pulsed  sinusoid),  th*.  condition dA/2jrc  »  1  is  also  sufficient  to  ensure  that  E(w)  will 
be  varying  much  more  quickly  than  Yp(w)  within  p.  If,  on  the  other  hand.  /3T  »  2ir  and  Yp(«) 
contains  fine  structure  within  P  (e.g.,  yj.(t)  is  a  pulse  train),  then  in  order  to  ensure  this  result 
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one  must  impose  the  more  stringent  condition  (2jr/T)  A/2irc  =  A/Tc  »  1.  Schematic  illustra¬ 
tions  of  the  two  extreme  cases  j3A/2»c  «  1  and  A/Tc  >>  i  are  shown  in  Fig.  ii(a-b).  [A  preeiac 
picture  of  a  single  cycle  of  E(w)  has  already  been  presented  in  Fig.  10. j  A  slight  change  in  A, 
i.e.,  a  change  of  the  order  of  \/Z  =  »c/cj,  will  have  the  effect  of  shifting  the  phase  of  the  E(w) 
modulation  with  respect  to  Yp(a>).  Whereas  in  the  A/Tc  »  1  case,  such  a  phase  change  will  be 
relatively  unimportant,  in  the  ^A/2»c  «  1  case,  it  may  be  extremely  important.  For  example, 
if  Ik]  is  close  to  unity  and  A  is  such  that  Yp(w)  is  centered  on  a  minimum  of  E(a>)  (often  re¬ 
ferred  to  as  a  "null”),  then  Y(w)  will  be  approximately  zero.  Inasmuch  as  A  -•  0  as  the  target 
approaches  the  horizon,  the  case  ^A/2»c  «  1  will  always  be  important  if  the  radar  is  required 
to  pro*/ide  coverage  at  low  altitudes.  Inasmuch  as  the  maximiun  value  of  A  is  2h,  whether  or 
not  the  case  A/Tc  »  1  can  occur  will  depend  on  the  parameter  2h/Tc. 

The  above  remarks  on  the  relation  of  the  free-spsce  signal  to  the  earth-modification  factor 
have  all  been  stated  In  the  frequency  domain.  Corresponding  remarks  may  be  made  in  the  time 
domain  t  (the  conjugate  of  w)  and  in  the  pathlength-difference  domain  A  (the  dual  of  w).  In  the 
time  domain,  y(t)  will,  in  general,  have  five  different  sections  corresponding  to  five  different 
intervals  on  the  time  axis.  Assuming  fur  oimplicity  that  y^it)  has  a  rectangular  envelope,  one 
can  illustrate  these  sections  for  the  two  extreme  cases  A/Tc  «  1  and  A/Tc  »  1  as  shown  in 
Fig.  12(a-b).  Whethar  or  not  there  will  be  any  overlapping  between  the  pulses  corresponding  to 
the  paths  DD,  II.  and  (DI,  ID)  depends  on  whether  A/Tc  >  1  or  A/Tc  <  1.  (If  1  <  A/Tc  <  2,  then 
only  two  of  the  three  pulses  will  overlap  at  any  one  time.  If  A/Tc  <  i,  then  all  three  will  over¬ 
lap.)  The  value  of  the  envelope  of  y(t)  in  the  overlapping  sections  will  depend  upon  the  frequency 
and  phase  of  the  carrier.  If  j3A/2»c  «  1,  then  y(t)  will  consist  almost  entirely  of  the  result  of 
three  overlapping  pulses  and  the  envelope  will  be  constant  during  this  overlap,  its  value  depv  d- 
ing  upon  the  precise  value  of  A.  If  A/Tc  »  i,  then  the  pulses  will  be  widely  separated  and  the 
precise  value  of  A  will  be  unimportant  ItfiT  »  Zr  and  conditions  are  such  that  0A/2»c  »  1 
and  A/Tc  «  1  simultaneously,  then  the  pulses  will  be  strongly  overlapping  and  the  envelope  will 
be  modulated  in  the  overlapping  section. 


F«g.  11.  Relation  of  lYp(u)  ^^fo  iE(«)l^: 
(o)  9A/2'4c  «  1 ,  (b)  »  1 . 
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Fig.  12.  Five  sections  of  y{t):  (o)  A/Tc  «  1 , 
(b)  A/Tc  »  1 . 
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Turning  attention  to  the  dual  domain  A.  assume  now  that  h.  R,  and  w  are  held  fixeo  and 
consider  the  relation  between  the  functions  Yp(A)  and  E(A),  The  dual  of  the  bandwidth  p  of 
Yp(w)  is  the  width  p'  of  Yp(A),  and  the  dual  of  the  time  duration  T  of  yj{t)  is  the  width  T'  of 
the  Fourier  transform  of  Yp(A).  Inasmuch  as  Yp  depends  on  A  only  through  the  product  g^s, 
if  s  is  assumed  independent  of  A  across  the  width  of  g^(A),  then  P'  gives  the  A-width  of  the 
square  of  the  antenna  function.  Similarly,  if  g^  is  assumed  independent  of  A  across  the  width 
of  s(A),  then  P'  gives  the  A-width  of  the  reflectivity  pattern  of  the  target.  The  dual  of  the  as¬ 
sumption  that  K(w)  is  slowly  varying  across  p  is  the  assumption  that  K(A)  is  slowly  varying 
across  p\  and  the  duals  of  the  parameters  pA/Ztc  and  A/Tc  are  the  parameters  p'u/Zry.  and 
w/T'c.  The  parameter  p'w/Zirc  gives  the  number  of  cycles  of  E(A)  in  P',  and  the  parameter 
w/T'c  determines  the  relation  between  the  period  of  E(A)  and  the  fine  structure  in  Yp{A). 

Assume  now  for  simplicity  that  the  earth  is  flat  and  that  h/H  and  H/R  are  sufficiently  small 

to  approximate  A  by  A  =  Zhy  (see  Sec.I-A).  Ignoring  the  factor  2h,  one  then  observes  that  (1) 

2 

Yp(A)  oc  g  (A)  s(A)  can  be  interpreted  as  the  complex  two-way  antenna  pattern  of  a  hypothetical 
antenna;  (2)  /}'  equals  the  angular  beamwidth  of  this  pattern;  (3)  the  variable  conjugate  to  A 
(the  dual  of  t)  is  the  distance  along  the  aperture  of  this  antenna;  (4)  the  Fourier  transform  of 
Yp(A)  [the  dual  of  y^d))  is  the  current  distribution  along  this  aperture;  (5)  T’  is  the  maximum 
extent  of  this  aperture.  Schematic  illustrations  of  the  two  extreme  cases  p'w/Zirc  «  i  and 
w/T'c  »  1  can  be  obtained  from  Fig.  H  by  merely  replacing  A,  u,  p  and  T  by  w.  A,  /?'  and  T', 
respectively.  The  case  p'w/Zitc  «  1  corresponds  to  a  situation  in  which  the  antenna  gain  pattern 
or  the  target  reflectivity  pattern  is  highly  directional  compared  to  the  E{A)  modulation,  and  the 
case  w/T'c  »  i  corresponds  to  a  situation  in  which  both  of  these  patterns  are  broad  and  smooth 
compared  to  the  E(A)  modulation. 

Although,  in  general,  the  e,xistence  of  this  duality  between  w  and  A  in  E  can  be  very  useful, 

.n  nclually  applying  this  duality  to  concrete  problems,  it  should  be  remembered  that  w  and  A 
■  '  as  duals  in  E  only  to  the  extent  that  they  occur  as  duals  in  K.  In  the  above  discussion,  this 
uu.iiity  in  K  has  been  ensured  by  assuming  that  K  is  slowly  varying  over  the  relevant  regions  of 
u)  and  A.  In  a  particular  application,  however,  the  assumption  that  K  is  slowly  varying  in  w  may 
be  valid,  whereas  the  assumption  that  K  is  slowly  varying  in  A  may  not  be  valid,  and  conversely. 

F.  Modifications  in  Received  Signal  Due  to  Motion  of  Antenna  and/or  Target 

In  the  above  discussion,  it  was  assumed  that  the  antenna  and  target  were  stationary.  -Atten¬ 
tion  will  now  be  given  to  some  of  the  effects  introduced  by  relaxing  this  assumption.  Assume  first 
that  (1)  the  orientations  of  the  antenna  and  target  are  fixed;  (2)  the  translational  velocities  of  the 
antenna  and  target  are  small;  (3)  A  is  small.  [A  more  precise  statement  of  assumptions  (2)  and 
(3)  will  be  given  in  Sec.  I-I.  J  Let 

R(t)  =  R  +  cR't  ,  where  R’  = 

A{t)  =  A  +  cA't  ,  where  A'  = 

Y(w,t)  =  component  of  received  signal  corresponding 
to  transmitted  component  XCw)  exp(iwt) 

Yp(w.t)  =  free-space  component  of  Y{w,t)  , 

F(a),t)  -  Yp(ai,  t)/X(w)  =  free-space  factor 

E((i),t)  =  Y(w,  t)/Yp(w,  t)  =  earth-modification  factor. 
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It  can  then  be  shown  (see  Sec.  l-I)  that 

Y(w,  t)  -  X(u,')  F(«.  t)  E(w.  t)  =  Yp(w.  t)  E{w.  t) 

where 


F(w.t)  =  exp(iw  ((1  -2R')  t  -(1  -R’)  ^l} 

wR  ^ 


E(w,t)  =  (l  +  K(«)  exp(-ia)(A't  + 

One  sees  that  2a;R'  is  the  Doppler  shift  along  DD,  2a)(R'  +  A’)  is  the  Doppler  shift  along  11. 
and  a)(2R'  +  A')  i.s  the  Doppler  shift  along  DI  and  ID.  The  term  u)A'  is  the  Doppler  beat  between 
D  and  1. 

Assuming  that  the  Doppler  shifts  can  be  ignored  in  the  evaluaiion  of  Z(u;)/^‘  and  K(j;).  one 
finds  that  ihe  complex  spectrum  corresponding  to  Y(w.t)  is  giver.'  by 


Y(aj) 


^  {X(rzw> 

wR  “ 


(1  -  R')  (2R/c) 
— — 


+  K‘(w)  X( 


cv 


1  -  2R'  -  2A 


-,)  exp  (— iw 


(1  -  R’)  (2R/c)  -r  2A/c 
1  -  2R’  -  2A’ 


+  2K(w)  X(^--^-^,)  exp(-iu,« 


(1  -  R’)  (2R/c)  +  A/c.-, 
1  -  2R’  -  A' 


(23) 


The  first  term  in  this  equation  represents  the  signal  due  to  DD  {i.e.,  the  free-spate  signal  Yj.4u;)|. 
the  second  the  signal  due  to  II,  and  the  third  the  signal  due  to  DI  and  ID.  As  in  the  static  case, 
Y(aj)  can  be  factored  into  the  products  Y(w)  =  X(a))  F(u))  E(u))  =  Yp.(u;)  E(a))  by  merely  defining  the 
free-space  factor  F  and  the  earth-modification  factor  E  by  the  ratios  F  =  Yp,/X  and  E  -  Y/Yj... 


(b) 

Fig.  13.  Effects  of  translation.  Reiotions  omong  the  energy  spectra  for  components 

DD,  D!,  ID,  ond  II:  (o)  u  A'/P  »  1 ,  (b)  u  A*T/2«  «  1. 

o  o 

Letting  denote  the  center  frequency  of  Vp(w).  one  can  illustrate  ih,.  two  extreme  ca.ses 
w^A’/^  »  1  and  a'^A'/(2ir/T)  =  w^A'T/2ir  «  1  as  shown  in  Fig.  13(a-b}.  When  »  1.  <h» 

signals  along  the  three  paths  DD.  II,  and  (DI,  ID)  will  Oe  separated  in  frequency  and  the  energy 
spectrum  will  be  given  by 
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$ 


Y{w)!^  .  i\X{j^^,)f-  4  \  K{w)f  !X(i-32#:r2^)l^ 


a)‘'R 
+  -li  K//..^  1^  i  VI 


b) 


_u2, 


tOA\ 


1  -  2R‘  -  A*"  * 

When  u^^'T/Zw  «  1.  the  Doppler  beat  will  be  small  even  with  respect  to  the  fine  structure 

in  Yp{w)  and  can  be  ignored,  in  ibis  case,  letting  q{u)  -  -  cl)A/c(1  -  2R')..  one  obtains 

Y(w)  =  YpU‘>)  E(w) 


where 


Yp(a,5  -  ^  exp(-i« 

<uR 

Eiw)  =  {1  +  K(ai)  expi-io) 


and 


where 


(26) 


!Y(w)P=  !Yp(u,)|^  |E(a.)i' 


Ci^  r\ 

iE(w)|^  =  (1  4  1K(w)1^  4  zl  Kiw)!  cosq(w)}^ 

In  this  case,  the  effect  of  the  earth  on  the  free-space  equation  is  identical  to  the  .  atic  case  ex¬ 
cept  for  the  factor  1  -  2R'  occurring  in  the  function  exp{-ia.'(A/c)/(l  -  2R')1.  Thus,  the  pre¬ 
vious  discussion  of  E  applie.A  without  change  except  for  the  contraction  factor  1  -  2R'.  In  many 
cases,  the  effect  of  this  coniraction  factor  on  E  can  be  ignored.  Assuming  that  Z{u))/u)  and 
K(w)  are  constant  across  X{a))  (and  approximating  1  -  2R'  by  1  in  the  amplitude),  one  finds  that 
the  equation  for  the  complex  time  function  y(t;  corresponding  to  Eq.  (23)  for  y(a>)  is  given  by 

y(t)  =  j  {x  ((1  -  2R')  t  -  (1  -  R')  X  [(1  -  2R'  -  2A')  t  -  (1  -  R’)  ^ 

w  R 
o 


+  2Kx  [(1  -2R'  -A')t-(i  -R‘)  ^  -  fl) 


(27) 


Relaxing  the  stationarity  assumptions  still  further,  assume  now  that  the  antenna  is  scanning 
and/or  the  orientation  of  the  target  is  changing.  If  these  motions  are  slow  with  respec;  lo  the 
round-trip  time  2(R  4  A)/c  for  II,  the  complex  time  function  y(t)  will  be  given  by  the  same  for¬ 
mula  as  that  shown  in  Eq.  (27)  except  that  K  =  K(a>  J  and  Z  =  mu.«t  now  be  replaced  by 

K  =  K(u)  ,  1)  and  t)  in  order  to  account  for  the  time  variations  in  ihe  antenna  factors  g(D,  w^) 

and  g(I,  Wq)  and  the  target  factors  s(D,  D,  w^)  and  s(D.  I,  w^).  The  corresponding  spectrum  Y(w) 
can  then  be  obtained  from  Eq.  (23)  by  replacing  Ziw)  and  K(w)  by  the  h'ourier  transforms  of 
Z(w  ,  t)  and  K(w^,  t),  and  by  cepl.acing  the  multiplications  of  X{w),  Z(a)),  and  K(w)  by  convolutions. 
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G.  Total  Received  Energy  and  Energy  Coverage  Diagram 

The  energy  w(Y)  in  tne  received  signal  Y  is  given  by 

i2(Y)  ^  \  |Y(w)i^dw  -  \  y^(t)  dt  =  A  (0) 

|Y(w)|^  dw  -  2  lyd)!^  dt  =  2A_{0)  .  (28) 

Although,  in  general,  actually  computing  J2(Y)  will  be  a  rather  complicated  task,  there  are  a 
number  of  important  special  cases  in  which  this  computation  is  particularly  simple  and  th«*  re¬ 
sults  are  both  praciually  significant  and  intuitively  revealing.  The  formulas  for  filY)  in  these 
special  cases  are  given  l>y 

n(Y)  n(Yj..)  !i:{q)i^  -  nivj  (i  +  !k|^  -t  2!k|  cosq)^ 

fi(Y)  =  nlYp)  aveiE(q)i^  =  n(Yp)  (1  +  IkI^*  +  4|K|^)  (29) 

where  q  ~  a;A/c(l  -  2H')  and  !  I-(q)|  has  been  described  in  detail  in  .Sec.  l-C.  In  both  of 

these  equations,  it  is  assumed  that  K(ai)  is  constant  across  ,3  The  first  of  these  equations  is 
valid  provided  that  both  of  the  following  conditions  are  satisfied:  (1)  the  Doppler  beat  oIqA'  can 
be  ignored  (i.e.,  «  1):  (2)  Klw)  is  constant  across  d  (i.e.,  DA'Zn  «  1).  The  second 

of  these  equations  's  valid  provided  either  one  of  the  following  conditions  is  satisfied:  (3)  the 
Doppler  beat  is  large  (i.i*..  >">  1);  (4)  K(a))  varies  much  more  quickly  than  Yj.(u;) 

across  /j  (i.e..  A/Tc  »  1).  The  second  equation  can  also  be  interpreted  as  giving  the  average 
energy  in  Y  for  a  situation  in  which  conditions  (1)  and  (2)  are  satisfied  and  in  which  the  A- 
coordinate  of  the  target  is  defined  probabilistically  with  all  A’s  within  a  given  cycle  of  E  having 
the  same  probability 

If  all  the  parameters  in  QiY)  are  held  constant  except  the  location  of  the  target,  the  locus 
of  points  which  satisfy  the  equation  n(Y)  =  C  (t  a  cons.ant)  defines  a  surface  in  space  referred 
to  as  a  coverage  contour.  3y  definition,  this  surface  has  the  property  that  whenever  a  target 
with  the  specified  cross  section  is  located  at  any  point  on  this  surface,  the  returned  energy  will 
b«'  equal  to  the  constant  C.  Although  these  contours  contain  no  new  information  and  are  merely 
t ransformaiions  of  the  energy  equation  from  which  they  are  derived,  insofar  as  the  radar's  per¬ 
formance  is  a  function  of  the  received  energy,  they  arc  useful  in  helping  one  to  visualize  the 
regions  of  .“ipace  in  which  a  specified  performance  can  be  expected  to  occur.  An  illustration  of 
the  effect  of  the  earth  on  a  vertical  slice  of  this  contour  is  shown  in  Fig.  14(a-b).  The  curves  in 
Fig.  14  have  been  computed  assuming  that  (1)  the  antenna  and  target  are  isotropic  (i.e.,  |g|  and 
Isj  are  independent  of  angle);  (2)  the  earth  is  flat;  (3)  F  -  -1;  (4)  2irc/u;  (1  —  2R')  =  X  =  10cm; 

•1  A  ^ 

(5)  h  15.000ft;  (6)  n{ V)/12{Y'j.)  R**  -  10~°  (R  measured  in  nautical  miles).  The  total  number  of 
lobes  on  one  side  of  the  antenna  is  given  by  {snaxA)/X  -  2h/X  and.  assuming  that  (h  +  H)/d  is 
small,  the  angular  spacing  between  the  lobes  is  given  by  X/2h.  The  amplitude  of  the  lobing 
modulation  in  this  figure  is  given  by  |  kI  -  H/(R  +  A)  -  (d^  +  (h  -  H)^j  V{d^  +  (h  +  and 

the  phase  of  the  modulation  is  given  by  -  ir.  The  modulation  in  the  received  signal  caused  by 
the  target's  moving  through  these  lobes  is  exactly  the  same  modulation  as  that  which  has  been 
desc  ribed  previously  in  terms  of  the  Doppler  beat  w^A'. 
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Fig.  14.  II lostroMon  of  }h«  eorth's  ftffsct  on  the  coverage  diogrwn:  (a)  detoiis  of  lobe 
structure  (only  five  bottom  lobes  shown),  loci  of  lobe  mcximo,  minimo,  ondoveroge. 

H.  Implications  for  Detection  and  Parameter  Estimation 

In  previous  sections,  an  analysis  has  been  given  of  how  the  signal  received  from  a  target 
is  modified  by  the  presence  of  a  jierfeclly  smooth,  homogeneous  earth.  In  this  section,  a  few 
remarks  will  be  made  on  the  implications  of  this  analysis  for  the  tasks  of  detection  and  param¬ 
eter  estimation.  In  discussing  these  tasks,  it  will  be  assumed  that  the  target  is  a  point  target 
and  reference  will  be  made  to  two  types  of  signal-processing  systems:  one  matched  to  the  free- 

A  A 

space  signal  Y..  (denoted  Y,.)  and  one  matched  to  the  earth -modified  signal  Y  (denoted  Y).  In 

*'  A 

general,  whereas  the  receiver  in  Y.,  w'ill  consist  of  a  filter  bank  of  two  dimensions,  one  .^'or 

**  A 

range  and  one  for  Doppler,  the  filter  bank  in  Y  must  have  additional  dimensions  in  order  to 
account  for  the  factors  K  and  A,  as  well  as  the  Doppler  beat  w  A'.  The  application  of  *o 

A  A  A  ®  A  ^ 

Yp  and  to  Y.  and  of  Y  to  Y,  will  be  denoted  by  Yp  •  Yp,  Yp  •  Y,  and  Y  •  Y,  respectively. 
The  extent  of  the  mismatch  Yp  ■  Y.  and  whether  or  not  the  deterioration  in  performance  re¬ 
sulting  from  this  mismatch  v/ill  outweigh  the  additional  complexities  involved  in  the  construction 

A  , 

of  Y,  will  depend  on  the  specific  situation.  In  the  event  that  /?A/27rc  «  1  and  w  A'T/Z*  «  1, 

A  A  o 

the  two  systems  Y  ana  Yp  differ  only  by  a  complex  constant.  Inasmuch  as  the  over-all  ampli¬ 
tude  and  phase  are  seldon'  matched  by  the  receiver  in  practice  a.nyway.  this  difference  is  a 
trivial  one.  If,  on  the  other  hand. either /J A/ 2xc  »  i  or  le  A'T/2a-  »  1  so  that  Y,.  resolves  the 

A  °  t 

components  of  Y,  and  Yp  •  Y  gives  the  appearance  of  three  distinct  targets  (/?A/2sc  »  1  implies 
the  components  will  be  resolved  in  range  and  u;^A'T/2jr  »  i  implies  the  components  will  be  re¬ 
solved  in  frequency),  'hen  the  mismatch  can  be  quite  significant  One  can.  of  course,  attempt  to 

A 

operate  further  on  Yp  Y  by  making  use  of  one's  a  priori  knowledge  about  the  effects  of  the 
earth,  but  these  further  operations  (whether  they  are  automated  or  take  place  in  the  human  brtiin) 

A.  n,  A 

conslilule  nothing  more  thon  an  attempt  to  t-xtend  to  Y  and  to  construct  Y  -  V. 

*  A  A 

In  the  event  that  ffA/Zre  «  1  and  u:-^A’T/2ir  «  i  (in  which  case  Y  and  Yp  are  essentially 
the  same),  the  change  in  detection  performance  caused  by  the  earth's  presence  will  be  deter- 
mined  by  the  energy  ratio  n(Y)/fl(yp,)  =(l+iKi  t2|K|  cosq) .  If  the  indirect  ray  is 
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-rat 


much  weaker  than  the  direct  r&y  ao  that  one  can  make  the  approximation  K  =  0,  then  this  ratio 


-ill  W**  •••-.la*-  -t«  — 

f*t*  <^/waaaA»»?.«A^  *04  a*J.  t|-» 


ir  the  hiuireci  ray  is  ox  sufficient  strength  compared  to 


thd  direct  ray  to  make  iKj  >  E  {e.g-,  the  antenna  gain  Is  greater  along  I),  then  this  ratio  will  be 
greater  ih;  n  unity  and  the  detection  performance  will  be  improved  for  all  vtaues  of  q.  If  0  < 

(k!  <  E,  then  whether  or  not  this  ratio  is  greater  than  unity  will  depend  on  q.  In  this  latter 
case,  -whether  oi‘  not  the  average  detection  performance  will  be  improved  will  depend  upon  the 
distribution  of  values  of  q  {either  probabilistic  or  time  varying),  and  upon  which  portion  of  the 
probability-of-detectiors  {P  .J-vs-signal-to-rtoise-ratio  {S/N)  curve  is  being  used.  If  the  values 
of  q  can  be  aaaumeU  to  occur  uniformly  so  that  ave  jEr  =  i  +  !kP  +  4lK‘i,  and  the  S/N  ratios 
for  these  values  occur  in  a  region  of  the  P^(S/K)  curve  where  the  curvature  is  non-negative  so 
that  ave  P^{S/N)  >  P^{aveS/N).  then  the  average  detection  performance  will  be  improved  for  all 
values  of  !  Kj  >  C.  If,  on  the  other  hand,  either  the  target  tries  to  hide  in  a  null  of  the  lobe  pat¬ 
tern  by  flying  a  path  of  appropriate  A,  or  the  S/N  ratios  occur  in  a  region  where  the  curvature 
ai  P^iS/N)  is  negative  sc  that  ave  P^{S/N)  <  P^(aveS/N),  then  the  detection  performance  may 
deteriorate.  {Note  that  in  the  case  of  horizontal  polarization  and  a  surface  consisting  of  smooth 
sea  water,  the  target  can  remain  in  a  null  by  merely  flying  close  to  the  surface.)  Inasmuch  as 
most  P^{S/N)  curves  have  positive  curvature  at  low  S/N  ratios  and  negative  curvature  at  high 
S/N  ratios,  the  c?mdition  ave  P^(S/N)  >  P^{aveS/N)  is  likely  to  be  relevant  to  early  warn,i.g  and 
the  condition  avoP^{S/N)  <  F^{aveS/N)  to  high  quality  tracking.  In  order  to  defend  against  this 
possible  deterioration  in  detection  performance,  one  can  make  use  of  one  or  more  of  the  follow¬ 
ing  devices:  (1)  a  high  center  frequency  {which  reduces  the  size  of  all  the  nulls  and  makes  it 
harder  for  the  target  to  hide);  iZi  a  variable  center  frequency  (which  permits  the  radar  "operator" 
to  try  to  track  the  target  in  frequency  and  force  it  to  remain  in  a  ntaximum);  (3)  a  broad-band 
signal  (which  averages  gut  all  but  the  lowermost  lobes);  (4)  a  shift  in  polarization  (which  shifts 
the  phase  of  the  lobe  pattern).  Which  of  these  techniques  will  be  the  most  useful  will  depend 
upon  the  target's  tactics  and  upon  the  specific  function  the  radar  is  intended  to  fulfill. 

In  tho  event  that  p^/Ztc  »  1  or  w^A'T/Ze  »  1  and  the  signals  alcmg  the  three  paths  DD,  II, 
and  (01,  ID'  are  resolved  by  Yp,  it  is  clear  that  the  detection  performance  dp  for  the  three 
cases  Y  Yp  •  Y.  and  Yp  *  Yp  will  satisfy  the  inequality  dp(Y  •  Y)  >  dp(Yp  •  Y)  dp{Yp  •  Yp) 
independent  of  the  values  of  IkI  and  w  and  independent  of  the  chojce  of  flight  path.  The  relation 

A  A 

dp{Y  •  Y)  ^  dp(Y„  •  Y)  follows  from  the  definition  of  the  matched  processing  system,  and  the  re- 

A  A 

lation  dp{Yp  •  Y)  ^  dp{Yp  •  Yp)  follows  from  the  fact  that  Yp  is  a  component  of  Y  which  is 
resolved  from  the  other  components,  (in  order  to  make  precise  statements  al>out  the  detection 

A  A 

performance,  one  must,  of  course,  specify  the  signal-processing  systems  Yp  and  Y  in  detail 
and  compute  probability-of-detection-vs-probability-of-false-alarm  curves.) 

In  addition  to  affecting  the  radar's  ability  to  detect  a  target  (i.e.,  to  vietermine  its  existence), 
the  presence  of  the  earth  also  aifects  the  radar's  ability  to  estimate  various  characteristics  of 
the  target.  According  to  Eq.  (23),  the  parameters  that  enter  into  the  earth -modified  signal  Y 
which  not  occur  in  the  free-space  signal  Yp  (and  estimates  of  which  would  occur  in  the  out¬ 
put  of  Y)  are  the  parameters  A,  A',  and  K.  Looking  at  Eq.  (9)  for  K,  one  observes  that  if  the 
location  of  the  target  and  the  reflective  properties  ot  the  surface  are  known,  tme  can  use  Y  to 
estimate  the  variable  8(D,  I)/s(D,  D)  and  thereby  obtain  information  or  the  vertical  scattering 
pattern  of  the  target.  In  most  practical  cases,  however,  this  estimate  is  likely  to  be  extremely 
noisy  (particularly  when  a  rough  reflecting  surface  is  conside'''ed).  Attention  will  therefore  be 
restricted  to  the  parameters  A  and  A'. 


21 


m 


<9 


With  the  ordinary  search  radar  visualized  in  free  space,  the  only  variables  describing  the 
target's  position  and  velocity  within  the  vertical  plane  that  can  be  measured  with  any  degree  of 
accuracy  and  reliability  are  the  ranee  H  and  the  range  rate  R'.  In  order  tc  deterinine  fuf-iher 
coordinates  in  this  plane,  one  must  employ  an  antenna  pattern  with  a  high  degree  of  vertical 
resolution  (i.e.,  a  high  frequency  or  a  large  vertical  aperture).  This  constraint  has  frequently 
proven  to  be  sufficiently  incompatible  with  the  search  function  to  necessitate  the  construction  of 
a  completely  separate  radar  (usually  referred  to  as  a  "height -finder")  whose  specific  function 
is  to  determine  these  further  coordinates.  If  the  earth's  surface  is  taken  into  account,  however, 
and  the  radar's  antenna  pattern  and  the  reflecting  properties  of  the  surface  are  such  that  one  is 
assured  of  having  both  a  direct  ray  and  an  indirect  ray,  then,  as  mentioned  above,  the  returned 
signal  Y  can  be  used  not  only  to  estimate  R  and  R',  but  also  A  and  A'.  If  the  radar's  own 
altitude  h  and  rate  of  climb  h'  are  known,  R  and  A  are  sufficient  to  determine  the  target's 
position  in  the  vertical  plane,  and  R'  and  A*  are  sufficient  to  determine  its  velocity  vector  in 
this  plane.  The  basic  idea  of  these  techniques  is  to  incorporate  the  earth's  surface  as  part  of 
the  antenna  and  thereby  eliminate  the  need  for  an  additional  radar. 

in  addition  to  requiring  that  both  rays  be  sufficiently  strong,  the  ability  to  make  high  quality 
measurements  of  A  and  A’  (like  R  and  R’)  will  depend  on  having  a  sufficiently  large  bandwidth 
P,  a  sufficiently  long  time  duration  T,  and  a  sufficiently  refined  data -processing  system.  In 
the  event  that  R/2h  is  large  and  the  information  on  the  target's  location  is  required  in  the  co¬ 
ordinates  R  and  H  rather  than  R  and  A,  the  measurement  of  A  will  have  to  be  extremely 
accurate  to  produce  an  accurate  measurement  of  H  since  dll/dA  will  then  also  be  large.  [For 
a  flat  earth,  H  =  (2AR  +  A^)/4h.]  Note  also  that  in  order  to  transform  the  velocity  coordinates 
R'  and  A'  to  R'  and  H',  one  needs  to  know  both  A  and  R,  and  therefore,  in  order  to  produce 
R'  and  H',  one  needs  not  only  a  long  time  duration  T,  but  also  a  Large  bandwidth  p.  One 
makeshift  technique  for  determining  H  which  has  been  used  on  occasion  when  P  has  been  too 
small  to  measure  A  directly  and  some  outside  intelligence  is  available  on  H',  is  to  measure 
R,  R',  and  A'.  This  technique  of  estimating  H  from  a  measurement  of  A'  is  referred  to  as 
"lobe  counting."  Its  weakness,  as  compared  to  a  direct  measurement  of  A,  is  Us  dependence 
on  outside  information  concerning  H'  (information  which,  under  most  conditions,  must  be  very 
precise)  and  the  long  time  duration  required  to  measure  A'. 

So  far.  it  has  been  assumed  that  only  one  target  is  present.  In  a  realistic  environnrient,  tne 
number  of  targets  is  unknown.  In  this  case,  the  existence  of  the  earth-reflected  ray  can  lead  to 
a  variety  of  ambiguities.  For  example,  if  three  distinct  signals  are  received  separated  by  equal 
intervals  (either  in  Doppler  or  range),  these  signals  can  be  interpreted  as  arising  either  from  a 
single  target  with  the  three  paths  DD,  II,  and  (DI,  ID)  or  from  three  targets,  each  with  only  a 
single  path.  Similarly,  if  six  equally  spaced  signals  S^,  . . . .  are  received,  these  signals  can 
be  identified  with  targets  T.  in  the  fashion  (T^  --  (S^.S^.S^),  T^  =  (S^,  S^,  Sj^)|:  [T^  =  (S^.S^.S^), 
T2  =  (S^.S^.S^)!;  (Tj  =  Sj,  T^  =  (Sj,S^,S-).  T^  =  -S^j;  and  so  foith. 

In  general,  the  earth's  presence  not  only  introduces  new  target  parameters  to  be  estimated, 
but  also  affects  the  estimates  of  the  old  free-space  parameters.  In  order  to  make  precise  state¬ 
ments  about  the  influence  of  the  earth  on  the  radar's  parameter-estimation  ability,  not  only  must 
the  processing  system  be  specified  in  detail,  but  one  must  select  criteria  for  evaluating  the 
estimates. 
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I.  Generalised  Eq^tion  and  Genera.' ized  Approach 

to  sifects  oi  Trajismtioa-Motitm 

Before  tunning  to  the  coropUcotions  intreduccu  by  surfHCc  roughness,  it  is  imponam  lo  note 
that  the  previous  discussion  is  extremely  aimplified  even  for  an  earth  which  ia  perfectly  smooth. 
In  this  sectitm,  the  previous  results  will  be  generalized  in  two  directions.  First,  the  static  equa¬ 
tions  of  Sec.  I-C  will  be  generalized  by  omittutg  a.ssumptions  (1).  (2),  and  (5)  of  that  section. 
Second,  a  teebniqui;  will  be  presented  for  incorporating  the  effects  of  translations  of  the  antenna 
or  of  target  which  can  ?>e  used  to  generalize  the  results  of  Sec.  I-F.  In  both  cases,  the  re¬ 
sults  stated  previously  will  be  derived  as  special  cases. 

1,  Generalization  of  Section  I-C 

Assume  now  that  assumptions  (1),  (2).  and  (5)  of  Sec.  I-C  are  omitted.  Let  p  denote  the 
polarization  h  or  v  and  t  the  path  D  or  1.  Decomposing  the  electric  vectors  along  D  and  I 
into  the  sums  of  their  h  and  v  components,  one  can  write  the  total  received  signal  as  the  sum 
of  sixteen  terms,  each  of  which  gives  the  received  signal  for  a  particular  path-polarization 
combination  (p^,  p2,  t^)  (where  "(Pj,  P2.  < *2^"  means  "p^  component  out  along  and 

component  back  along  t^"].  Denoting  the  transmitting  and  receiving  antenna  functions  by  g^  and 
g^.  and  inserting  the  polarization  arguments  into  g^.  g^.  s.  rmd  p  (but  omitting  the  argument  w), 
define 

Z{pj.  P2»  *2^  "  2  ®r^P2'  *2^  Pa*  *1'  *2^ 

(tip,  f)  =  1  for  I  =  D  pip,  t)  =  p(p)  for  f  =  I 

L(«  =  length  of  i  .  (30) 

Continuing  to  omit  the  argument  w  wherever  possible,  one  then  has,  as  the  generalization  of 
Eq.  (9). 


w 


Pi'  P2  x) 

I'.f^^in  (D.I) 


^<Pl'P2**l' V  P^Pi*  V'  P^P2'^2^  exp{-i(w/c)  (L(l^)  +  LCf^))} 

L{fj)  LH^) 


where 


=  XFE  =  YpE 


F  = 


(31) 


Pi'  P2  ^ 


Z(Pl.  P2.  D,  D) 


E  = 


R^  exp(-2i«A/c)  2  Z(p.,  p,.  1, 1)  p(p.)  p(p-) 

_ P^.Pzinfh.v) 

(R  +  A)^  2  Z{p.p2.D,D} 

p^.  P2  in  (h.  v} 

R  6x0  (“ ici/^/c)  £  ^2*  ^1*  ^2^  ^^^2* 

Pj.P2in(h.v} 

*1^2^*^  (D,  1} 

(R  +  A)  2  Z(p^.  P2.  D.  D) 

Pj,  P2  in  (h.  v) 
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Since  the  tranomliiiiig  and  receiving  .antennas  are  no  longer  assumed  identical,  the  reciprocity 
theorem  is  no  longer  applicable  and  one  can  no  longer  equate  the  signal  received  over  the  pato 
DI  with  the  signal  received  over  the  path  ID  (i.e.,  Z(p^,  is  not  necessarily  equal  to 

Z(py,  p^,  fj)|.  Also,  since  reflection  from  the  earth's  surface  {even  when  this  surface  is  per¬ 

fectly  smooth)  will  depolarize  any  wave  which  is  not  pure-horizontal  or  pure-vertical  [due  to  the 
differences  between  p(h)  and  p(v)|,  the  target  scattering  length  s  cannot  be  defined  merely  in 
terms  of  the  polarization  of  the  antenna  as  was  done  in  Sec.  I-C  and  as  is  usually  done  in  the 
discussion  of  the  free-space  case. 

If  one  now  restores  assumptions  (1)  and  (2)  of  Sec.  I-C,  the  equation  for  F  becomes  the 
same  as  in  Sec.  I-C  and  (defining  K  as  in  Sec.  I-C  and  dropping  the  polarization  notation)  the 
equation  for  E  becomes 


E  =  (1  +  K  exp(-iwA/c)|2 


exp(-2ia>A/c)^®^ 


1  + 


,I)s(D.D)  , 

2 - ^ 

s  (D.  1) _ . 


[1  -t  K  exp(-ia,'A/c))' 


(32) 


The  role  of  assumption  (5)  in  Sec.  I-C  was  to  reduce  the  second  factor  in  this  equation  to  unity. 
[Note  here  the  asymmetry  in  s  and  g  with  lespect  to  the  assumptions  required  for  writing  E 

as  the  square  of  1  +  K  exp(-iwA/c).  The  reason  for  this  asymmetry  is  that  whereas  one  always 

2  2  2  2 
has  g  (I)  g  (D)  =  [g(D)  g(I))  ,  the  equality  s(1, 1)  s(D,  D)  =  s  (D,  I)  demands  a  special  assumption.) 


2.  Algorithm  for  Generalizing  Section  I-F 

Attention  will  now  be  turned  to  presenting  a  technique  for  generalizing  the  equations  of 
Sec.  1-F  describing  the  effects  of  antenna  or  target  translations.  T.his  technique  does  not  take 
account  of  the  effects  of  the  motion  on  the  various  reflection  processes  involved  or  of  the  changes 
in  the  angles  of  orientation,  but  is  merely  a  nonrelativistic  algorithm  for  taking  account  of  the 
time-varying  nature  of  the  travel  times  of  the  signals. 

Ignoring  all  attenuation  and  distortion,  one  can  write  the  received  signal  for  a  given  trans¬ 
mitted  comjKtnent  expdwt)  and  a  given  path  as  exp{iw{t  -  T(ty  where 

is  the  delay  suffered  by  the  signal  which  travels  along  (f^,  ^2)  arrives  back  at  the  radar  at 
time  t.  Inasmuch  as  r(fj,  will  not,  in  general,  be  linear  in  t,  the  effect  of  t's  dependence 
on  t  cannot  be  described  merely  in  terms  of  a  Doppler  shift  of  w.  Also,  one  cannot  expect 
t(<^.  ^2**^  symmetric  in  and  (since  the  path  (D,  1)  will  be  geometrically  distinct  from 

the  path  (I,  D)).  or  t^;  t)  to  be  a  function  merely  of  the  relative  motions  of  the  antenna  and 
target  a'ong  D  and  I  (since  there  are  three  bodies  involved  rather  than  only  two).  In  order  to 
compute  t(Ij,  t),  consider  a  single  photon  that  travels  over  the  path  (f^,  t^)  and  arrives  back 
at  the  radar  at  time  t,  and  let 

t  =  time  at  which  photon  is  transmitted 
t^  =  time  at  which  photon  hits  target 

A(l,;t .)  =  ti  -  tjj  =  time  it  takes  photon  to  go  out  along  a.s  a  function 
of  the  time  at  which  it  hits  the  target 

B{<2;t)  =  t  -  tj  =  time  it  takes  photon  to  get  back  along  f2  ss  a  function 
of  the  time  at  which  it  arrives  back  at  the  radar 

L(f:x.  y)  =  distance  along  f  between  radar  at  time  x  and  target  at  time  y. 

Thus 
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A{t^;  t  j)  =  -  t  j)  =  ~  L  -  A(/j;  t^),  t ^  I 

J  M#2;t.tj)  =  I  L  ‘  -B{#2;t)| 

t{^|,  ^2'  ^  ~  *o  ”  ^  1^  ^  ^2^ 

=  A  (l^;t  -  B(#2;tH  +  B{#2;t)  (33) 

The  first  two  equations  are  functional  equations  for  A(#^;tj)  and  8(^2;  0  in  terms  of  the  arbitrary 
distance  functions  L(<^:x.  y)  and  L(/2,  x,  y).  and  their  solutions  in  terms  of  these  functions  can  be 
approximated  by  iteration.  F'or  example,  for  A{/j;tj)  one  obtains  the  sequence 

J  !-(<,;•,.<,>  ■ 

=  i  L  1.1,1 

etc. 

Once  A(fj;tj)  and  B(<2;  12)  have  been  determined.  t(#^.  <2:*)  can  be  computed  by  use  of  the  third 
equation. 

In  order  to  see  how  Eq.  (22)  of  Sec.  I-F  is  derived  from  this  general  algorithm  (and  to  make 
assumptions  (2)  and  (3)  of  Sec.  I-F  more  explicit),  let  v(f)  =  velocity  of  radar  along  t  at  t  =  0, 

V(f)  -  velocity  of  target  along  i  at  t  -  0.  and  M<;0)  =  1,(#;0.  0).  Assume  that  (a)  L(I;x.  y)  can  be 
approximated  by  the  first  two  terms  of  its  Taylor  series  expansion  around  (x  -  0,  y  -  0);  (b) 
A(fj;tj)  and  B((2;t)  can  be  appro*;  nated  by  A^^'(f^;  t  j)  and  >ccms  in  v^{f)/c^ 

and  \^{t)/c^  can  be  ignored.  Witi.  these  assumptions.  T(t^,  t^:\)  can  be  approximated  by 

T(f  #2:  t)  =  I  (v«i)  +  V(f  j)  +  v(f2)  +  Va^)}  t  +  ^  ( 1  -  ^  (v(f,)  +  V((j) 

+  V(f2)|}Mf2:0)  +  I  (1  -  I  v(fj)l  L(/j;0)  (34) 

Since,  according  to  this  approxima'ion.  T(t^,f^;i)is  linear  in  t.  the  function  cxpfiw  (t  -  r(fj. 
can  be  rewritten  as  exp(i  (w(#j.  t^)  t  -  0(1^.  f2))^-  where  w(fj.  and  arc  independent  of  t. 

Letting  r'  denote  the  time  derivative  of  t.  one  has 

=  «  |1  -  r'(fj,(2:t)| 

Letting  and  vary  over  D  and  I,  one  obtains  for  the  frequencies  w(fj.  l^}  the  usual  Doppler 
shifted  frequencies 

u.(D.D)  -  u)(l  ~  I  fv(D)  4  V(D)!} 
w(I.I)  -  w(l  -  I  (v(I)  ^  VU)1) 
aj(D.  1)  -  w(I,D)  =  I  |w{U.  D)  t  wd,  D) 

wfi  -  i  |v(D)  +  V(D)  +  v(l)  4  Vd))}  .  (36) 
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Writing  R  =  MD;0)  and  A  =  L{I;  0)  -  L(D;0).  and  assuming  the  validity  of  the  approximations. 


wR  .IvU)  -v<D)|  "  lV(l)  -  V<Di|,  ,  0 

c  ^  c 

one  can  approximate  t  le  phases  <p{iy  by 

<i>{D.D)  =  (i  -  i  fvtD)  +  V{D)1} 

^  (1  _  1  {v(D)  +  V(D))}  +  ^ 
ip{D.  I)  =  D)  =  j  (v>{D.  Di  +  <0(1.  D) 

Writing  R'  for  [v(D)  +  V(D))  /c.  and  ,2k’  for  (v(I)  4  V{I)  -  [v(D)  4  V(D)l)/c,  then  leads  to  liq.  (22) 
for  y(w.  t)  ip  Sec.  I-F. 

n.  GENERAL  REMARKS  ON  INFLUENCE  OF  A  ROUGH  EARTH 

In  Sec.  I,  the  received  signal  was  described  under  the  assumption  that  the  earth's  surface 
const.tmed  a  smooth,  homogeneous  sphere.  For  certain  iyp<  s  of  radars,  terrains,  and  radar 
probJf;  ns,  the  results  based  on  this  assumption  will  constitute  a  useful  approximation.  In  ether 
cases,  however,  these  results  will  be  inappropriate  and,  at  best,  can  only  s«>rve  as  a  starting 
point  for  further  analyses. 

A.  Ba&ic  Effects  of  Surface  Roug^ess 

For  the  radar  designer,  roughness  (either  in  shape  or  in  the  properties  of  the  earth  sur¬ 
face  material)  has  two  effects:  (1)  it  corrupts  the  idealized  picture  of  forw'ard  scattering  pre¬ 
sented  in  Sec,  I  and  thus  further  modifies  the  rharacter  of  the  signal  received  from  the  target; 

(2)  it  spreads  out  the  .source  of  the  backscatiered  signal  (clutter)  and  thus  gives  importance  to 
the  earth's  role  as  a  competing  target.  Insofar  as  the  complexity  of  the  earth's  surface  usually 
requires  that  these  effects  be  treated  statistically,  to  the  extent  that  one  is  only  interested  in 
gathering  information  about  the  target  and  not  about  the  eaith's  surface,  both  effects  must  be  re¬ 
garded  as  potential  sources  of  confusion,  and  both  must  be  evaluated  in  order  to  predict  a  radar's 
performance  The  second  effect,  which  leads  to  a  form  of  interference,  tends  both  to  increase 
the  total  level  of  signal  again.st  wiiich  the  target  signal  must  compete  and  to  introduce  a  new  fluc¬ 
tuation  spectrum  into  this  signal.  7'he  lirst  eifect,  which  may  be  viewed  as  a  statistical  corrup¬ 
tion  of  the  earth-modification  factor  E  discussed  in  Sec.  I.  tends  to  wash  out  the  smooth -surface 
interference  lobes  (replacing  them  by  an  irregular  fine  structure),  and  to  introduce  a  new  fluc¬ 
tuation  spectrum  into  the  target  signal  The  influence  of  the  roughness  on  the  total  received 
energy  from  the  target  will  depend  upon  the  specific  characteristics  and  orientations  of  the  an¬ 
tenna  and  tariret,  as  well  as  upon  the  specific  characteristics  of  the  surface.  Of  the  two  effects, 
the  one  which  has  received  most  attention  by  radar  designers  is  the  interference  effect.  On  the 
whole,  the  corruption  of  the  earth-reflected  target  signal  has  been  regarded  as  an  effect  of  the 
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second  order.  Since  the  clutter  Signal  may  often  be  as  much  as  50  db  above  the  internal  noise 
and  one  may  not  be  able  to  see  the  target  at  all,  corruption  or  no  corruption,  this  attitude  is 
quite  understandable.  Except  for  the  brief  description  in  Sec.  Ill  of  forward  scattering  from  the 
ocean,  this  report  will  be  concerned  exclusively  with  the  clutter  oroblem. 

a  "State  of  the  Art" 

Despite  the  fact  that  a  considerable  amount  of  effort  has  been  spent  in  studying  scattering 
from  the  earth's  surface,  at  present  there  is  no  unified  body  oi  knowledge  which  is  adequate  for 
radar  design.  More  specifically,  one  finds  that  (t)  although  there  is  a  considerable  amount  of 
experimental  data  available  on  certain  aspects  of  the  problem,  many  of  these  data  are  incon¬ 
sistent  and  incomplete;  (2)  although  there  are  many  theoretical  studies  of  scattering  from  rough 
surfaces,  there  is  no  theory  available  which  starts  out  with  a  realistic  descripiion  of  tfie  earth's 
surface  and  obtains  results  which  are  both  consistent  with  the  experimental  da<a  (a  oilficuli  task 
since  the  data  are  themselves  sometimes  inconsistent!)  and  sufficiently  determinate  to  be  of  use 
to  tne  radar  designer.^  The  basic  reasons  for  this  state  of  affairs  appear  to  be  the  following. 
First,  the  surface  being  considered  is  extremely  complex.  In  experimental  work,  this  com¬ 
plexity  usually  results  in  an  inadequate  measurement  program  for  determining  the  surface 
charactertistics  Thus,  a  substantial  amount  of  data  exists  in  which  certain  aspects  of  the 
scattered  field  are  described  in  great  detail,  but  the  characteristics  of  the  surface  which  has 
caused  this  field  are  described  in  no  detail  at  all.  In  theoretical  work,  this  complexity  results 
either  in  equations  which  are  unusable  for  concrete  computations  or  in  the  replacement  of  the 
actual  surface  by  a  model  which  is  extremely  artificial.  In  addition,  this  complexity  see.ms  to 
have  caused  a  great  gap  between  theoretical  and  experimental  work.  Not  only  have  certain 
theories  been  proposed  with  only  token  reference  to  the  data,  but  many  experiments  appear  to 
have  been  designed  without  regard  for  the  results  which  have  been  obtained  theoretically.  Sec¬ 
ond,  the  systems  used  for  m*asuring  the  scattered  signal  are  often  quite  complex.  Thus,  one 
often  finds  experimental  results  in  which  it  appears  that  the  experimenter  has  been  unable  to 
obtain  sufficiant  informati  in  on  the  measuring  device  to  allow  him  to  separate  those  effects 
which  are  due  to  the  scattering  prooerties  of  the  surface  from  those  which  are  due  to  the  meas¬ 
uring  device.  Third,  from  the  scientific  viewpoint,  the  earth-scattering  problem  is  basically 
uninteresting  (i.e.,  its  solution  vill  have  no  implications  for  physical  theory),  and  the  main 
motivation  for  studying  this  pi  oblem  has  been  the  desire  to  improve  radar  and  communication 
pe-formance.  Thus,  very  few  experiments  have  been  designed  with  a  view  toward  understanding 
the  scattering  phenomenon  as  a  whole.  For  e.xample,  practically  no  work  has  been  done  to  ob¬ 
tain  a  picture  cf  the  whole  scattered  field.  Experimenters  have  studied  either  backscattering 
or  forward  .’scattering  (depending  upon  the  specific  application  in  mind),  but  not  both.  Finally, 
the  reader  should  al.so  be  aware  cf  the  fact  that  even  if  there  were  an  adequate  theory  for  de¬ 
scribing  the  scattered  field,  the  radar  designer  would  still  be  faced  with  the  problem  that  there 
is  no  adequate  radar-design  theory  for  minimizing  the  effecis  of  surface  roughness.  This  ques¬ 
tion  will  be  di.scussej  furthci  in  See.  U-D. 


fA  good  'ntrodvctofy  biblttigrophy  of  thaoraticol  pop«n  on  scottaring  from  rough  v;rfoc*s  can  b«  founo  In  the 
paper  by  beckmotinP 
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In  general,  the  types  of  interference  with  which  the  target  signal  will  have  to  compete  are 
(1)  internally  generated  noise  (e.g.,  thermal  noise  in  the  receiver,  (2)  externally  generated 
noise  (e.g.,  other  radars),  and  (3)  clutter  (e.g.,  birds,  storms,  earth).  The  distinguishing 
feature  of  clutter  is  that  it  represents  energy  transmitted  by  the  radar  and  reflected  from  the 
environment,  whereas  the  other  forms  of  interference  represent  energy  generated  outside  the 
radar  transmitter.  Accordingly,  the  nature  of  the  clutter  signal  will,  in  general,  be  markedly 
different  from  the  nature  of  the  other  interfering  signals,  and  the  techniques  that  will  be  useful 
for  its  elimination  will  be  markedly  different.  For  example,  whereas  increasing  the  trans¬ 
mitted  energy  is  a  fundamental  technique  for  overcoming  noise  interference,  it  is  an  utterly 
useless  technique  for  overcoming  clutter  interference:  the  target -to-clutter  ratio  will  remain 
fixed.  On  the  other  hand,  whereas  for  many  forms  of  noise  the  idea  of  resolving  the  target  and 
noise  is  useless  because  the  noise  is  too  extensive  (i.e.,  the  noise  exists  at  all  times  and  all 
frequencies},  in  the  case  of  clutter,  since  the  signal  represents  the  return  from  another  well- 
defined  target  with  definite  boundaries  in  range,  angle,  and  velocity,  it  may  be  possible  to  re¬ 
solve  the  target  and  clutter  and,  subsequently,  to  identify  and  suppress  the  clutter  signal.  There 
are  certain  aspects  of  the  clutter  and  noise  problems  which,  of  course,  are  similar.  For  ex¬ 
ample,  insofar  as  the  target  and  clutter  signals  are  not  resolvable  in  range,  angle,  or  velocity, 
one  may  find  it  useful  to  employ  a  matched  filter  against  clutter  in  much  the  same  way  as  one 
employs  a  matched  filter  against  noise.  In  general,  however,  the  techniques  which  have  evolved 
and  which  are  actually  employed  in  the  two  cases  are  quite  different. 

Regardless  of  the  detailed  shape  and  electromagnetic  characteristics  of  that  portion  of  the 
earth's  surface  over  which  the  radar  is  operating,  there  are  certain  fundamental  properties  of 
the  earth  target  as  compared  with  the  aircraft  target  that  are  valid  for  most  portions  of  the 
earth's  surface  and  which  lead  directly  to  a  large  number  of  techniques  for  the  elimination  of 
earth  clutter.  Briefly,  these  properties  are:  (1)  the  earth  target  is  located  at  a  different  posi¬ 
tion  in  space  than  the  airc'aft  tar/’et;  (2)  the  earth  target  has  a  different  relative  velocity  than 
the  aircraft  targc«.;  (3)  the  earth  target  is  much  larger  than  the  aircraft  target;  (4)  aircraft  tar¬ 
gets  which  occur  above  the  horizon  line  occur  at  different  angles  than  the  earth  target;  (5)  the 
relative  position  and  velocity  of  the  earth  target  are  knowable  beforehand,  whereas  the  rtlavive 
position  and  velocity  of  the  aircraft  target  are  probabilistic;  (6)  the  aircraft  target  is  elevated 
above  the  earth  ana  therefore  may  have  in  image  below  the  earth's  surface,  whereas  the  earth 
(as  a  target)  will  not.  Additional  oroperties  which  may  be  useful  for  distinguishing  between  the 
earth  target  and  the  aircraft  target,  such  as  the  de^>endence  of  the  reflection  characteristics  on 
frequency,  polarization,  and  incidence  angle,  will  depend  on  the  detailed  characteristics  of  these 
targets.  (It  should  also  be  no'cd  that  if  the  altitude  of  the  aircraft  ^s  of  the  same  order  of  magni¬ 
tude  as  the  heights  of  the  irregularities  on  the  surface,  or  if  the  radar  has  sufficient  resolution 
to  resolve  individual  scattering  elements  on  the  surface,  then  certain  of  the  above  statements 
will  no  longer  be  appropriate  and  the  features  that  will  be  useful  for  distinguishing  between  the 
aircraft  and  the  surface  may  depend  even  more  heavily  on  these  detailed  characteristics.) 

Despite  tl.e  fact  that  many  land  surfaces  are  poore.>'  reflectors  than  the  surface  of  the  ocean, 
the  interference  caused  by  land  clutter  is  usually  more  difficult  to  overcome  (or  even  to  describe) 
than  that  caused  by  sea  clutter.  Hot  only  is  the  total  energy  in  the  clutter  signal  at  small  graz¬ 
ing  angles  often  greater  for  land  surfaces  due  to  the  greater  slopes  encountered,  but  also,  the 
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statistical  properties  of  the  signal  are  often  much  more  complicated.  The  basic  reasons  for  this 
added  complexity  are:  (1)  the  electromagnetic  properties  of  land  surfaces  tend  to  be  less  homo¬ 
geneous;  (2)  the  variety  of  surface  irregularities  tends  to  be  greater,  both  in  size  and  in  shape; 
(3)  these  irregularities  tend  to  be  distributed  less  uniformly.  Whereas  in  many  cases  the  sea 
clutter  can  be  treated  as  a  .'simple  random -noise  phenomenon,  the  land  clutter  must  often  be 
treated  as  a  multiple-signal-plus-random-noise  phenomenon.  Some  results  on  the  characteris¬ 
tics  of  land  clutter  can  be  found  in  the  report  by  Katzin,  Wolff  and  Kat..in.^  A  description  of  sea 
clutter  is  given  in  Sec.  III. 

D.  Techniques  for  Reduction  of  Earth  Clutter 

As.suming  that  the  function  of  the  radar  being  designed  is  to  gather  information  on  targets 
other  than  the  earth's  surface,  one  must  regard  the  signal  returned  from  this  surface  as  a  form 
of  clutter  and  de.sign  the  radar  in  such  a  way  that  the  effects  of  this  clutter  will  be  minimized. 

In  principle,  for  any  given  system  of  c<.  nstraints  imposed  by  the  state  of  component  technology 
and  for  any  given  radar  function,  it  should  be  possible  to  determine  the  system  which  will  best 
satisfy  that  function  in  the  presence  of  clutter.  In  praciice,  however,  the  problem  usually  be¬ 
comes  so  complicated  that  the  optimum  system,  is  never  actually  determined.  This  state  of 
affairs  is  due  primarily  to  the  following  facts.  First,  as  already  indicated,  the  statistical 
characteristics  of  the  clutter  signal  depend  not  only  on  the  nature  of  the  .scattering  surface,  but 
also  on  a  large  complex  of  factors  describing  the  behavior  of  the  radar.  Thus,  in  trying  to 
optimize  the  system  with  r^sp'^ci  to  a  given  radar  function,  one  is  faced  w'ith  an  optimization 
problem  in  which  not  ortiy  mu^i  many  variables  be  optimized  over  simultaneously,  but  one  in 
which  each  variable  has  an  effect  on  the  interfering  signal  as  well  as  on  the  target  signal.  Sec¬ 
ond,  as  also  indicated  previously,  one's  knowledge  of  now  the  statistics  of  the  clutter  signal 
vary  with  the  various  radar  parameters  is  still  very  limited.  Finally,  thc.so  stati.stios  will,  in 
general,  be  nonstationary.  Thu.s,  for  the  system  to  be  optimum,  it  will  have  to  be  adaptive;  the 
basic  design  of  the  radar  will  have  to  contain  a  pLan  for  sampling  the  clutter  .signal  and  adjusting 
the  various  radar  parameters  according  to  the  results  of  this  sampling.  In  general,  it  is  clear 
that  the  problem  of  determining  an  optimum  anticlutter  system  is  very  complex  and  that  no  one 
has  yet  determined  such  a  system.  On  the  other  hand,  due  to  the  importann?  cf  the  clutter  prob¬ 
lem,  considerable  effort  has  been  spent  in  developing  specific  anticlutter  techniques  and  in  de¬ 
veloping  optimization  theories  that  are  valid  within  certain  limited  contex*s.  Roughly  speaking, 
these  efforts  can  be  divided  into  two  classes,  according  to  whether  or  not  they  make  use  of  the 
specific  reflective  charactf  "istics  of  the  clutter  source  (i.e.,  the  cluttei  's  dependence  on  fre¬ 
quency,  polariz.ation,  and  grazing  angle).  The  purpose  of  the  remarks  in  this  section  is  to  out¬ 
line  some  of  the  clutter-rejection  schemes  which,  at  lea.st  in  basic  conception,  arc  independent 
of  these  characteristics.  Some  techniques  which  canitalize  on  these  cliaracteristics  for  tlie 
special  case  of  the  ocean  can  be  derived  from  the  description  of  these  characteristics  given  in 
Sec.  III. 


1.  Elevated  Target  Indicators  (.ETI) 

Since  the  target  and  the  earth's  surface  will  always  be  located  at  different  positions  in  space, 
it  is  theoretically  possible  to  eliminate  the  effects  of  clutter  by  designing  a  radar  with  a  high  de¬ 
gree  of  spatial  resolution  and  ignoring  all  signals  which  correspond  to  positions  on  the  surface. 
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Inasmuch  as  the  surface  will  exist  at  all  a^amuths  and  all  ranges  oat  to  the  horizon  range,  aside 
from  requiring  that  the  radar  have  sufficient  range  resolution  to  be  able  to  separate  returns 
ftom  insiJc-  and  oui&iue  the  horizon  range,  the  basic  resolution  requirement  will  be  on  the  alti¬ 
tude  variable  K,  the  height  of  the  target  above  the  surface.  Any  system  which  can  be  designed 
to  filter  out  s  gnals  that  arise  from  scatterers  with  zero  elevation  will  be  effective  in  eliminat¬ 
ing  the  clutter  signal. 

The  most  obvious  way  of  obtaining  such  a  system  is  to  employ  a  transmitted  waveform  with 
a  high  degree  of  range  resolution  and  an  antenna  pattern  with  a  very  narrow  vertical  beam,  and 
to  filter  out  the  zero-height  signals  by  a  simple  gating  process.  A  somewhat  more  elegant  sys¬ 
tem  which  make.s  smaller  demands  on  the  vertical  "perture  of  the  a.ntenna  is  to  replace  the  nar¬ 
row  vertical  beam  by  a  broad  vertical  beam  combined  v/ith  an  adjustable  vertical  null.  Targets 
with  positive  altitudes  can  then  be  detected  by  scanning  the  null  through  the  region  of  elevation 
angles  corresponding  to  the  position  of  the  clutter  and  comparing  the  observed  return  with  the 
position  of  the  null.  If  no  target  of  positive  altitude  is  present  at  range  R,  the  return  fr  om 
range  R  will  be  zero  when  the  null  is  pointing  in  the  direction  of  range  R  on  the  surface.  If  a 
target  of  positive  elevation  is  present  at  range  R,  the  return  will  not  be  zero.  In  order  to  be 
able  to  examine  all  ranges  R  on  a  single  pulse,  one  can  restrict  the  occurrence  of  the  n:  to 
the  receiving  antenna  only  and,  making  use  of  phased-array  techniques,  perform  the  nul 
scanning  operation  in  the  receiver.  One  disadvantage  of  this  system,  as  compared  with  the 
narrow-beam  system,  is  that  whereas  the  narrow-beam  system  enables  one  to  simultaneously 
eliminate  the  clutter  and  measure  the  target's  altitude,  in  the  null-scanning  system  one  can  only 
determine  the  target's  altitude  when  it  is  at  a  range  where  there  is  no  clutter.  In  order  to  detect 
a  target  in  clutter  and  simultaneously  measure  its  altitude  with  a  null-scanning  system,  one  needs 
to  employ  two  nulls,  the  positions  of  which  can  be  adjusted  independently. 

A  system  which  makes  .no  demands  at  all  on  the  angular  >'esolution  of  the  vertical  antenna 
pattern,  but  demands  an  exceedingly  fine  range  resolution,  is  based  on  the  use  of  the  forward¬ 
scattering  properties  of  the  surface.  Assuming  that  conditions  are  such  that  the  direct  and  in¬ 
direct  rays  are  both  sufficiently  strong,  and  making  use  of  the  fact  that  the  altitude  H  is  posi¬ 
tive  if  and  only  if  the  pathlength  difference  A  is  positive,  one  can  filter  out  signals  of  zero 
altitude  by  filtering  out  signals  of  zero  pathlength  difference.  One  simple  way  of  achieving 
this  (but  by  no  means  the  ideal  way)  is  to  compare  the  range  autocorrelation  function  of  the  total 
received  signal  with  the  autocorrelation  function  of  the  transmitted  waveform.  Whereas  the 
autocorrelation  functicn  of  the  clutter  component  will  (a.ssuming  appropriate  statistics)  lend  to 
be  the  same  as  the  autocorrelation  function  of  the  transmitted  waveform,  the  autocorrelation 
function  of  the  target  component  will  have  additional  structure  due  to  the  various  paths.  For 
example,  if  the  transmitted  waveform  is  a  simple  pulsed  sinusoid  of  pulselength  T,  whereas 
the  autocorrelation  function  of  the  clutter  component  will  have  a  peak  at  r  ~  0  and  will  tend  to 
vanish  by  the  time  T  -  T,  the  autocorrelation  function  of  the  target  component  will  (assuming 
the  distortion  of  H  is  sniaJl)  have  peaks  at  t  =  0,  t  -  ±A/c,  and  t  =  ±2A/c.  Thus,  to  filter  out 
signals  of  zero  altitude,  one  need  .nerely  set  a  threshold  on  the  level  of  tne  autocorrelation 
function  in  a  r-region  away  from  t  =  0.  Since  the  peaks  at  t  =  ±A/c  and  t  ±2A/c  enable  one 
not  only  to  detect  the  existence  of  the  target,  but  also  to  measure  the  value  of  A,  this  system, 
like  the  narrow-hcamwidth  system,  allows  one  to  simultaneously  eliminate  the  clutter  and  meas¬ 
ure  the  target's  altitude. 


On  the  whole,  KTI  clulter-rejeciion  techniques  have  received  very  little  attention.  Aside 
from  the  gross  mapping-out  procedures,  which  have  been  incorporated  into  a  va.'iety  cl  sy.s- 
ic-ins,  most  of  the  work  has  been  theoretical.  Very  little  effort  has  been  made  to  test  these 
ideas  empirically  or  to  develop  the  appropriate  equipment. 

2.  Point  Target  Indicators  (PTI) 

In  the  ET!  techniques,  the  clutter  is  eliminated  by  employing  a  radar  with  sufficient  spatial 
resolution  to  resolve  the  target  from  the  earth's  surface,  and  then  by  "gating  out"  the  clutter. 
-Another  set  of  techniques  for  clutter  elimination,  based  on  the  differences  between  the  target 
and  the  surface  with  respect  to  their  size,  makes  use  of  spatial  resolution  to  'Thin  out"  the 
clutter.  If  the  target  is  a  point  target  and  tiie  surface  constitutes  a  more  or  less  continuously 
extended  source  of  edioes  which  cannot  be  resolved  from  the  target,  the  effects  of  the  clutter 
can  be  reduced  by  spreading  the  clutter  energy  out  over  an  increased  numlx.'i‘  of  resolution 
"boxes."  The  most  widely  used  technique  of  this  sort  consists  of  shortening  the  radar  pulse 
(either  directly  or  by  pulse  compression)  or  narrowing  the  beamwidth,  llu-reby  decreasing  tlic 
illuminated  area  (i.e.,  that  area  on  the  surface  which  contributes  to  the  received  signal  at  a 
given  instant  of  observation).  As  long  as  the  size  of  the  resolution  volumi-  remains  larger  than 
the  target  so  that  the  target  itself  does  not  appear  extended,  a  decre.ise  in  the  area  of  illumina- 
tiini  will,  on  the  average,  result  in  an  increase  in  the  larget-io-i  latter  ratio. 

Another  PTI  technique  makes  use  of  pulse-to-pulse  frequency  jumping  and  coherent  integra¬ 
tion.  Whereas  the  reflective  properties  of  a  small  target  vary  slow  ly  vvith  changes  in  fretjuency, 
the  reflective  properties  of  a  large  complex  target  like  the  earth's  surface  vary  rapidly  and 
randomly  w'ith  such  changes.  Thus,  if  the  frequency  jumps  are  chosen  properly,  whm  eas  the 
target  signal  will  tend  to  integrate  coherently,  the  clutter  signal  will,  like  noise,  ti-nd  to  integrate 
incoherently  Assuming  that  the  source  of  the  clutter  signal  <an  be  represented  us  a  eoUection 
of  statistically  independent  random  scatterers,  one  can  show  tliat  the  clutter  signal  will  be  sure 
to  be  decorreiated  provided  the  frequency  jump  is  greater  than  or  equal  to  the  bandwidth  of  the 
envelope  of  the  transmitted  signal. 

3.  Moving  Target  indicators  (MTI) 

In  addition  to  taking  advantage  of  the  differences  in  the  positions  of  the  target  and  the  earth's 
surface  as  is  dene  in  the  ETI  techniques,  one  can  attempt  to  eliminate  the  clutter  by  making  use 
of  the  differences  in  their  velocities.  Since,  in  general,  the  target  will  have  a  nonzero  velocity 
relative  to  the  car'h's  surface,  it  is  theoretically  possible  to  eliminate  the  effects  of  clutter  by 
designing  a  radar  which  is  responsive  to  velocity  vlifferences.  These  velocity-discrimination 
techniques  can  he  conveniently  divided  into  two  categories:.  Doppler  MTI  techniques  and  volume 
MTI  techniques.  In  the  Doppler  techniques,  the  velocity  discrimination  is  based  on  changes  in 
the  range  of  a  target  of  the  ord{?r  of  a  wavelength  of  the  t.-ansmitted  frequency.  In  the  volume 
techniques,  the  veloc-ity  discrimination  is  based  on  changes  in  the  position  of  the  target  of  the 
order  of  the  dimensions  of  the  illuminated  volume.  Inasmuch  as  the  velocities  are  usually  too 
small  to  cause  the  target  to  move  from  one  resolution  box  to  another  between  pulses,  volume 
MTI  systems  must  usually  be  based  on  the  changes  in  position  which  occur  between  scans  of  the 
antenna.  The  Doppler  systems,  whuh  hai’c  received  the  greater  attention  of  the  two,  can  be 
further  subdivided  actording  to  whmlier  they  are  coherent  or  incoherent  and  continuous  or  pulsed. 
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In  the  coherent  systems,  the  discrimination  is  based  on  information  concerning  the  Doppler- 
shifted  frequencies  themselves,  and  the  clutter  is  rejected  by  suppressing  those  frequencies 
which  cori  esponci  to  the  v<>loci!y  of  the  earth  relative  to  the  rauuj-.  In  the  incoherent  systems, 
the  discrimination  is  based  on  information  concerning  the  beat  frequencies  between  the  Doppler- 
siiincii  frequencies,  and  the  clutter  is  rejected  by  suppressing  those  signals  which  do  not  contain 
the  appropriate  beat  frequencies.  Although  the  incoherent  systems  are  simpler  in  that  they  do 
not  require  information  on  the  velocity  of  the  earth  relative  to  the  radar,  they  are  of  use  only 
to  the  extent  that  the  locations  of  the  clutter  source  are  predictable.  (Inasmuch  as  the  decision 
as  to  the  presence  of  a  target  is  based  on  the  observation  of  a  beat  frequency  between  the  target 
and  clutter,  if  the  system  is  used  in  regions  where  there  is  no  clutter,  it  will  suppress  the  tar¬ 
get  signal.)  The  pulsed  systems  differ  from  the  continuous  ones  in  that  the  pulsed  systems  only 
provide  samples  of  the  information  provided  by  the  continuous  ones  and  thus  are  afflicted  by 
ambiguities  (leading  to  the  so-called  "blind -speed"  and  "fold-over"  problems). 

Independent  of  the  particular  type  of  Doppler  system  employed,  the  extent  to  which  the  sys¬ 
tem  will  be  successful  will  depend  on  the  location  and  extent  of  the  clutter  spectrum.  In  a  con¬ 
tinuous  system,  the  ability  to  detect  targets  in  clutter  will  depend  on  the  relative  characteristics 
of  the  clutter  spectrum  and  ihe  target  spectrum.  In  a  pulsed  system,  the  performance  will  also 
depend  on  the  relation  of  the  clutter  spectrum  to  the  pulse-repetition  frequency.  In  general,  the 
fluctuations  entering  into  the  clutter  signal  will  be  of  two  types:  those  which  result  from  varia¬ 
tions  in  the  scattering  elements  on  the  surface  (the  so-called  "internal  clutter  process")  and 
those  which  result  from  variations  in  the  radar.  If  the  radar  is  airborne,  an  important  item 
in  the  latter  category  is  the  motion  of  t.he  antenna.  Assuming  that  the  antenna  is  bei.ng  both 
scanned  and  translated,  one  must  consider  not  only  ihe  spectral  broadening  caused  by  the  internal 
clutter  variations,  but  also  the  broadening  caused  by  variations  in  the  illumination  of  the  scatter¬ 
ing  elements  and  by  the  differential  Doppler  shifts  associated  with  scattering  elements  at  differ¬ 
ent  angles. 

4.  Resolving  the  Clutter 

In  all  the  techniques  previously  discussed,  the  earth's  surface  has  been  viewed  macroscopi- 
cally.  Another  set  of  techniques  which  can  be  used  to  reject  the  clutter  is  based  on  the  fact  that, 
in  many  cases,  the  source  of  the  clutter  signal  will  not  actually  be  continuous,  but  will  consist 
of  a  collection  of  closely  spaced  discrete  targets.  If  the  radar's  spatial  resolution  can  be  made 
sufficiently  fine  to  resolve  these  discrete  targets,  and  the  extent  of  the  aircraft  target  is  suf¬ 
ficiently  limited  to  permit  the  aircraft  to  fit  in  between  these  discrete  targets,  then  the  prob¬ 
lem  of  determining  whether  or  not  an  aircraft  is  present  is  reduced  to  a  pure  miiltiple-target- 
identification  problem.  Assuming  that  this  problem  cannot  be  solved  by  ETI  or  MTI  techniques, 
how  difficult  it  will  be  to  solve  will  depend  on  the  degree  to  which  the  detailed  scattering  prop¬ 
erties  of  the  aircraft  target  differ  from  those  of  the  discrete  earth  targets. 

5.  Theoretical  Work  on  Optimum  Anticlutter  Waveforms 
and  Optimum  Anticlutter  Filters 

In  recent  years,  increasing  thought  has  been  given  to  the  problem  of  determining  the  optimum 
waveform  and  the  optimum  filter  for  purposes  of  clutter  rejection.  Although  these  efforts  have 
not  yet  culminated  in  a  comprehensive  theory  of  anticlutter  design,  some  of  the  results  are  of 
considerable  interest  and  should  be  of  definite  use  to  the  radar  designer. 
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One  line  of  research,  exemplified  by  Applebaum  and  Howells,^  V/esterfield,  Prager  and 
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Stewart,  and  Fowle,  Kelly  and  Sheehan,^  is  based  on  the  use  of  the  time-freoyency  ambiguity 
function.  These  analyses  start  by  assuming  that  the  receiver  consists  of  a  bank  of  filters  matchec 
(in  frequency  and  tlrne)  to  an  aruiirary  waveform  u{t)  competing  against  white  Gaussian  noise. 
Letting  t  be  the  delay  variable  and  <p  the  Doppler  shift  variable,  one  measures  the  "goodness" 
of  u(t)  for  purposes  of  clutter  rejection  by  the  extent  to  which  the  volume  under  the  ambiguity 
function  'if{T,<p)  -  |  J  u(t)  u'*'(t  +  t)  exp{-2iriv>t)  dt|^  is  prevented  from  overlapping  the  volume 
under  the  probability  surface  describing  how  the  clutter  energy  Is  distributed  in  t  and  <fi.  Al¬ 
though  the  maximum  value  of  ♦  and  the  total  volume  under  ♦  are  independent  of  the  waveform 
design  and  depend  only  on  the  total  energy  in  u(t),  the  volume  in  the  central  spike  cf  ♦,  and  the 
manner  in  which  the  remaining  volume  is  distributed  in  the  (t,  (/>)  plane,  are  determined  by  the 
detailed  structure  of  u(t).  If  the  clutter  energy  is  distributed  uniformly  throughout  the  (t,  tp) 
plane,  then  all  waveforms  become  equivalent.  If  the  clutter  energy  is  concentrated  in  a  localized 
region  and  the  target  of  interest  lies  sufficiently  far  outside  this  region,  then  a  "good"  waveform 
is  one  in  which  all  the  volume  under  ♦  is  concentrated  in  the  central  spike.  If  the  clutter  energy 
is  concentrated  in  a  localized  region  and  the  target  lies  inside  this  region,  then  a  "gcod"  wave¬ 
form  is  one  in  which  ❖  has  a  narrow  central  spike  with  as  much  volume  as  possible  existing 
outside  this  region.  When  looked  at  in  this  way,  the  problem  of  anticlutter  waveform  design  be¬ 
comes  quite  similar  to  the  problem  of  anticlutter  antenna  design.  In  the  lost  case,  for  example, 
ihe  basic  idea  for  both  antennas  and  waveforms  is  to  "buy"  target -to-clutter  ratio  at  the  "cost" 
of  ambiguity.  In  general,  this  trade  can  be  effected  by  the  introduction  of  periodicity.  Making 
the  aperture  periodic,  i.e.,  chopping  up  the  aperture  and  separating  the  pieces  in  space,  enables 
one  to  decrease  the  width  of  the  main  lobe  but  results  in  strong  sidclobes.  Making  the  waveform 
periodic  enables  one  to  decrease  the  width  of  the  central  spike  in  ♦  but  results  in  "blind-speed" 
ambiguities  and  "second-time-around"  ambiguities.  For  a  detailed  understanding  of  this 
ambiguity-function  approach  and  for  specific  quantitative  results,  the  reader  is  referred  tu  the 
above-mentioned  references. 

A  second  line  of  theoretical  work  oriented  toward  the  development  of  a  clutter-rejecUon 
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theory  is  based  on  Dwork's  matched  filler  theorem  for  colored  Gaussian  noise  and  is  exem¬ 
plified  by  Urkowitz^^  and  Manasse  Letting  y(w)  be  the  Fourier  transform  of  the  received 
signal  and  N(w)  the  power  spectium  of  the  noise,  Dwork  has  shown  that  the  transfer  function  'i(w) 
of  the  linear  filter  which  maximizes  the  output  peak  S/N  ratio  is  given  by 

Hlw;  =  exp{-iwT)  (38) 

where  T  is  a  conveniently  chosen  time  delay  (assumed  hereafter  to  be  zero).  The  peak  S/N 
ratio  obtained  with  this  filler  is  given  by 


N  opt 


Assume  now  that  (a)  the  duration  of  the  signal  transmitted  fay  the  radar  is  sufficiently  short  with 
respect  to  the  movements  of  the  radar,  target,  and  clutter  that  everything  can  be  regarded  as 
fixed  during  the  duration  of  ihe  signal;  (b)  the  signal  received  from  the  target  is  merely  a  de¬ 
layed  and  aue.nuated  version  os  the  transmitted  signal;  (c)  the  statistics  of  the  clutter  signal  are 
Gaussian  and  the  power  spectrum  of  the  clutter  (cbtai.ned  from  the  fluctuations  of  the  clutter 
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signal  with  range)  is  proportional  to  the  energy  spectrum  of  the  transmitted  signal.  As  will  be 
seen  in  Sec.  Ill,  the  last  of  these  assumptions  (the  second  part  of  which  is  identical  to  the  assump¬ 
tion  mentioned  above  that  the  autocorrelation  function  of  the  clutter  is  proportional  to  the  auto¬ 
correlation  function  of  the  transmitted  pulse)  is  reasonably  well-satisfied  for  the  case  of  sea 
clutter.  With  these  assumptions,  the  above  theorem  implies  that  the  linear  filter  which  maxi¬ 
mizes  the  peak  signal-to-clutter-plus-noise  ratio  has  a  transfer  function  given  by 


H(w)  = 


bIxmPTn^ 


where  is  the  receiver  noise  power  per  cycle,  X(w)  is  the  Fourier  transform  of  the  trans¬ 
mitted  signal,  and  B  is  a  constant  of  proptortionality  describing  the  total  clutter  power.  The 
peak  signal-to-clutter-plus-noisc  power  ratio  obtained  with  this  filter  is  given  by 


(  S  ,  ^  r“  |X(m1|^ 

C  +  N  opt  2w  J_^  B|X((j)!"  +  N^/2 


db) 


(41) 


where  A  is  the  amplitude  of  the  returned  target  signal. 

If  the  clutter  power  is  zero.  H(w)  reduces  to  the  complex-conjugate  filter  H(u)a  X*((i?), 
and  (S/(C  +  N)l^pj  becomes  (S/(C  +  Nll^p^  =  <®/N^opt  ®  jCl  |X(w)|^  dw.  Aside  from 

A  and  N^,  the  quantity  (S/N)^p^  depends  only  on  the  energy  in  X(h>).  If  the  noise  power  is  zero 
(the  case  considered  by  Urkowitz),  H(a))  reduces  to  the  inverse  filter  H(w)  «  l/X(u),  and  (S/(C  •*- 
N)|opt  becomes  fS/(C  +  N)|^pj  *  (S/O^p^  «  (A^/B)  dw.  Thus  {S/C)^pj  is  completely  inde¬ 
pendent  of  X(w)  and  is  infinite.  This  last  result  can  be  interpreted  in  more  familiar  terms  by 
noting  that  the  use  of  the  Inverse  filter  Is  theoretically  equivalent  to  using  a  Dirac  delta  function 
as  the  transmitted  signal  and  decreasing  the  illuminated  area  to  zero.  This  equivalence  can  be 
demonstrated  formally  by  representing  the  antenna  and  surrounding  space  as  a  fixed  linear  filter 
(referred  to  as  the  "space  filter"  in  the  Introduction)  and  reversing  the  order  of  application  of 
this  filter  and  the  filter  H(w).  Assuming  that  the  construction  of  the  inverse  filter  is  limited  to 
a  finite  frequency  interval  Wj<  |w|  <  w^  and  that  the  energy  outside  this  interval  is  attenuated 
to  zero,  Urkowitz  has  shown  that  (S/Cl^p^  «  (w^  -  « j).  In  other  words,  the  ability  to  eliminate 
the  clutter  depends  linearly  on  the  bandwidth  over  which  the  inverse  filter  can  be  constructed. 

Returning  to  the  general  case  (considered  by  Mnnasse*^)  and  assuming  that  both  clutter  and 
noise  are  present,  one  sees  that  H(w)  is  a  compromise  between  the  complex-conjugate  filter 
X'f'lw)  and  the  inverse  filter  l/X(w),  and  that  (S/(C  +  N)1  ,  depends  on  both  the  total  energy  in 
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X(w)  and  on  the  shape  of  |X(w)|  .  In  regions  of  w  where  |X(w)r  «  the  contribution  to 
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(S/(C  +  N))jjpj  is  given  by  (AVeN^)  f  |X(w)|  ®nd  thus  depends  only  on  the  total  energy  of 
X(w)  in  those  regions.  In  regions  of  w  where  |X(w)|^  »  N^ZB,  the  contribution  to  (S/(C  + 
N)|opt  is  given  by  (A^/2eB)  f  dw  and  thus  is  entirely  independent  of  X(w)  and  depends  only  on 
the  extent  of  those  regions.  It  is  clear,  therefore,  that  If  one  wants  to  increase  the  value  of 
[S/(C  N)|^p^  by  changing  the  waveform  (subject  to  a  fixed  total  energy),  one  should  take  the 

superfluous  energy  in  the  |X(w)|^  »  N^/ZB  regions  and  put  it  in  the  |X(w)|^  «  N^/2B  regions. 
More  precisely,  Manasse  has  shown  (through  variational  methods)  that  in  order  to  maximize 
(S/(C  4  N)I^p^,  the  transform  X(w)  should  be  chosen  such  that  |X(w)|^  is  flat.  If  X(w)  is  con¬ 
strained  to  be  zero  outside  the  interval  W|  <  |w|  <  W2  and  to  have  a  total  energy  (1. .  it  will 
maximize  (S/(C  +  N)|  ^p^  when 
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X(ul) 

|2  ..  irfl 

for  w  j  ^  1  w  j  a>2 

X(w) 

otherwise. 

(42) 

The  choice  of  a  flat  jX(u))|^  makes  the  clutter  powet  spectrum  have  the  same  shape  as  the  noise- 
power  spectrum  and  ehminaies  the  need  for  compromise  m  the  optimum  filter  Hfu,)  [H(u  l  now 
being  give-n  by  “  X'M.w)  “  l/Xfw)).  The  value  of  |S/<C  +  ob’-'^med  with  sucli  a  wave  ¬ 

form  is  given  by 

( — ^ — ) 

C  +  N  opt,  opt 

Large  values  of  fl  will  have  the  effect  of  reducing  the  noise  interferen  .-e,  and  large  values  of 
have  the  effect  of  reducing  the  clutter  interference.  For  a  pulse  ratiar  whuh  is 
peak  power  limited,  this  implies  that  the  transmitted  pulse  should  h.ive  a  large  time-bandwidth 
product  {i.e.,  the  pulse  should  be  "coded").  'I'wo  techniqin*s  for  ac-hieving  such  a  signal  (and 
also  satisfying  the  requirement  for  a  flat  energy  spectium)  are  to  use  long-duration  pulses  of 
broad-band  white  noise  or  long-f-uration  pulses  with  a  linearly  swept  carrier  frequency  of  large 
deviation.  •• 

Unfortunately,  although  both  of  these  approaches  to  the  clutter-rejectior.  problem  are  capable 
of  making  substantial  contributions  to  a  general  anticlutter  theory,  at  present,  they  both  suffer 
from  some  serious  limitations  and.  in  no  sense,  can  be  regarded  as  constituting  complete  theories 
in  themselves.  Aside  from  the  fact  that  both  approaches  disregard  the  specific  reflective  charac¬ 
teristics  of  the  clutter,  in  neither  approach  is  any  consideration  given  to  the  problem  of  choosing 
an  optimum  antenna  patti-rn.  if  the  clutter  is  three  dimensional,  ns  in  the  case  of  rainstorms, 
the  latter  omission  may  not  be  too  important  in  the  sense  that  the  total  opiimization  process  may¬ 
be  factorable  and  one  may  be  able  to  optimize  over  the  antenna  pattern  in  a  separate  operation. 
However,  in  the  case  of  earth  surface  clutter,  where  range  resolution  is  a  function  both  of  the 
transmitted  waveform  and  the  vertical  antenna  pattern,  this  omission  may  be  quite  impoitant. 

In  addition  to  these  two  limitations,  which  apply  equally  well  to  the  ambiguity-function  ap¬ 
proach  and  the  colored-matched-filter  approach,  the  ambiguity-function  approach  suffers  from 
its  a  priori  assumption  about  the  nature  of  the  receiver.  The  optimum  waveforms  which  can  be 
determined  through  this  approach  are  only  optimum  -with  respect  to  the  given  receiver.  No 
optimization  is  performed  over  the  receiver-waveform  pair.  In  general,  one  may  regard  the 
target-in-clutter  problem  either  as  a  target -in-noise  problem  or  as  a  multiple-target  problem. 
Insofar  as  the  former  philosophy  is  adequate,  the  noise  against  which  the  target  signal  is  com¬ 
peting  is  the  receiver  noise  plus  the  clutter  noise  and.  in  view  of  the  Dwork  theorem,  the  matched- 
filter  bank  assumed  in  the  ambiguity-function  argument  is  obviously  not  optimum.  Insofar  as  the 
latter  philosophy  is  adequate  and  one's  objective  is  to  detect  each  target  and  then  reject  the  clut¬ 
ter  targets  on  the  grounds  of  certain  parameter  values  <e.g.,  range  and  velocity),  the  matched- 
filter  bank  is  still  inadequate  bi  l  ause  it  is  not  properly  matt  hed  to  the  multiple-target  situation. 
Thus,  in  either  case,  it  is  clear  that  th«*  assumed  receiver  in  the  ambiguity-function  approach 
is  not  optimum,  and  therefore,  that  the  optimum  waveforms  which  result  from  this  approach  will 
not  really  be  optimum,  .\noiher  limitation  of  this  approach  is  that  it  incorporates  the  effects  of 
velocity  in  such  a  way  as  to  be  valid  only  for  narrow-band  signals  For  signals  of  appreciable 
bandwidth,  one  must  take  account  of  the  dependence  of  the  Doppler  shift  on  frequency  over  the 
bandwidth  of  the  signal. 
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The  colored-matched>filter  approach,  on  the  other  hand,  although  it  optimizes  over  both 
the  receiver  and  the  waveform,  suffers  from  the  fact  that  it  totally  ignores  the  velocity  informa¬ 
tion.  Moreover,  the  results  of  this  approach  u'ill  only  be  useful  to  the  extern  thai  the  clutter 
can  be  regarded  as  continuous  and  the  cross  sections  of  the  individual  scatterers  that  constitute 
the  source  of  the  clutter  signal  can  be  regarded  as  infinitesimal  in  comparison  to  the  cross  sec¬ 
tion  of  the  target.  Once  the  clutter  is  resolved,  the  statistics  of  the  clutter  signal  will  no  longer 
be  Gaussian  and  the  matched-filter  theorem  will  no  longer  be  applicable.  Also,  to  the  extent 
that  the  irregularities  on  the  reflecting  surface  have  an  ordered  structure  and  are  not  random, 
the  power  spectrum  of  the  clutter  will  no  longer  be  merely  a  replica  of  the  energy  spectrum  of 
the  transmitted  sigi'al. 

In  the  event  that  the  forward-scattering  conditions  ensure  a  strong  reflected  ray,  both  o'. 
the  above  approaches  also  need  to  be  modified  to  include  the  earth-modification  factor  E  de¬ 
scribed  in  Sec.  I.  In  the  ambiguity-function  approach,  the  presence  of  the  reflected  ray  implies 
that  the  receiver  must  consist  of  a  multidimensional  filter  bank  to  account  for  K,  A,  and  A'  as 
well  as  for  t  and  <p,  and  the  ambiguity  function  ♦  will  be  a  function  of  K,  A,  and  A'  ns  well 
as  of  T  and  <p.  In  the  colored-matched-filter  approach  (still  ignoring  velocity),  the  existence 
of  a  reflected  ray  implies  that  the  ideal  receiver  consists  of  a  filter  bank  of  the  form 
E*(w,K.  A)  exp(-iu)T)/{N^/2  +  B|X(aj)|^).  where  there  is  a  separate  matched  filter  for  each 
value  of  A  and  K. 


III.  SCATTERING  HIOM  THE  OCEAN 

In  Sec.  II,  the  problem  of  surface  roughness  was  considered  without  reference  to  the  specific 
reflective  characteristics  of  the  surface.  In  the  present  section,  attention  will  be  devoted  to  the 
reflective  properties  of  the  ocean.  This  surface  has  been  chosen  for  detailed  consideration  be¬ 
cause  (a)  it  is  a  surface  of  great  practical  importance;  (b)  despite  its  complexity,  it  is  simpler 
than  most  other  portions  of  the  earth's  surface;  (c)  it  has  received  considerable  attention,  both 
experimentally  and  theoretically.  In  accordance  with  the  fact  that  the  clutter  phenomenon  is 
usually  of  greater  concern  to  the  radar  designer  than  the  corruption  of  the  forward-scattered 
signal,  the  main  concern  will  be  with  the  clutter  phenomenon.  The  discussion  will  be  oriented 
toward  the  problem  of  predicting  the  sea-clutter  signal  for  a  narrow-band,  pulsed,  airborne, 
early-warning  radar  of  broad  vertical  bcamwidth  and  wavelength  \  in  the  region  1  to  100  cm. 
Inasmuch  as  the  early-warning  function  tends  to  focus  interest  on  the  longer  ranges,  special 
attention  will  be  given  to  grazing  angles  o  in  the  region  0'  to  20".  A  photograph  illustrating  the 
appearance  of  sea  clutter  on  the  PPI  scope  of  an  airborne  early-warning  radar  is  shown  in 
Fig.  15. 

In  CT-di-r  to  understand  the  sea  clutter  phenomenon  and  to  be  able  to  predict  the  properties 
of  the  clutier  signal  for  a  variety  of  radar  configurations,  it  is  necessary  not  only  to  make  quan¬ 
titative  measurements  of  the  received  signal,  but  also  to  make  quantitative  measurements  of 
the  ocean  surface  responsible  for  this  signal.  If  the  electromagnetic  properties  of  the  sea  are 
known  and  .«uch  factors  as  the  spray  and  foam  can  be  ignored  (so  that  the  sea  surface  actually 
constitutes  a  continuous  interface  between  the  air  and  the  water),  the  scattering  properties  of 
the  sea  surface  will  be  determined  by  specifying  its  shape  as  a  function  of  time,  l.etting 
z(x.  y.  t)  denote  the  height  of  the  surface  at  the  point  (x.y)  at  the  time  t,  one  can  give  a  com¬ 
plete  statistical  description  of  this  shape  by  specifying  the  set  of  probability  density  functions 
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F}g.  15.  AppearatK«  of  sea  clutter  of»  PPI  scope  of  airborne  early-worning 
rodcr.  Strobe  is  set  at  135  tm,  30  nm  short  of  horizon. 

on  the  set  of  variables  {zlx,  y,  t))  .  Thus,  in  the  ideal  clutter  experiment,  along  with  the 

data  on  the  received  signal,  data  would  be  presented  on  these  probability  density  functions  de¬ 
scribing  the  ocean  surface.  (If  the  spray  and  foam  cannot  be  ignored,  then  of  course,  even  these 
data  would  be  inadequate  and  one  would  need  to  supplement  it  by  a  probabilistic  description  of 
various  spray  and  foam  parameters.)  Unfortunately,  in  most  of  the  clutter  experiments  reported 
in  the  literature,  the  amount  of  data  on  the  sea  surface  is  extremely  limited.  With  few  excep¬ 
tions,  the  most  information  that  is  available  is  an  estimate  of  wind  velocity,  wave  heigh;  and 
direction,  and  a  qualitative  description  of  the  ocean's  appearanc  e.  In  some  cases,  no  data  at 
all  are  available.  Aside  from  the  fact  that  the  estimates  of  these  parameters  are  often  made 
very  crudel_ ,  and  the  parameters  being  estimated  are  often  defined  very  crudely,  it  is  obvious 
that  such  a  technique  is  basically  inadequate  for  describing  as  complex  a  surface  as  th  it  of  the 
ocean.  Although  the  hydrodynamic  constraints  on  the  shape  of  the  water  surface  tend  to  reduce 
the  randomness  of  the  suf-face  and  thus  make  some  of  the  information  in  the  above-described 
probability  density  functions  redundant,  these  c'onstraints  are  not  sufficiently  strong  to  enable 
one  to  replace  these  functions  by  a  few  simple  constants  such  as  wave  height  and  direction. 
Similarly,  although  the  momentary  local  wind  velocity  will  undoubtedly  be  correlated  with  cer¬ 
tain  properties  of  t.he  surface  (e.g..  the  ripple  structure),  it  is  an  insufficient  statistic  for  de¬ 
termining  the  surface  as  a  whole.  In  general,  the  condition  of  the  sea  surface  in  a  given  area 
and  at  a  given  time  will  be  a  function  of  the  history  of  the  wind  velocity  over  a  wide  region  of 
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the  ocean  surface  (as  well  as  the  currents,  tides,  ant!  the  proximity  and  tcpogr.iphy  of  the  ocean 

floor).  For  example,  it  is  well  known  that  the  occurrence  of  "swells"  in  a  given  deep  water  urea 

is  not  due  to  the  momentary  local  wind  velocity,  bat  rather,  to  wind  conditions  which  existed  in 

regions  far  away  and  at  an  earlier  time.  These  waves  have  outrun  the  wind  which  created  them 

and  are  in  the  process  of  decay.  As  far  as  the  waves  which  are  generated  locally  are  concerned, 

their  characteristics  are  determined  not  only  by  the  momentary  local  wind  velocity,  but  also  by 

the  length  of  time  during  which  the  wind  has  been  blowing  and  the  length  of  ocean  over  which  it 

has  been  blowing.  In  describing  data  on  the  clutter  signal,  the  symbol  "f"  will  be  used  to  denote 

the  "condition  of  the  ocean  surface."  In  the  ideal  experiment,  t  would  consist  of  the  set  of  ail 

probability  density  functions  previously  mentioned,  in  practice,  however,  it  must  be  interpreted 

as  the  wave  height  or  wind  velocity,  or  even  more  vaguely,  as  the  "degree  of  roughness."  For 

some  detailed  results  on  the  actual  characteristics  of  the  ocean  surface,  the  reader  is  referred 

13 

to  works  of  hydrodynamicists  and  oceanographers. 

The  discussion  in  this  section  will  be  divided  into  three  parts.  In  Sec.  III-A,  an  attempt  will 
be  made  to  provide  a  brief  summary  of  the  major  experimemal  results  on  the  clutter  signal. 

When  considering  any  of  the  statements  in  this  summary,  the  leader  should  be  aware  of  the  fol~ 
lowing  facts.  F'irst,  the  data  on  sea  clutter  are  inconsistent  and  incomplete  and,  in  the  usual 
sense  of  the  word,  there  is  no  such  thing  as  a  "typical"  resuit.  Although  a  certain  class  of  re- 
sult.s  may  occur  more  frequently  than  any  other  single  cla.ss,  there  is  no  class  which  occurs 
more  frequently  than  the  total  of  all  other  classes  combined,  unless  one  defines  this  class  so 
generally  as  to  be  almost  meaningless.  Second,  in  view  of  these  inadequacies,  in  order  to  pre¬ 
sent  a  general  summary  of  the  data,  one  must  either  ia)  violate  the  empirical  facts  !?nd  ignore 
large  portions  of  the  available  data;  (b)  state  the  summary  in  such  vague  terms  as  tc  make  the 
summary  almost  meaningless;  or  (c)  perform  a  massive  statistical  analysis  of  thousands  of 
curves,  weighting  each  curve  according  to  the  results  of  a  comprehensive  analy.«is  of  the  de¬ 
tailed  conditions  under  which  the  curve  was  obtained.  Understandably,  no  one  to  date  has  per¬ 
formed  the  analysis  required  by  (c).  In  this  report,  the  writer  has  done  his  best  to  compromise 
between  evils  (a)  and  (b).  The  summary  has  been  obtained  by  looking  at  a  large  number  of  ex¬ 
perimental  studies  and  represents  a  syntnesis  of  the  writer's  impressions  of  the  mein  results 
Specific  references  will  frequently  be  omitted,  but  many  of  the  remarks  will  be  illustrated  by 
concrete  experimental  data.  In  Sec.  III-B,  some  theoretical  concepts  will  be  introduced  for  use 
in  interpreting  the  experimental  data.  Since  there  is  no  unified  theoretical  model  available  which 
is  consistent  both  with  a  realistic  description  of  the  ocean  surface  and  with  the  clutter  data,  the 
best  the  radar  designer  ^an  do  in  trying  to  predict  the  ctutter  signal  for  an  l.ypothesized  radar 
is  to  combine  the  results  obtained  from  a  study  of  the  empirical  findings  with  results  computed 
from  various  fragmentary,  phenomenologic  il  models.  The  concepts  ehosen  for  distussion  have 
been  selected  for  their  popularity  with  researchers  who  have  actually  "dirtied  their  hands"  in 
the  sea-clutter  problem  and  have  been  forced  tc  make  concrete  recommendations  on  raour-tk sign 
problems.  In  .Sec.  III-C,  a  brief  summary  wii!  be  gi.'en  of  soi  of  the  results  on  forward  scat¬ 
tering  from  the  ocean. 

Those  who  are  familiar  w.ih  the  subject  of  sea  clutter  will  note  that,  with  few  exceptions, 
the  results  presented  in  Secs  !lI-.\  and  B  represent  only  a  very  modest  advance  over  those  pre¬ 
sented  by  Goldstein  in  the  years  194*5  tc  195!.*'*”^*’  This  us  not  meant  to  imply  that  no  important 
work  has  been  done  since  that  time  t)n  ihe  con'rary.  th*.  picture  of  s«-a  clutter  has  been  extenu“d 
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artd  refined  in  a  variety  cf  ways.  For  example,  (a)  data  are  now  available  at  much  lower  fr-t  • 
quenciea,  (b)  more  work  has  been  done  with  coherent  spectra;  (c)  the  dependence  on  polariza¬ 
tion  has  been  explored  more  fully;  td)  m<»re  careful  data  analysis  techniques  have  been  worked 
out;  (e)  specific  theor  etical  mv.dels  have  been  studied  in  more  detail.  Nevertheless,  if  one  con¬ 
fines  one'e  attentron  to  the  airborne  early-warning  problem,  tne  above  statement  still  holds: 
a  significant  portion  or  the  knowledge  which  is  now  useful  to  the  radar  designer  can  be  found  in 
the  writing's  of  (Joldsieiu. 

Finally,  the  reade*^  should  be  aware  of  the  fact  that  the  references  cited  in  connection  with 
Ihis  discussion  arx?.-  in  no  sense,  intended  to  constitute  a  bibliography.  A  complete  htbliograpny 
on  this  subject  would  be  longer  than  the  discussion  itself. 

A.  Erqpei'lmentjtl  Detst  on  Seo  Clutter 

1.  Average  Power 

For  a  point  target  in  ft'ee  space,  the  dependence  of  the  received  power  on  the  characteris¬ 
tics  of  the  tari;ot  is  described  through  use  of  the  -arget  crosa  sectioA  6  “  jsl^.  Aside  from  the 
characteristics  of  the  target  proper,  this  tVirameler  depends  only  on  the  frequency,  the  polariza¬ 
tion,  and  the  angle  of  orientation;  it  is  indepenavnt  of  the  transmitted  power,  the  transmitted 
pulselen.gth,  the  antenna  function,  and  ;;;e  range.  This  independence  has  the  important  advan¬ 
tage  cf  allowing  one  to  extrapolate  results  obtained  cu  the  target  characteristics  with  cne  radar 
configuration  tc  those  of  another,  prc'dded  only  that  the  frequency,  polarization,  and  angle  are 
known,  one  wants  to  achieve  a  similar  inftependence  for  e.  parameter  desci'lbin;»  the  back- 
scatterutg  properties  of  tae  ccean  surface,  the  notion  of  cross  section  is  inadequate.  Since  the 
ocean  is  an  extended  target,  <r(oc<;aai  will  depend  on  how  much  cf  the  surface  is  illuminated 
(i.e.,  how  much  of  the  surface  contributes  to  the  returned  signal  at  a  piven  instant  of  observa¬ 
tion)  and,  therefore,  will  be  e  function  of  the  antenna  pattern,  th®  transmitted  waveform,  and 
the  range.  In  order  to  eliminate  this  dependence,  it  is  nece«cary  to  determine  how  (T(occ'in) 
varies  witn  the  illuminated  area  and,  making  use  of  this  k.iowledgs,  to  define  a  new  parameter 
describing  the  properties  of  the  ocean  surface  proper. 

l.et  A  denote  the  area  of  illumination  and  a  the  pulse -to-pulse  average  of  (/xocean).  It  has 
been  shown  experimentally  that,  over  a  wide  variety  of  conditions,  a  varies  lineat-ly  with  A: 

V  =  a^A  (44) 

(These  data  have  been  obtained  by  (a)  varying  the  transmiited  pulse^Argth,  (fa>  varying  the  antenna 
pattern,  and  (c)  varying  the  range  or  height  while  ktsping  the  grazing  angle  constant.}  Thus,  a 
natural  definition  for  the  sought-after  parameter  is  =  v/A,  the  cresj  section  per  unit  areu. 
This  parameter  plays  the  same  role  for  the  ocean  ss  the  ordinary  cross  seciion  a  plays  for 
the  point  target;  aside  from  the  properties  n?  the  target  itjiclf,  it  depends  only  on  the  frequency, 
the  polarization,  and  the  angle  of  orientation.  Letting  9  denote  the  horizontal  faeamwidth, 
the  vertical  beamwidth,  T  the  transmitted  pulselength,  a;id  a  the  graxing  angle,  ent  can  ap¬ 
proximate  t.be  area  of  illumination  A  (sec  Fig  I6{a-b)j;  by 

c  r 

A  =  R@  seen  when  lano  < 

7 

A  =  R  ©♦  CSC  a  when  '.ana  >  .  /45) 

V  i  /  C* 
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A  brief  description  wtii  now  oe  given  of  the  d'>pendence  of  the  cross  section  per  unit  urea  on 

the  grazing  angle  a,  the  transmitted  wavelength  X,  the  polarization  p,  the  "condition  of  the  sea 

surface"  s.  and  the  azimuth  with  respect  to  the  structure  on  the  sea  surface. 

For  all  but  tne  calmest  of  seas,  the  form  of  the  function  7  idb)  vs  logo  most  often  etuoun- 

o  “ 

tered  is  (hat  yhown  in  Fig.  17  (or^{db)  =  10  ‘Og^Qe^j.  Starting  at  very  small  o,  the  parameter 

increases  fairly  rapidly  with  n  until  a  '.>quals  the  "critical  angle"  Oj.  It  tnen  inei eases  much 

more  slowly  {if  at  all)  until  the  second  "critical  angle"  O2  is  reached.  Beyond  o^.  it  increases 

very  rapidly  until  just  before  a  =  r/Z,  whereupon  it  tends  to  flatten  out  again  The  total  range 

of  slopes  for  logo  vs  logo  {i.e..  exponents  for  or^  vs  a)  is  about  0  to  20.  In  the  region  a  <  a  .. 

the  slope  is  usually  between  2  and  5,  and  in  the  region  <  o  <  o^.  the  slope  is  usually  between 

0  and  2.  The  total  measurable  range  in  from  very  small  a  to  a  rr/2  may  cove:-  as  much  as 

80  db  with  being  as  great  as  10  to  20 db  at  a  =  ir/2.  The  actual  range,  of  course,  is  infinite 

in  decibels  since  a  must  equal  zero  at  the  horizon  (The  lower  measurable  limit  on  o  is  de- 
o  o 

termined  by  the  am.ount  of  receiver  noise.) 

Letting  (7^^  be  the  ai  ?rage  level  of  ir^  between  a .  and  and  letting  ir^ivv)  and  rr^dih)  be  the 

values  cf  a  for  vertical  and  horizontal  polarizations,  one  can  make  the  following  genei-al  state- 

^  -  1  •  i 

ment  about  the  variables  a  a,  ,  a  ,  and  (T  (vv)/<r  (hh):  thev  all  tend  to  be  monotomcallv  de- 

1  4  00  o  o  ‘ 

creasing  with  (a)  decreasing  X  and  (b)  increasing  t.  The  variable  cr^lir '2)  ali  >  appears  to  be 
monotonically  decreasing  with  increasing  t.  Presumably,  the  same  result  would  hold  for  de¬ 
creasing  X,  hut  very  few  data  are  available  on  o^{ie/Z)  for  low  frequencu  s.  The  value  of  ir<vv) 
ff(hh)  is,  of  course,  always  equal  to  unify  at  a  ■  ?r/2  since  vv  and  hh  are  identical  a’  »  2.  for 
all  wavelengths  in  the  region  1  cm  <  X  <  100  cm  and  all  but  the  smallest  ?,  one  finds  that  «  ,  >  20°. 


Typical  values  of  and  0^  for  X  -  10  cm  and  a  moderate  to  rough  sea  are  -  ISdh  and 


Rest  riding  atlemion  ‘o  cr  c^.  one  finds  that  “  X"”.  0  ^  n  4  In  terms  of  o  j(X)  and 

Cf  (X).  the  dependence  of  a  on  X  can  usuallv  be  described  bv  the  relations  a  a  x"^.  0  <  n<  2. 

and  «  x'*,  0  <  n  1.  'I'hus  the  smaller  n  in  <t^  ^  X*’’  tend  to  occur  when  a  is  sufficiently 

large  to  make  o  >  o  j{X).  and  the  larger  n  tend  to  <'ccur  when  this  condition  is  violated.  IT  one 

holds  X  fixed  and  examines  the  variation  of  a  and  a  .  with  T.  one  finds  that  a  ini-reases 

00  1  00 

rapidly  with  t  for  relatively  small  f  but  then  tends  to  "saturate,"  the  saturation  level  occurrint 
at  smaller  ?  for  smaller  X.  If  or  identifies  ?  with  the  wind  speed  W,  one  finds  a  i  VV'*. 

0<  3.  If  one  identifies  ?  wit  .  the  wave  height  H.  o  ^  appears  t(-  vary  with  H  as  it  ^  »  fi"*'. 

0<  n<  1.  Thus,  for  a  fixed  a,  the  smaller  n  in  a  x”'’  tend  to  occui  when  t  is  large,  and 
the  larger  n  when  f  is  small. 

For  large  X  and  small  f,  cr^,j(vv)/iT^{hh)  is  much  greater  than  unity,  sometimes  reaching  a 
value  of  more  than  30  db.  As  X  is  decreased  or  C  isincuased,  (/^(vvl/a^^dih)  decrease-s;  how¬ 
ever,  the  rate  of  decrease  diminishes  as  (r^^lvvl/er^lhh)  decreases.  .At  small  X  and  large  r. 

(7^{ vvl/oTothh)  approaches  unity  and  sometimes  even  become.s  a  few  lieiibels  less  than  unity  In 
general,  a  (hh)  is  found  to  be  moi-e  sensitive  than  a  (vv)  (o  variations  in  r .  .Also,  there  is  t  vi- 

dence  that  7  (vv)/7  (hh)  decreases  as  a  becomes  extremelv  small.  The  i  ross-polat  i/.at  ion 
0  0 

cross  sections  7  (vh)  and  7  (hv)  (which  theoreticaUv  should  be  the  same  heiau.se  of  rtciprm  it\) 
00 

are  usually  smaller  than  either  7^^(hh)  or  o^(vv). 

Finally,  assuming  that  the  direclion  of  the  wind  and  tlie  wave  strueture  on  the  surfate  ar<- 


approximately  the  same,  7^  (upwind)  .s  almo.st  always  a  few  decibels  higher  than  eolu-r  7^ 

(downwind)  or  7^  (crosswind),  and  (i-rosswind)  is  usually,  but  not  always  a  little  greater  than 

7  (downwind), 
o 
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In  tho  above  remarks,  the  dependence  of  the  average  backscattered  power  on  the  char¬ 
acteristics  of  the  surface  has  been  described  in  terms  of  the  average  cross  section  per  unit 
atc-u  cr^y  Another  parametei'  which  has  occasionally  been  used  and  which  measures  ihe  surface's 
(levialicn  liom  perfect  roughness  is  the  parameter  f^^  which  is  defined  as  the  ratio  of  tin-  powei- 
actually  observed  to  the  powe  that  would  be  observed  if  the  surface  scattered  isotropically.  It 
can  be  shown  that  f^^  and  are  rehited  liy  the  equation  2  sin  (n  )  r,(o). 

Somi  specific  results  illustrating  some  of  th»'  above  remarks  are  shown  in  I'igs,  18  through 
24.  I'igure  18  illustrates  the-  dependence  of  on  X  and  a  for  small  values  of  o.  I'igure  19 
illustrates  the  dependenci'  of  cx^  on  X.  u,  and  f  for  small  values  of  a.  h'igure  20  illustrates 
the  dependence  of  rr^  on  a  and  t.  for  large  values  of  ir.  h'igure  21  illustrates  thi*  flependence 
of  (T^  on  p,  and  o  ovei-  a  wide  range  of  a.  ^■igure  22  is  a  coru-elogram  for  the  pola  ri/.ai  ion 
ratio  (x^(vv)/(T^{hh)  taken  at  two  different  X  with  variable  f.  I'igure  23  shows  how  ix^  varies 
with  p  and  t  for  small  a.  l''igur-e  24(a-o)  shows  how  fy(o)  ix^j((r)/2  sino,  whii'h  was  observed 
to  be  cemstant  over  the  region  1‘  <  o  <  I',  varies  with  X  and  T.  The  smooth  curves  in  l''igs.  18 
and  22  are  theorrdical  and  are  discussed  in  .Sec.III-B-3  In  addition  to  the  soui’ces  cited  in 

connection  with  these  figures,  further  data  on  the  average  backscattered  power  can  bo  found 

14  23  24  2‘3 

in  Kerr.  Kai/in,  firtinl  and  X'aitlee,  and  VVillse,  Schlesingirr  and  .lohnsoti  (iis  well  as 

many  other  s). 

2.  Kluctuat  ions  — 

In  the  previous  sr.'clion,  the  clutter  signal  was  de!i<‘rib<-{l  in  ti-rms  of  ihe  avertige  powr-r 
pariiiTurier-  rx^.  In  this  serf  ion,  attention  will  be  given  to  how  ihe  eluiier  signtil  varies  (ti)  as  a 
function  of  range  within  ii  given  sweep  and  (b)  as  a  funr  tion  of  lime  at  ii  given  timge 

In  genet'iil,  ther'e  appear'  to  lx*  very  few  datii  on  the  fluctuiitions  in  I'ange.  However,  what- 
evi-r  data  therrr  are  indicate  thiit  if  one  <  onfin«'S  one's  samples  to  the  small  angle  region  tano  < 
<I>K/(cT/2)  where  the  area  of  illumination  is  defined  by  ’ho  transmitted  waveform,  the  auto- 
correliition  function  of  the  r  eturned  range  sweep  is  approximately  the  same  as  the  rrutocor r  eia - 
tiori  firnctiorr  of  lh<'  ir'irnsmitted  waveforrri.  In  addition,  it  has  bi-en  observed  thirt  for  small  «, 
the  signal  tenils  to  resolvr-  its<'lf  into  distinct  point -tar-get -like  echoes  with  substantial  regions 
of  /er-o  sigtiitl  occur'fing  between  these  echoes.  This  "spikyness"  lends  to  incr»-ase  with  (a) 
cU'creasing  grazing  angle  and  (b)  rict'reasing  area  of  illumination.  Also,  (c)  this  spikyness  is 
r-nitch  mrrre  pr'onounced  with  hor-i'/.onial  polarization  th;m  with  vertical  polarization.  There  is 
good  evidence  that  the  occur-r'cnce  of  these  spikes  is  related  to  the  occurrence  of  well-defined, 
sleep-crested  waves  on  Ihe  ocean  surfac*'. 

.Although  there  ar-e  considerably  more  data  on  flui-tuations  in  time  thiin  on  fluctuiitions  in 
range,  there  are  still  fewer  d:;l!i  lh;m  on  the  more  eii.sily  measured  Iii  general,  it  has  been 

found  that  these  fluctuations  lend  to  resolve  themselves  into  three  main  regions:  (a)  Ihe  "fast" 
flttctuaiions  (presumably  due  to  the  beating  of  signals  of  different  frequencies  returned  from 
different  elements  on  the  sea  surface  moving  with  different  velocities);  (b)  the  "slow"  fluctua¬ 
tions  (piesumably  due  to  the  changes  in  amplitude  caused  by  the  changes  in  shape  and  orientation 
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Fig.  21 .  a  (a, p,  t)  dofa.  X  =  24  cm.  (Affer  Macdonald.  ) 


Fig. 22.  ao(X,P,  t)  dofo.  0.7®.^  a  1 .3*. 
^fter  Goldstein.-^)  Nomben  on  curve  give 
H  in  feet  («»«  Sec.  lil-B-S). 


Fig.  23.  Oq(p,  £}  doia.  A  =  4.8cm,  a=4* 
(After  Sowing,  et  al.^l) 
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Fig.  24.  fo(X,f)doto.  p  =  hh,  a  inragion  1 .0^  < 
a  5.0®.  (After  Davies  ond  Macforicna .^2) 


of  thi'  individual  seatluring  elements,  the  passage  of  these  elements  through  the  area  of  ilhimtn- 
atton.  and  the  growth,  decay,  and  shadowing  of  these  elements):  (c)  the  '  very  slow"  fluctuations 
(presumably  due  to  changes  in  meteorulogic  and  oceanographic  conditions).  In  making  measure¬ 
ments  of  cTq,  the  general  practice  is  to  present  results  which  are  averages  over  the  fast  fluctua¬ 
tions,  but  not  over  tlie  very  slow  fluctuations.  The  extent  to  which  the  slow  fluctuations  ar(> 
averaged  out  will  depend  on  the  detailed  ciiaracteristics  of  the  device  used  to  measure  and  on 
the  amount  of  averaging  done  on  these  measurements  before  presentation  of  the  final  results. 

If  the  transmitted  wavelength  is  sufficiently  short  (e.g.,  tOcm).  the  fast  and  slow  flui  tuations 
are  easily  d  stmguishable;  however,  if  the  transmitted  wavelength  is  long  (e.g..  100cm),  thi-se 
two  regions  will  tend  to  merge 

With  regard  to  the  fast  fluctuations,  there  are  three  facts  on  which  tiiere  is  more  or  less 
general  agreement:  (a)  the  fi-st  and  second  probability  density  functions  are  often  similir  to 
those  for  a  narrow-band  (laussian  noise  process;  (b)  the  width  of  the  power  spectrum,  or  equiva¬ 
lently,  the  inverse  of  the  correlation  time,  varies  approximately  linearly  with  the  Iransmitti.-d 
frequency:  (c)  if  one  assumes  that  tnese  fluctuations  are  due  primarily  to  the  relative  motion  of 
the  scattering  elements  on  the  siirfai-e.  the  average  relative  veloc  ity  of  these  eli-ments  is  of  the 
order  of  0.5  to  7.0  knots.  There  is  considerable  uncertainty  as  to  what  general  statements  can 
bi  made  about  the  precise  location  and  shape  of  the  spectrum  or  exactly  how  the  spectrum  vanes 
with  polari2.ation.  grazing  angle,  azimuth,  and  sea  condition.  There  appears  to  be  some  evi¬ 
dence  that  the  width  of  the  spectrum  is  independent  of  polarization  and  that  it  is  monotonically 
increasing  with  sea  roughnc’ss  and  grazing  angle*. 

With  regard  to  tiie  slow  flucttiatiuns.  about  all  that  can  be  said  iii  general  is  that  (a)  ton- 
sidcrable  energy  may  exist  in  the  region  0.001  ti>  lOeps  even  fo;  high  iiansmitted  frequencies; 
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Fig. 25.  incoherent  power  spectrum  os  a  function 
of  X.  (A(iopted  from  Goldstein. 


Coherent  power  spectrun  o;  o  function 
of  f.  X  «  3.2  cm,  p  =*  hh,  2*  ^  8*.  (After 

Kovoiy,  et  ^  .26) 
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Fig.  27.  Spectrum  v/idth  os  o  funcJion  of  X.  (See  text  for  detoils 
and  source  of  dafo.) 
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•r*. 

(b)  the  alow  fluctuations  often  appear  to  be  much  more  periodic  than  the  fast  fluctuations  and  the 
density  functions  often  appear  to  be  less  Gaussian  when  slow  fluctuations  are  included;  (c)  the 
slow  fluctuations  contain  relatively  more  energy  when  the  polarization  ia  horizontal  than  when  ^ 

it  is  vertical  and  when  the  clutter  is  spiky  than  when  it  is  continuous;  (d)  the  slow  fluctuations 
do  not  scale  linearly  with  frequency.  Frequently,  the  slow  fluctuations,  rather  than  being  re¬ 
garded  as  a  subject  for  study,  have  been  regarded  as  a  source  of  error  in  the  measurements  of 
the  fast  fluctuation  parameters.  Practically  no  attention  at  all  has  been  paid  to  the  statistics  of 

the  very  slow  fluctuations.  • 

Illustrations  of  the  pulse -to-pulse  fluctuation  rates  are  shown  in  Figs.  25  to  28.  (In  Figs.  25 

to  27,  it  has  beer,  assumed  that  the  fluctuations  are  due  to  Doppler  beats  and  the  frequency  scale 

is  plotted  in  terms  of  scatterer  velocities.)  Figure  25  shows  incoherent  power  spectra  for 

X  =  1.3,  3.2,  and  9.2  cm  plotted  as  a  function  of  v  =  fX/2  (f  =  observed  frequency  in  cycles/sec, 

X  =  transmitted  wavelength).  Since  these  spectra  are  incoherent  (i.e.,  the  returned  signal  has 

not  been  compared  to  any  reference  signal),  they  contribute  information  only  on  the  relative 

velocities  of  the  scatterers,  not  the  absolute  velocities.  Goldstein  found  these  spectra  to  be 

roughly  Gaussian  in  shape  and  to  be  independent  of  the  sea  condition.  Figure  26  shows  the  coher  - 

26  27 

eat  spectrum  for  X  =  3.2  cm.  According  to  Kovaiy,  et  al..  and  Hicks,  et  ^.,  the  core  spectrum 
(again  found  to  be  Gaussian)  was  obtained  for  all  wind  directions  and  calm  to  moderate  seas,  and 
the  core-plus-whitecap  spectrum  was  obtained  for  upwind  and  downwind  runs  on  rough  sea  states 
with  whitecaps.  The  half-power  width  of  the  total  spectrum  at  a  grazing  angle  of  4*  was  found  to 
vary  from  about  2  to  about  5  knots,  depending  upon  the  sea  condition.  There  was  some  evidence 
that  the  width  decreased  with  a  decrease  in  grazing  angle.  The  average  downwind  velocity  of  the 
scatterers  was  found  to  be  of  the  order  of  0  to  7  knots.  The  result  of  convolving  the  core  spectrum 
with  itself  to  obtain  the  equivalent  incoherent  spectrum  led  to  results  consistent  with  those  of  Gold¬ 
stein.  Figure  27  is  a  summary  of  estimates  of  the  root -mean -square  spectral  width  as  a  function 
of  the  transmitted  wavelength.  These  estimates  were  computed  by  Pike  from  pulse-to-pulse  cor¬ 
relation  data  obtained  at  the  M.  I.  T.  Radiation  Laboratory,  the  Naval  Research  Laboratory,  and 
the  M.I.T.  Lincoln  Laboratory .t  Of  the  80  data  samples  considered,  Pike  was  able  to  fit  about 
50  of  them  to  within  observational  errors  by  assuming  that  ^he  corresponding  power  spectrum 
was  a  mixture  of  the  following  components:  (a)  a  spike,  (b)  not  more  than  two  narrow-band 
Gaussian  spectra,  each  centered  cn  the  frequency  of  the  spike,  and  (c)  white  noise.  Of  these 
50,  about  30  could  be  fittec  by  white  roise  plus  a  single  Gaussian  component,  in  those  cases  in 
which  two  Gaussian  components  occurred,  the  energy  associated  with  the  broader  component  was 
almoBt  always  small  compared  to  that  of  the  narrower  component.  The  results  shown  in  Fig.  27 
refer  to  the  width  of  the  main  Gaussian  component  only.  To  what  extent  t.ie  white-noise  component 
resulted  from  receiver  noise  and  to  what  extent  from  the  clutter  signal  itself  was  not  determined. 

If  one  assumes  that  it  wa.s  all  receiver  noise,  then  the  chief  difference  between  the  widths  plotted 
and  the  widths  of  the  total  clutter  signal  occurs  in  the  elimination  of  the  spike  componem,  cxclu  ■ 

Sion  of  the  spike  necessarily  leading  to  greater  widths.  According  to  Pike's  an-^lysis,  this  spike 
contained  relatively  little  energy  except  at  136  cm,  At  this  wavelength,  the  energy  in  the  spike 

t'h*  peJnh  st  1.3,  3. 2,  onrf  9. 2  cm  ore  boMd  on  d-oto  ^»fcln«d  by  Ooid«f«in  end  oiiociatM  at  the  M,  L  T.  Rodio- 
Hon  the  points  ot  24  and  136  cm  ore  based  cn  date  obtoined  by  MoedoneSd  ond  associates  at  the  Na^-ai 

Research  Laboratory;  the  points  at  70  cm  ore  based  on  dotn  obtained  by  McGinn  ond  associates  at  Lincoln  Le^ 
oratory.  (The  Goidtfein  ond  Moedonald  doto  were  obt^ned  by  Rtke  and  McGinn  thr^^h  private  communication.) 
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Fig.  28.  Autocorrelation  function  for  long  time  tomples.  A  =  4.8  cm, 
0.5*<  a-<  4.CP,  lO'minute  somptes.  (After  Boring,  et  ol. 21) 


Fig.  29.  First  order  proJsabllHy  denrity  fonctlcn  on  tfie  power.  A  =  4,8  cm, 
a<  4.0^,  1 0*mlnut8  samples.  Abscisra  gives  sigrrat  power  In  arbitrary 
linear  units;  ordinote  gives  fraction  of  slgno!  received  in  unit  power  rortge. 
(After  Boring,  gt  oj.2^) 


fit 


exceeded  half  the  total  energy  in  the  function  for  about  one  third  of  the  samples.  The  data  at 
1.3,  3.2,  and  9.2  cm  are  for  vertical  polarization  and  a  grazing  angle  of  0.6*;  tlie  data  at  24  cm 
combine  both  vertical  and  horizontal  polarization  and  cover  grazing  angles  ranging  from  about 
5°  to  20';  the  data  at  70  cm  are  for  horizontal  polarization  and  !i  grazing  angle  of  about  i*;  and 
the  data  at  136  cm  are  for  horizonttil  polarization  and  grazing  angles  in  the  region  1*  to  14*.  As 
a  result  of  his  analysis,  Pike  has  concluded  that,  at  wavelengths  less  than  about  50  cm,  the 
width  of  the  main  Gaussian  component  increases  with  sea  roughness  but  is  independent  of  graz¬ 
ing  angle.  At  wavelengths  greater  than  about  50  cm,  he  has  concluded  that  this  width  increases 
with  an  increase  in  grazing  angltt,  but  is  independent  of  sea  roughnes.s.  (This  observtttlon,  to¬ 
gether  with  the  increase  in  the  width  at  about  50  cm,  has  led  Pike  to  conjecture  that  a  different 
scattering  mechanism  takes  over  at  this  wavelength.)  'Icmembering  that  the  spectrum  widths 
reported  by  Hicks,  et  (2  to  5  knots)  refer  to  the  total  width  at  the  half-power  points,  one 

notes  that  the  results  in  Fig.  27  at  3  cm  are  consistent  with  those  of  Hicks.  The  results  at 
136  cm  are  somewhat  different  from  those  reported  by  the  Naval  Research  Laboratory,  despite 
the  fact  that  both  sets  of  results  are  based'oh  the  same  deta.  Whereas  the  136-cm  results  in 

28 

Fig.  27  vary  from  about  1  to  9  knots,  the  root -mean -square  widths  reported  by  Ament,  el  aj_., 
are  confined  to  the  region  1  to  4  knots.  Also,  unlike  Pike.  Ament  found  the  spectral  width  to  be 
Independent  of  grazing  angle.  (In considering  these  differences,  it  should  be  noted  that,  whereas 
the  results  in  Fig.  27  refer  only  to  the  main  Gaussian  component,  the  results  of  Ament  refer  to 
all  the  components  mixed  together,  only  the  white  noise  being  eliminated.)  In  contradistinc¬ 
tion  to  the  results  reported  by  Hicks  at  3  cm,  Ament  found  the  width  at  136  cm  to  be  greater  for 
crosswind  than  for  upwind  or  downwind.  Figure  28  shows  the  appearance  of  the  slow  fluctua¬ 
tions  for  \  -  4.8cm  aiul  <r  in  the  region  0.5*  to  4.0*  when  long  observation  intervals  are  em¬ 
ployed.  The  upper  curve  is  an  expatuled  version  of  t)te  initial  drop  pictured  in  the  lower  curve 
and  corresponds  to  the  fast  fluctuations.  The  width  D  of  the  fast  component  was  found  to 
vary  inversely  with  frequency  and  wind  speed,  but  was  independent  of  polarization.  The  width 
!•:  of  the  slow  rantlom  component  was  found  to  be  independent  of  frequency,  polarization,  wave 
height,  and  wind  speed.  Tlie  nttio  13/A.  measuring  the  relative  power  in  the  two  components, 
was  found  to  be  ten  times  larger  on  horizontal  polarization  than  on  vertical  polarization  (0.7  on 
till  and  0.07  on  vv).  The  periodic  tail  component  (descrllied  by  the  parameters  G  and  C)  oc¬ 
curred  mainly  for  large  wave  tieights  and  high  wind  spt»ecl8.  l-’lgure  29(a-b)  shows  the  first 
probability  density  function  for  the  pulse-to-pulse  fluctur.i ions  in  power  as  a  function  of  polari¬ 
zation  for  X  4.8  cm  and  long  observation  intervals.  To  the  extent  that  the  process  is  Gaussian, 
it  can  be  shown  that  these  histograms  should  lie  linear.  The  histogram  for  vertical  polariza¬ 
tion  is  typical  of  the  results  usually  obtained  when  tlie  observation  interval  is  much  shorter  and 
the  alow  fluctuations  are  eliminated.  The  fact  that  the  graph  for  horizontal  polarization  has  a 
variable  slope,  whereas  the  graph  for  vertical  polarization  does  not,  is  consistent  with  the  fact 
that  the  slow  fluctuations  contain  relatively  more  energy  in  the  horizontal  case  and  Ihnt  the 
spikyness  effect  is  more  pronounced  in  the  horizontal  case.  Additional  data  on  both  the  first 
and  second  probability  density  functions  of  the  puIse-to-pulse  fluctuations  can  be  found  in  Kerr. 

.An  illustration  of  the  spikyness  effect  in  the  range  fluctuation  data  and  its  strong  dependence 
on  polarization  is  shown  in  Fig.  30(a-b).  These  spikes  were  found  to  be  correlated  with  the  oc¬ 
currence  of  well-defined,  steep-crested,  ocean  waves.  Some  Interesting  results  on  the  spiky- 

29 

ness  phenomenon  can  also  lie  found  in  a  paper  by  Macdonald. 
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Fig. 30,  Afomm^nc*  of  ipiky  ciutror  on  th«  A-'tsope.  X  =  4.S  cm,  b««Tnt>>dth  3>  2.3*, 
puiM^ongth  =  0.2 p«ec,  antonno  olevrtJoo  =  52  ft;  (a)  p  =  w,  pic^tr^ tokon  nitandan, 
fr)  p  =  W»,  pcturis  takon  only  whon  w  ocHo  wos  vJjjbJe.  (Afi»r  Boring,  «t  of .2*) 


B.  Concepts  for  Interpreting  Experimental  Data  on  Sea  Clutter 

1.  Hanciom-Sciittorcr  ('onctpt 

One  idea  fri-qui-nlly  used  in  inlerpreling  the  data  on  the  clutter  signal  is  that  the  source  of 
the  cluller  can  he  regarded  as  a  v.ollection  of  a  large  numher  of  random,  similar,  ind(>pendcnt , 
uniformly  distributed,  point  scatterers.  This  is  the  standard  assumption  used  in  phenomenologi¬ 
cal  models  for  deSitribing  targets  that  are  very  complex  and  which  contain  no  dominant  reflector. 
Aside  from  greatly  simplifying  certain  computations,  this  assumption  leads  to  the  predictions 
that  (a)  the  probability  density  functions  are  Gaussian:  (b)  the  average  power  varies  linearly 
witli  the  itrea  of  illumination:  (c)  the  autocorrelation  fum-tion  of  the  fluctuations  in  range  is  eqinil 
to  the  autoeorridation  function  of  the  transmitted  pulse.  Together  with  the  assumption  of  a  narrow¬ 
band  transmitted  signal,  it  allows  one  to  treat  tin-  clutter  signal  as  narrow-band  Giiussian  noise. 
This  concept  will  obviously  be  inappropriate  for  describing  the  clutter  when  the  radar  has  suf¬ 
ficient  resolution  to  resolve  the  irreguhi cities  on  the  o«'«*an  surfai-e  or  whi-n  these  irregiihirities 
are  correlated  over  distances  which  are  appreciable  with  respect  to  the  dimensions  of  tlu-  illu¬ 
minated  area,  l-'or  discussions  tif  clutt«-r  signals  as  narrow-band  Giiussitm  noise  processi>s.  the 

reader  is  refi-rred  to  Kerr,^"^  l.awson  tuid  L'hlenbei-k.^^  Kelly  and  l.erner.^^  and  McGinn  anti 
32 

Vikt .  For  :i  general  background  in  relevant  statistical  theory,  the  reader  is  referred  to 
Cram^r.^^  IJavenptirt  tintl  Hoot.^^  tint!  Hendat 


2.  Uoppler-lmiige  Concept 

A  concept  which  frequently  has  been  used  for  interpreting  tin*  dtitii  on  pulse-to-pulse  fluctua¬ 
tion  rates  (already  employed  in  presenting  the  data  In  Sec.  Ill-A)  is  that  the  pulse-to-pulse  fre¬ 
quency  spectrum  of  the  clutter  can  be  interpreted  as  the  Doppler  image  of  the  velocity  spectrum 
of  the  scatteri-rs.  If  u  denotes  the  velocity  spectrum  of  th»*  sciitterers.  f  the  fluctuation  rate 
(in  cycles  per  second)  of  the  clutter  signal,  and  X  the  transmitted  wavelength,  the  fluctuation 
spei'trum  ('(D  is  given.  ac<-ording  to  this  assumption,  by  I  (f)  -  u(fX.  2).  The  fact  that  the  spec¬ 
trum  of  the  fast  fhictuat itms  scales  linearly  with  tlu*  transmitted  frequency  at  the  higher  fre¬ 
quencies  woidd  appear  to  lend  strong  support  to  th  s  interpretation  for  the  fast  fluctuations  tit 
these  frequencies.  However,  in  view  of  the  fact  that  (a)  at  small  grazing  angles,  the  irregultiri- 
ties  on  the  ocean  surface  which  backscatter  most  effectively  at  one  wavelength  are  not  likely  to 
constitute  the  most  effective  scatterers  for  other  wavelengths,^  and  (b)  different  irregularities 
are  likely  to  travel  with  different  velocities,  this  confirmation  of  the  Doppler-image  interpre¬ 
tation  must  be  regarded  with  considerable  skepticism.  In  general,  although  the  Doppler-image 
inierpreialion  of  the  fast  fluctuation  spectrum  may  be  correct  for  a  wide  range  of  transmitted 
frequencies,  one  has  no  right  to  assume  that  this  spectrum  will  scale  linearly  with  the  trans¬ 
mitted  frequency  unless  one  can  show  that  the  important  scattering  elements  at  the  different 
frequencies  travel  with  the  same  velocity.  In  this  connection,  it  is  important  to  note  that  the 
Doppler-image  coni’ept  has  frequently  led  (along  with  other  considerations)  to  tlie  belief  that. 

at  the  higher  frequencies,  the  scattering-ittotions  are  similar  to  the  orbital  motions  suffered  by 
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a  small  particle  of  water  (see,  for  example,  Hicks,  el  a).  )  and,  at  the  lower  frequencies,  these 


tAt  imall  grazing  ongUs,  whar*  no  portion  of  tho  lurfoco  It  normol  to  tho  incoming  onorgy,  Irrogulorltlot  which 
or*  largo  with  rotpoct  to  tho  tronsmlttod  wavolongth  will  rofloct  tho  ono^y  forward.  On  tho  othor  hand,  Irrogu- 
larltlot  which  oro  small  with  rotpoct  to  tho  tronsmlttod  wovolongth  will  couso  littio  tcattorlng  In  any  dlroctlon. 
Thus  tho  Irrogulorltios  which  will  bo  most  offoctivo  ot  smoll  grazing  onglot  (por  Irrogularlty)  will  bo  thoio  whoso 
tizo  It  of  tho  ordor  of  tho  tronsmlttod  wovolongth. 
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motions  are  similar  to  the  phase  velocities  of  the  ocean  waves  (see.  for  example.  Ament,  et  al.  ). 

At  extremely  low  frequencies  (li  to  tSMcps),  the  identification  of  the  clutter  spectrum  with  the 
phase  velocities  has  received  strong  support  from  the  work  of  Crombie^^  and  Stutt,  Fricker, 
Ingalls  and  Stone^^  (see  diffraction-grating  concept  in  Sec.  IlI-B-5.) 

3.  Interference  Concept  — 

Another  concept  which  frequently  has  been  used  to  interpret  certain  aspects  of  the  clutter 
data  is  that  the  rougti  ocean  surface  can  be  regarded  as  a  smooth  ocean  surface  with  scatterers 
elevated  above  it,  Insofar  as  the  return  from  each  scatterer  can  then  be  thought  of  as  consist¬ 
ing  of  a  direct  and  reflected  ray  as  discussed  in  Sec.  1.  this  concept  has  proved  useful  in  inter¬ 
preting  the  results  on  the  average  power  at  small  grazing  angles,  as  well  as  the  occurrence  of 
spikyness  at  small  grazing  angles.  According  to  this  concept,  the  variation  of  with  a, 
and  p  is  interpreted  in  terms  of  'he  variation  of  the  interference  pattern  with  these  variables, 
and  the  variation  of  with  t  is  interpreted  in  terms  of  the  he  ight  distribution  of  the  scatterers. 
Thus,  fur  example,  assuming  that  the  pularization  is  horizontal  (so  that  the  reflection  coefficient 
r  1).  the  critical  angle  o  j  below  which  falls  off  rapidly  with  decreasing  a  is  interpreted 
as  the  angle  at  which  tin-  bottom  lobe  in  the  interfer«-nce  pattern  becomi'S  too  far  removed  from 
the  surface  to  be  avertiged  out  by  the  height  distribution  of  the  scatterers.  The  fact  that 
increases  with  X  then  follows  directly  from  the  fact  that  the  height  of  the  bottom  lobe  increases 
with  X.  The  differences  between  horizontal  and  vertical  polarization  are  interpreted  in  terms 
of  the  differences  in  the  interference  patterns  implied  by  the  difference's  in  the  reflection  co¬ 
efficients  for  the  two  pol.irizations.  In  general,  the  return  on  vertical  polarization  tends  to  be 
stronger  because  the  phase  tag  of  the  reflection  coefficient  is  such  that  more  energy  exists  at 
points  close  to  the  surfiice.  As  X  decreases  or  the  scatterers  are  raised  higher  above  the  sur¬ 
face  (corresponding  to  tin  increase  in  sea  roughness),  more  of  the  interference  pattern  tends  to 
be  averaged  out  and  its  value  tit  points  close  to  the  surface  becomes  less  crucial.  For  very 
small  X  or  very  lanigh  setis,  the  phase  of  the  interference  pattern  (determined  by  the  phase  of 
the  reflection  coefficient)  lu-comes  Irrelevant  anti  dep-nds  only  on  the  average  field  strength 
(iletcrmined  by  the  mtignitude  of  the  reflection  coefficient).  I'intilly.  the  spikyness  effect  which 
ticcurs  at  small  grazing  angles  (where  the  area  of  illumination  is  actually  thi‘  largest)  Is  t  x- 
plained  by  the  fact  that  when  the  bottom  lobe  is  sufficiently  fur  removed  from  the  surface,  only 
those  scatterers  with  exceptional  ht.'ight  will  be  illuminated.^ 

In  lioMstein's  original  discussion'''  of  ptjsslble  theories  for  explaining  the  behavior  of 
he  considered  the  interference  concept  only  in  connection  with  the  assumption  that  the  scattering 
elements  were  droplets  of  water  actually  septirated  from  the  water  surface.  This  theory  (the 
so-called  "droplet  theory")  proved  itseful  for  explaining  the  polarization  and  grazing  angle  de¬ 
pendence  offfjj,  but  failed  when  the  frequency  dependence  was  considered  unless  it  was  assumeil 
that  the  drops  were  comparable  in  size  to  the  transmitted  wavelength.  (Goldstein  also  pointed 
out  that  the  polarization  dependence  was  greatest  in  calm  weather  when  the  drops  were  least 
likely  to  occur.)  In  the  opinion  of  many  workers  in  this  field,  it  is  possible  to  retain  the  forw.'ird- 
scattering  postulate  without  assuming  that  the  scattering  elements  are  drops.  There  is  no  reason 

tTo  the  txttnt  that  tht  rtductlon  In  mo  clutter  achieved  by  using  horizontal  polarization  and  large  X  con  be  ex¬ 
plained  by  the  interference  concept,  these  clutter-rejection  techniques  should  be  regarded  as  special  coses  of  the 
ETI  concept  discussed  In  Sec.  Il-D.  According  tc  the  Interference  concept,  these  techniques  are  effective  in  re¬ 
jecting  the  clutter  because  they  txill-out  targets  t  very  srtxsll  altitudes. 


to  believe  that  forward  scat'ering  does  not  play  an  important  lole  even  when  the  scattering  ele¬ 
ments  are  irregularities  on  the  surface  itself.  A  more  serious  problem,  perhaps,  independent 
of  tl;c  scatterers'  identity,  is  the  assumption  that  the  forward-scafered  energy  can  be  represented 
in  the  same  way  as  for  a  surface  which  is  smooth.  This  is  obviously  an  extremely  artificial  a.s- 
sumption  and.  at  best,  represents  only  a  very  crude  approximation.  (In  defense  of  this  approxi¬ 
mation,  however,  it  should  be  noted  that  in  most  of  the  models  which  constitute  serious  attempts 
at  constructing  real  physical  theories,  the  question  of  multiple  reflect. ons  is  ignored  entirely.) 

In  order  to  illustrate  the  application  of  this  concept  to  the  experimental  data,  assume  that 
(a)  the  scatterers  are  staf tstically  independent  and  statistically  identical;  (b)  the  scatterer-s  are 
distributed  uniformly  over  the  illuminated  area  with  a  density  »):  (c )  the  earth-inodification  fai  - 
tor  i;  and  the  cross  section  a  of  a  single  scatierc-r  are  independent.  Denoting  the  ensemble 
average  by  a  bar  and  making  use  of  Kq.  (9)  for  M,  one  can  show  that 
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where 
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ff  1 


|h'i^  1  1  4  K  exp(-i^-^)r' 

=  1  +  +  4iK!“  +  2Ke  fK^  exp(-^iu,'^)  4  2K(1  4  jKi^)  exp(-i^'^)l 

Assuming  fur<hermore  that  (d)  K  and  ^  can  be  approximated  by  K  I"  and  ^  2ol!  (II  being 

the  height  of  the  sc -itteser).  one  can  rewrite  <he  equation  for  ‘K'"  as 

=  14  Irl"^  4  4|r|^  4  2Re  (r^  cxp(-4U  )  4  ani  4  irj^)  e.xp(-2ini  o?) 

where  4  -  uiall/c  ZiroU  \  and  the  averaging  operation  applies  to  the  variable-  M.  .-Xssumc-  now 
that  the-  set  of  probability  density  functions  (p^dU)^  v»n  II  corresponding  to  the-  set  of  sea  con¬ 
ditions  f  is  a  one-parameicT  family  in  which  each  member  p^  is  detet  mined  by  the  mean  \;ilne 
H.  and  the  higher  order  moments  n'*  are  related  to  H  by  ecjuations  of  the-  form  m''  *^ii**' 

ind«-pe'(d<‘iit  of  T-  Then  |K|^  will  depend  on  A.  o .  and  r  only  insofar  .is  they  detc-rmine  the 
prc'duc’  i  Ziroll/X.^  This  depi-ndence  can  be-  described  exphciily  by  expanding  expt  -lifl  and 
exp(— 2i^)  in  power  series,  averaging  each  term  with  respc-ci  to  II.  and  tlien  replacing  ii*'  by 

H^ll*'.  .Assuming  that  C  is  smtil!  and  that  the  sc-aiterers  are  c  lose  to  the  sui  luce  with  respei  t 

to  Viiriations  in  the  inierfc-reitee  pattern,  one  can  approximate-  ;  1:;“  by  the-  fir.'-.i  fc-w  terms  ol 
this  sc-ries  doing  to  the  other  extreme  and  assuming  tnal  i  is  large  .ind  that  the-  se.itic-rc-r 
h»-ighls  covc-r  a  number  of  lobes,  one  oblams  iK|  1  ■*  .  T  *4  I''*'  .-Xssuming  lii.ii  I'  1 

(horizontal  polarization',  one-  has 


tNot#  fhor  ‘he  variable  f  h  nothing  Ksore  than  the  clowical  Koyleigh  roughnew  porometer  often  used  for  decid- 
it^  whether  or  not  o  surface  is  "rough.*’ 
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6  +  2  cos  44—8  cos  24 


lEl"  =  6 


v.'hen  4  is  large 


|e1^  =  16  8^4**  when  4  is  small 


lEl"  = 
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(1  +  44“)  (1  +  164*') 


when  pj(H)  •  4r  j 


I  E|^  =  2  (3  -  4  exp(-  nj^)  +  exp(-4ir4‘^)j 
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when  p.<H)  ^  —  exp(--^) 
‘  »n  '  wH  ' 


0  H  » 
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Returning  to  a  general  T  for  use  in  vertical  polarization  and  assuming  that  p^(H)  -  (l/H)  expl-H/M), 
0  <  H  <  <*>.  one  obtains 


- T  .  ,  2|  r  Hcos  2«>_  +  44  sin2«>„) 

iEr 1  +  irr  +  4|rr  +  — - e — r - ^ 

1  +  164*^ 

4|ri(l  +  Irp)  (cos  (Pp  +  24  sin«>p) 


1  +  44^ 


(51) 


If  attention  is  restricted  to  horizontal  polarization  and  the  variation  of  or^  with  X.  a,  and  r  is 
ignored,  these  results  imply  that  cr  is  a  function  of  4  -  2»ali/X.  For  small  4  .  one  olituins 

-4  _  °  _o 

(Tq  «  4  •  For  large  4  .  one  has  •  These  results  are  substantially  consistent  with  the 

data  except  for  (a)  the  behavior  of  above  the  critical  angle  a ^  and  (b)  the  behavior  of  with 

X  and  C  above  the  critical  angle  a  j  before  "saturation."  In  addition  to  these  two  major  dis¬ 
crepancies.  the  first  of  which  is  important  only  if  one  is  interested  in  grazing  angles  a  >  20*. 
one  should  also  note  that  (c)  whereas  these  results  predict  that  «  (x/H)”.  n  1.  the  data 
indicate  that  n  may  be  anywhere  between  0  and  1,  and  (d)  whereas  these  results  predict  that 
<T^  »  n  =  0,  for  CT  >  a  j.  the  data  indicate  that  n  may  be  anywhere  between  0  and  2.  In  gen¬ 
eral,  there  are  two  attitudes  which  one  may  assume  toward  these  discrepancies:  one  can  con¬ 
clude  that  these  results  are  inappropriate  or  one  can  regard  these  discrepancies  as  constituting 
data  on  aV  (hi  the  whole,  tho.se  who  have  been  concerned  with  this  concept  have  assumed  th  it 
the  discrepancies  constitute  data  on  oV 

In  the  discussion  of  the  experimental  data,  (r^in)  for  0  <  was  characterized  by  the  plateau 
level  and  the  critical  angle  .^n  equivalent  characterization  can  be  effected  in  terms  of 
the  interference  concept  by  use  of  the  parameters  a\  and  H.  In  general,  the  most  effective  way 
of  characterizing  a  ^^(o)  curve  in  terms  of  these  parameters  is  to  overlay  the  1  E|  curve  and 
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to  evaluate  and  H  by  slidinj  the  1  curve  along  the  abscissa  and  ordinate  until  one  obtains 
a  oest  fit.  The  results  of  applying  this  procedure  fusing  Eq.  tsO)}  to  the  uatu  in  Fig.  to  aTe  shown 
on  Fig,  18  by  the  smooth  curves.  (These  fits  are  bt*Uer  than  average..)  Similar  fits  (using  Eq.  (49)) 
were  obtained  by  Goldstein.  A  slightly  different  technique,  but  one  which  constitutes  an  essen¬ 
tially  equivalent  test  of  the  interference  mechanism,  has  been  employed  by  Katzin.^^  The  values  of 
and  H  obUoined  from  the  fits  shown  in  Fig.  IS  are  given  by  H(X  -  70cm)-  4.9ft,  H(X  -  10  cm)  - 
2.4ft,  {T^(\  =  70cm)  =  -52db,  and  o'(\  -  10cm)  =  -44 db.  Inasmuch  as  (El  ^  -  6  for  large  one 
finds  that  o'  is  related  to  a  by  a'  =  o  /6.  The  relation  of  H  to  qi  will  deoend  on  the  precise 
definition  of  a^.  Defining  the  normalized  critical  angle  ( ^  in  the  1  E|  ^  curve  by  {  ^  -  0.5,  one  ob- 
tains  =  0.0b X/H  (corresponding  in  Fig.  18  tologa^  -  —1.4  for  the  70-cm  cur\o  and  iogOj  =  -1.9 
for  the  10-cm  curve).  In  general,  in  the  writer’s  experience,  the  values  of  Tl  obtained  from 
this  fitting  process  are  related  to  the  reported  wave  height  H  by  the  inequality  0.2  <  fi/H  ^  2. 
Inasmuch  as  the  observed  exponents  in  the  relations  a  and  “  (1/Tl)”  are  sometimes 
less  than  unity,  H  will  sometimes  exhibit  a  small  dependence  on  X.  increasing  wuth  an  increase 
in  X,  and  the  '■arlat'.on  of  H  with  H  will  sometimes  be  less  than  linear.  .According  to  the  model, 
one  would  expect  that,  in  general.  11  would  ’tic  correlated  with  the  distribution  of  the  heights  of 
the  large  waves  and  the  characteristics  of  the  macrostructure,  whereas  a'^  would  be  correlated 
with  the  state  of  the  fine  structure  superimposed  oti  these  waves  and  with  the  momentar’.  local 
wind  conditions. 

The  theoretical  curve  in  Fig  22,  originally  computed  by  Goldstein,  is  obtained  f*om  the 
above  model  by  assuming  that  <7^  is  independent  of  polarization  and  computing  the  ratio  cr^ivv),' 
<T^^(hh)  =  !  E{vv)  [2/j  K(hh)  j  ^  as  a  function  of  I!  {using  Eq.  {51)|.  Considering  the  scatter  in  the 
data,  the  fi*  is  reasonably  good.  .Although  the  values  of  H  required  to  fit  these  data  appear  un¬ 
reasonably  small  in  comparison  to  the  values  of  H  obtained  feom  the  cri’ical  angle  data,  it 
should  be  noted  that  all  values  of  the  theoretical  pfdarization  ratio  are  approximately  the  same 
foi  n  >0,5  ft.  Unfortunately,  if  one  attempts  to  make  the  same  sort  of  comparison  using 
Macdonald’s  data,  the  results  are  not  nearly  so  satisfying.  SjMJcifically,  one  finds  that  t  le  val¬ 
ues  of  H  reqiured  to  fit  the  polarization  ratio  data  (a)  increase  with  a  decrease  in  grazing  angle. 

(b)  increase  with  an  increase  in  wavelength,  and  (c)  are,  without  question,  considerably  smaller 
than  those  required  to  fit  the  critical  angle  data  on  horizontal  cr  vertical  polarization  laKen  alone. 
(The  increase  of  H  with  X  obtained  from  these  data  appears  to  occur  at  approximately  the  same 
rate  as  the  increase  of  H  with  X  obtained  from  the  critical  angle  da: a.)  The  problem  implied  by 
the  -flative  smallness  of  H  and  the  increase  of  H  with  a  decrease  in  grazing  angle  can  be  re¬ 
stated  as  follows;  if  the  values  of  above  the  critical  angle  are  the  result  uf  the  scattcrers 
being  distriby;t=d  ovi-r  a  whole  lobe  of  the  imerference  pattern  so  that  the  interference  effect  is 
averaged  out.  then  why  .are  ihc  results  on  vertical  polarization  still  greater  than  the  results  on 
horizontal  polarization  when  the  grazing  angle  is  larger  than  the  critical  angles  for  both  polari¬ 
zations?  if  there  is  any  difference  at  all  between  the  two  pclarizatlons  above  the  critical  angler, 
the  results  on  vertical  polarization  should  be  smaller.  As  far  as  the  writer  knows,  these  limi¬ 
tations  in  the  interference  concept  for  explaining  the  polarizatioit-ratio  data  were  first  pointed 
out  by  Macdonald. 

One  further  point  concerning  the  interference  concept  involves  its  relation  to  tne  pulse-to- 
pulse  fluctuation  data.  According  to  the  interference  assumption,  one  component  in  these  fluctu¬ 
ations  should  arise  from  the  Doppler  beats  beiwcen  the  direct  and  indirect  rays.  Inasmuch  as 
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the  rates  of  these  beats  will  he  proportional  to  a.  these  rates  shculii  be  slow  for  smaii  a  and 
increase  as  a  increases.  Also,  to  the  extent  that  the  identity  of  the  scatterers  is  independent 
of  frequency,  one  would  expect  these  rates  to  scale  linearly  with  frequency.  The  results  by 
Boring,  al.,^^  which  indicate  that  the  width  of  the  slow  fluctuation  component  is  independent 
of  frequency  at  small  grazing  angles  implies  that  this  source  of  fluctuations  is  relatively  un- 
impoi  tant. 

4.  Facet  Concept 

In  the  three  concepts  mentioned  above,  no  attention  was  paid  to  the  problem  of  specifying 
the  detailed  reflective  properties  of  the  scatterers.  Thus,  the  above  concepts  are,  at  best, 
only  capable  of  providing  a  framework  for  interpreting  the  data.  Some  work  on  which  takes 
account  of  the  interference  concept  and  goes  on  to  consider  a  specific  model  for  the  scattering 
elements  themselves  is  contained  in  the  work  of  Katzin.  Assuming  that  the  ocean  can  be  re¬ 
garded  as  a  collection  of  flat  plates  (or  "facets")  of  varying  izes  and  slopes,  Katzin  draws  the 
following  conclusions.  At  small  grazing  angles,  where  nont  f  the  facets  are  normal  to  the 
incoming  wave,  the  facets  which  backscatter  most  effectively  a.e  those  having  a  perimeter  of 
about  X/2  and  the  backscattcring  of  a  facet  increases  ubout  as  the  square  of  its  slope  (thus  mak¬ 
ing  the  wave  crests  of  special  importance).  At  large  grazing  angles,  the  facets  which  back¬ 
scatter  most  effectively  are  those  which  are  normal  to  the  incoming  wave  and  which  are  large, 
the  angular  variation  of  being  determined  primarily  by  the  slope  distribution  of  the  facets. 
When  the  grazing  angle  is  small,  is  proportional  to  wind  speed,  but  when  the  grazing  angle 
is  large,  is  inversely  proportional  to  wind  speed.  At  small  grazing  angles,  the  frequency 
dependence  of  <r^  is  determined  by  the  size  distribution  of  the  facets.  As  far  as  the  writer  knows, 
Katzin  has  examined  the  implications  of  this  model  only  for  the  average  power  in  the  clutter  sig¬ 
nal.  Some  results  relating  the  facet  concept  to  irregularities  actually  occurring  on  water  sur- 
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faces  can  be  found  in  the  papers  by  Schooley. 

5.  Diffract  ion -Grating  Concept 

One  further  concept,  based  on  the  Doppler-image  notion,  which  has  proved  extremely  useful 
for  interpreting  the  frequency  shift  in  the  coherent  clutter  spectrum  at  extremely  low  frequencies 
(13  to  25Mcps)is  the  diffraction-grating  concept  first  applied  by  Crombie**^  and  later  applied 
by  Stutt,  et  al.^^  Assuming  that  the  wave  trains  which  will  contribute  most  strongly  to  the  clut¬ 
ter  signal  are  those  which  arr  traveling  toward  or  away  from  the  radar  and  have  a  spacing  L 
equal  to  one-half  the  electromagnetic  wavelength  X,  they  have  been  able  to  predict  the  Doppler 
shift  in  the  returned  signal  to  a  high  degree  of  precision  by  use  of  the  classical  hydrodynamic 

I  /y 

formula  v  =  (gL/2jr)*'  *  relating  the  velocity  v  of  a  gravity  wave  to  its  wavelength  !.  (g  being 
the  acceleration  due  to  gravity).  According  to  this  model,  the  Doppler  beat  frequency  between 
the  transmitted  and  received  signals  is  given  by  f  =  2v/X  =  (g/»X)*/^.  More  recently.  Ranzi^^ 
has  shown  that  this  concept  is  applicable  to  frequencies  as  high  as  4i5Mcps,  but  as  one  would 
expect,  is  not  applicable  at  a  frequency  of  6800  Meps. 

C.  Fonnuti  Seatitring  from  the  Ocean 

According  to  the  previous  discussion,  the  radar  designer  needs  to  concern  himself  with  the 
effects  of  surfact  roughness  on  forward  scattering  for  two  reasons;  (i)  in  order  to  be  able  to 
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Incorporate  theae  effects  into  t?se  earth>modifle&tloA  E  needed  for  describing  the  signal 

returned  from  a  target,  and  (2)  in  order  tn  he  able  to  inocrporate  these  effects  into  varioui 
models  used  for  describing  the  clutter  signal  (such  as  the  "interference  concept"  described  in 
Sec.  lU-B'S).  Inasmuch  as  the  height  of  a  clutter  backscatterer  above  the  mean  surface  is  or¬ 
ders  of  magnitude  smaller  than  the  height  of  an  elevated  aircraft  target,  the  effects  in  the  twc 
cases  may  be  radically  different.  Unfortunately,  to  the  writer's  knowledge,  there  are  no  data 
available  which  can  be  used  with  confidence  in  the  clutter  models.  Aside  from  the  usual  dif¬ 
ficulties  encountered  in  scattering  experiments,  in  this  case  one  is  faced  with  the  additional 
problem  of  having  to  explore  the  scattered  field  at  heists  above  the  surface  which  are  of  the 
same  order  of  magnitude  as  the  heights  of  the  irregularities  themselves. 

As  far  as  the  effects  of  surface  roughness  on  the  target  signal  are  concerned,  although  a 
certain  amount  of  relevant  data  exists,  much  of  it  suffers  from  the  same  sorts  of  inadequacies 
a'j  one  finds  in  the  clutter  data  (e.g.,  inadequate  specification  of  the  sea  surface).  In  addition, 
since  most  of  the  forward-scattering  experiments  employ  only  one-way  pr-opagation,  to  the  ex¬ 
tent  that  one  is  concerned  with  the  radar  problem  and  not  the  communication  problem,  one  is 
faced  with  the  further  difficulty  of  having  to  convert  the  statistics  of  the  one-way  signal  to  those 
of  the  two-way  signal.  Although,  theoretically,  all  the  information  which  is  needed  for  describ¬ 
ing  the  two-way  signal  can  be  obtained  from  the  statistics  of  the  one-way  signal,  in  actual  prac¬ 
tice.  the  amount  of  information  available  on  the  one-way  signal  is  often  insufficient  for  this  pur¬ 
pose.  For  example,  even  if  one  assumes  that  the  polarisation  of  the  antenna  is  pure -horizontal 
or  pure-vertical,  since  the  slope  of  the  ocean  surface  is  variable,  in  order  to  obtain  sufficient 
information  through  a  one-way  experiment  to  describe  the  two-way  signal,  one  will  have  to  take 
data  on  all  four  polarization  combinations  hh,  w.  hv,  and  vh.^  In  addition  to  this  need  for  in¬ 
creased  polarisation  information,  one  in  also  faced  with  a  need  for  increased  statistical  informa¬ 
tion  in  the  form  of  higher  order  moments.  For  example,  in  order  to  obtain  the  average  value  of 
the  target  signal,  since  the  II  path  involves  two  reflections,  data  are  required  on  the  average 
value  of  the  square  of  the  reflection  coefficient  as  well  as  on  the  average  value  of  the  reflection 
coefficient  itself. 

Aside  from  the  problems  encountered  in  trying  to  determihe  the  two-way  signal  from  the 
results  of  one-way  experiments,  the  forward-scattering  situation  differs  from  the  clutter  situa¬ 
tion  in  that,  whereas  the  clutter  problem  can  be  treated  as  a  random  noise  problem  (see  the 
"random-scatterer  concept"  (Sec.IIl-B-l)I.  the  forward -scattering  problem  must  be  treated  as 
a  signal-plus-noise  problem.  Thus,  in  the  forward-scattering  case,  the  statistics  tend  to  be¬ 
come  more  complicated. 

Until  the  last  few  years,  the  usual  means  for  describing  the  effects  of  roughness  on  the 
forward -scattered  energy  was  to  assume  that  this  rouf^ess  could  be  accounted  for  by  replac¬ 
ing  the  smooth  surface  reflection  coefficient  F  by  an  "effective  reflection  coefficient"  P  and 
evaluating  i  P  |  by  measuring  the  amplitude  of  adjacent  maxima  and  minima  in  the  interference 
pattern.  Defining  K  by  K  =  R  J)rgp(I)g^(I)/(R  +  A)g^(D)g^(D).  and  R  by  R  =  (?/r)K.  one  can 
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cv?.lu6t6  of  tlis  G!t0“vs?2'^  irots^is  s£t?tK^n*o^ifiCwtion fwCtc**  S  ^  1  K 

exp{-l«Vc)by  the  equations  i^l  =<i£l^ax  ~  !®lmin>/<l®imax  ^  1^1  =  irnRi/lKl. 

The  results  obtained  in  this  manner  show  the  following  characteristics.  For  very  smooth  seas, 
the  results  on  both  horizontal  and  vertical  polarization  follow  the  smooth  surface  theory  tc 
within  experimental  error.  As  the  sea  becomes  rougher,  however,  whereas  the  results  on 
vertical  polarization  still  cluster  around  the  theoretical  curve  for  a  smooth  surface  and  evidence 
only  a  mild  degree  of  scatter,  the  results  on  horizontal  polarization  fall  below  the  theoretical 
curve  U.e.,  below  unity)  and  exhibit  a  high  degree  of  scattei',  A  typical  result  illustrating  these 
characteristics  is  shown  in  Fig.  31(a>b).  For  a  summary  of  this  earlier  work,  the  reader  is  re¬ 
ferred  to  Kerr.*^  Besults  similar  to  those  shown  in  Fig.  31  have  been  obtained  more  recently  by 
Macdonald  for  a  wavelength  of  24  cm  (unpublished  data). 

In  the  last  few  years,  substantial  advances  have  been  made  in  this  field  through  the  efforts 
of  the  Applied  Physics  Laboratory  (APL)  at  Johns  Hopkins  University  and  the  Electrical  Engineer¬ 
ing  Research  Laboratory  (EERL)  at  the  University  o?  Texar.  This  work  has  been  well  summa¬ 
rized  by  Beard^®  and  will  be  discussed  here  cmly  briefly. 

Let 

Yq  3  complex  modulation  of  signal  received  along  direct  path 
Yj  B  complete  modulation  of  forward-scattered  signal 
Yt  =  Yjj  +  Yj  =  total  received  signal 
Yj  =  average  value  of  Yj 

n-  t  I  r  I  m 

29  =  I  Yj  -  T?ji  a  power  in  the  fluctuations 
8  a  root-mecn-square  wave  height 
#  s  grazing  angle 
X  3  transmitted  wavelength. 

The  model  used  by  APL  and  EERL  for  processing  their  data  assumes  that  the  components  of  the 
fluctuating  signal  Yj  -  Yj  are  statistically  independent  random  variables  with  Gaussian  density 
functions  of  mean  zero  and  variance  9^.  (In  other  words,  the  statistics  of  the  forward-scattered 
signal  are  assumed  to  be  equivalent  to  those  which  would  arise  if  one  took  the  backscattered  sig¬ 
nal  as  described  by  the  random-scatter  concept  mentioned  in  Sec.IR-B-l  added  a  sine  wave.) 
The  steady  signal  Yj  is  referred  to  as  the  coherent  component  and  the  fluctuating  signal  Yj  —  Yj 
(assumed,  like  the  clutter  signal,  to  arise  from  a  large  number  of  independent  random  scatterers) 
is  referred  to  as  the  incoherent  component. 

The  results  (Gained  on  the  normalized  cc^erent  and  incoherent  terms  lYrl/iY^ri  and 
lYj^r  I  are  shown  in  Figs.  32  and  33.  respectively.  These  graphs  cc^ftaln  data  for  both  horizontal 
and  vertical  polarization  and  for  wavelengths  of  5.3,  3.2  and  0.9  cm.  Some  of  these  data  were 
obtained  by  measuring  the  total  field  Y.^.  and  extracting  information  on  Yj  throu]0  the  use  of  the 
maxima  and  minima  in  the  interference  pattern,  and  some  by  measuring  the  forward- ncsttered 
field  Yj  directly  through  the  use  of  antenna  patterns  with  a  high  degree  vertical  resoluti<m. 

The  distance  between  the  transmitting  and  receiving  antennae  was  5425  ft  and  the  range  of  values 
of  8  encountered  (obtained  from  step-gauge  recordings)  was  0.t2  to  0,?3ft.  The  curve  in  Pig.  32 
is  the  "classical"  coherent  theoretical  curve 

y  »  exp(-8t^^)^Vx^) 

Thie  curve  is  based  or  the  assumption  of  &  Gaussian  distributicri  for  the  hei|^t  of  the  water  surface 
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and  has  been  derived  by  many  writers  (some  of  whom  arc  referred  to  Ui  the  paper  by  Beckmann*). 
It  takes  account  merely  of  the  phase  di^rences  in  the  indirect  rays  caused  by  variations  in  the 
height  of  the  surface.  (The  more  general  results  which  have  been  obtained  frequently  include 
this  result  as  a  multiplicative  factor.)  Assuming  that  the  data  on  the  coherent  term  are  obtained 
from  a  total  field  measurement,  one  sees  that  these  data  differ  from  those  shown  in  Fig.  31  as 
follows:  whereas  in  Fig.  31  the  urta  were  obtained  by  measuring  |  Yj^  +  and  |Yjj  + 

the  data  in  Fig.  32  were  obtained  by  measuring  |Yjj  +  and  |Yp  +  Apparently, 

when  the  processing  iti  done  according  to  the  latter  procedure,  the  difference  in  the  effects  of 
roughness  on  the  two  polarizations  tends  to  disappear. 

One  should  also  note  the  approximate  equivalence  of  the  following  three  angles:  (1)  the 
angle  at  which  the  theoretical  coherent  curve  begins  to  fall  below  the  coherent  data;  (2)  the  an¬ 
gle  at  which  the  experimental  incoherent  curve  reaches  its  m;:ximum:  (3)  the  critical  angle 
in  the  backscattering  cross-section  data  (see  Secs.IlI-A  and  -B).  The  equivalence  of  the  first 
two  angles  led  Beard  to  suggest  that  a  common  mechanism  was  responsible  for  the  changes  in 
behavior  which  occur  at  this  angle  in  both  the  coherent  and  Incoherent  terms  of  the  forward- 
scattered  energy  and  that  multiple  reflections  might  be  important.  The  equivalence  of  both  these 
angles  with  the  critical  angle  in  the  backscattering  data,  and  the  success  of  the  interference 
concept  in  explaining  this  critical  angle,  lend  strong  support  to  this  statement.  Assuming  that 
the  independent-random-scatterer  model  is  correct  for  both  the  incoherent  forward-scattered 
energy  and  the  backscattered  energy  (all  of  which  is  incoherent),  one  sees  that  the  only  difference 
between  them  should  result  from  the  angular  variation  in  the  power  cross  section  of  an  Individual 
scatterer.  If  one  assumec  that  the  Interference  concept  can  be  applied  to  the  Incoherent  forward- 
scattered  energy  in  the  same  manner  as  it  has  been  applied  to  the  backscattered  energy  (l.e., 
assume  the  incoherent  forward-scattered  energy  arises  from  a  configuration  consisting  of  ran¬ 
dom  scatterors  elevated  above  a  smooth  reflecting  surface),  then  the  data  on  the  Incoherent 
forward -scattered  energy  can  be  processed  (after  appropriate  normalization)  in  exactly  tlie 
same  way  as  the  data  on  the  backscattered  energy.  The  differences  in  the  power  cross  section 
of  an  individual  scatterer  in  the  forward  and  backward  directions  will  result  in  differences  in 

the  value  of  <rL. 

o 

One  very  surprising  result  obtained  on  the  Incoherent  forward-scattered  energy  is  that  this 
energy,  like  tiie  coherent  energy,  appears  to  come  mainly  from  the  first  Fresnel  zone.  Inas¬ 
much  as  the  Fresnel  zones  are  defined  in  terms  of  phase  differences  and  ths  incohsrent  energy 
results  form  a  random  phase  addition,  this  result  constitutes  a  remarkable  coincidence. 

In  addition  to  the  data  on  |Y||  and  |  Y|  -  ?||  ,  this  project  also  obtained  considerable  data 
on  the  power  spectrum  of  the  fluctuations  in  |Y,j.|  >  |Yq  -i  Y^I.  Two  important  results  obtained 
on  this  spectrum  were  (i)  the  spectrum  was  independent  of  polarization  and  (2)  the  spectrum  was 
very  closely  correlated  with  the  spectrum  of  the  ocean  waves. 

For  fiirther  details  on  this  work,  the  reader  is  referred  to  the  paper  by  Beard^^  and  to  the 
references  cited  in  his  paper. 

IV.  SPHERICAL-EARTH  FORMULAS 

In  predicting  or  evaluating  the  performance  of  airborne  radars,  one  is  frequently  faced  with 
a  problem  in  which  the  earth's  curvature  is  a  significant  fact  that  cannot  be  Ignored  without  in¬ 
troducing  serious  errors.  It  has  often  been  assumed  that  Just  because  the  altitudes  of  the  radar 
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and  target  are  small  with  respect  to  the  raaius  of  the  eerth,  this  curvature  is  of  minor  impor¬ 
tance  and  need  not  be  considered.  Actually,  of  course,  no  matter  how  small  the  altitudes  are, 
there  will  always  be  a  horizon,  end  the  curvature  will  always  be  important  in  regions  near  that 
horizon. 

Ac  Functions  to  be  Computed 

Assume  that  the  energy  scattered  forward  to  the  target  can  be  repreoented  by  a  single  in¬ 
direct  ray  reflected  at  the  point  of  specular  reflection.  Let  the  variables  h^,  R,  M^,  1^2, 
dj.  d2.  Tj,  r2.  hj,  h^.  0^,  0^  ‘I'r  ^1*  defined  as  in  Fig.  34,  and  let 

A  -  -t  M2  -  R,  r  =  r^  +  r2.  d  =  d^  +  ^2*  0  =  0^  +  02,  a  =  modified  earth's  radius,  f  =  diver¬ 
gence.  A  =  area  of  illumination,  and  R^  =  horizon  range.  Assume  also  that  the  locations  of  the 
radar  antenna  and  target  are  constrained  to  a  fixed  vertical  plane  and  that  their  orientations  art- 
fixed.  The  variables  which  need  to  be  evaluated  in  order  to  describe  the  total  received  signal 
are  then  R.  A,  y^.  0^,  V2*  l/.  •  'L  «.  A.  and  R^.  The  first  eight  are  needed  in  order  to 

describe  the  target  signal,  the  next  three  in  order  to  describe  the  clutter  signal,  and  the  last  in 
order  to  determine  the  range  at  wh'ch  the  geometric  situation  becomes  degenerate  The  two 
ve'-tical -plane  coordinate  systents  which  ars  most  frequently  used  for  specifying  the  relative 
locations  of  the  radar,  target,  and  earth,  anc  which  serve  as  independen.  variables  for  this 
evaluation,'  are  (hj,  h2,  d)  or  h2,  Rl-  Inasmuch  as  R  and  d  arc  related  by  the  relatively 
simple  (a  ’  e.xact)  formula 

R^  ^  (a  4  4  («  4  h2)^  -  2(-  4  b^)  <«  4  h2)  cos  <  j)  (52) 

wnich  of  Jhese  coordinate  systems  is  actually  used  for  this  purpose  is  of  little  importance.  In 
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addition  to  these  functions,  all  of  which  are  required  in  order  to  describe  the  received  signal  in 
terms  of  the  relative  locations  of  the  radar,  target,  and  earth,  one  may  also  need  to  evaluate 
certain  functions  in  order  to  describe  the  location  of  a  target  in  terms  of  the  received  signal. 
Two  functions  of  potential  interest  in  this  category  are  the  height-finding  functions  h2(hj,  R) 
and  h2(h^,  A.  R). 


B.  Functions  Rjhj,  hj),  «(hj,R),  a(hj,R),  yj(hj,h2,R),  and  A(hj,R) 

The  exact  equations  for  these  functions  are  a^  follows: 

Ro(hj.h2)  =  U.h^  +  +  (2.h,  +  (53) 


«{hj,  R)  =  -in“* 


hj  R^-hj 

■r  *  2R(.  +  h  j) 


«(hj,  R)  =  sin 


yj(h^,h2.R)  =  £in"* 


.+  h2^  -hf  +  (-1)‘  R^ 
2R(  .  +  hj) 


A(hj,  R)  =  cos 


(54) 


(55) 


(56) 


(57) 


In  this  last  equation,  @  denotes  the  horizontal  beamvvidth  of  t!ie  antenna  in  the  direction  of  the 
clutter,  c  denotes  the  velocity  of  light,  and  T  denotes  the  transmitted  pulselength.  This  equa¬ 
tion  assumes  that  A  is  defined  by  the  pulseiength  rather  than  the  vertical  beamwidth  and  that 
the  illuminated  region  consists  of  a  section  of  a  lune.  The  first  factor  is  the  projection  of  0  on 
the  surface  of  the  earth. 

For  many  purposes,  Eqs.  (53)  to  (57)  can  be  replaced  by  the  following  approximations: 


Ho(hj.h2)  =  (2«h^)*/^  +  (2-h2)*^^ 

h^/- 

small 

(58) 

^  (i  * 

hj/« 

,  h^/R  small 

(59) 

hj/ « 

,  h^/R  small 

(60) 

ri<h,,h,,R).  (-.,■•*(  ‘)  [) « 

h./« 

,  ih2  -hj^l/R^  small 

(61) 

A(h^,R)  =  Re  ~  sec  a 

hj/* 

,  RA  ,  cT/4R,  0,  a  small 

.  ibZ) 

«  V  I'’ 


C.  Functions  <^hj,h2tR)>  <('(hj,h2,R),  /Sj(hj,h2,R),  and  $(hj,h2,R) 

These  functions,  all  of  which  involve  the  indirect  rty  to  the  target,  are  much  more  difficult 
to  compute  than  those  considered  in  Sec.  IV-B.  A  numbe’'  of  approaches  to  these  functions  will 
be  considered. 

1.  Exact  Equations 

2  2  2  i/2 

Letting  =  1  +  —  and  =  (a  sin  ^  +  2ah.  +  h.  )  ,  one  can  derive  the  following  exact 

equations: 


dj  =  ^ 

hj  =  («  +  h^)  cos  ©£  -  « 
r^  =  («  +  h^)  sin 

Mi  =  (r.^  + 

A  =  +  (h'^  +  -  (r^  +  (h'^ 

h!  \  ,  /h',  +  h!. 


*  =  tan  *  (-?[)=  tan  (-V^) 


®  = 


(«  +  hj)  (<■  +  h^)  Bin© 


•(M^  +  M2)  eo8<!> 


f  _JMiI^|j 

P  ^  .(M^  +  M2)  sin^^Jj 


•1/2 


(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 


4B^B2  cos^e^  -  (4B^B2^  +  4B^^B2  cog©)  cos^e^  +  -t  cosO 


-43^3^)  coa^e^  +  ain^O  +  -tB^Bg  cos©)  cos©^ 


+  (Bj^B2^  sin^  B  ~  Bz  "  ^®1®2  ^  ° 


©2  =  C03 


-i 


3^  sin^©j  +  (B^  C08©J  -i)  [82(34  -f  1  -2B4  cD8©j)  -Bj^  sin^U^l 


1/2 


B^{3^  +  i  -  2Bj  cos  ©4) 


M.  =  C  - •sin^ 

«i  -  *‘"''(7/5;) 

R  =  (M^  m|  +  2MjM2  cos  2*) 

h^r 

'  (h^ThjnTThp 


(71) 

(72) 

(73) 

(74) 

(75) 

(’6) 
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Equationti  (!>2)  and  (63)  to  (71)  are  sufficient  to  compute  A,  and  $  to  any  desi/ed  accuracy. 

Starting  with  h^,  h2.  and  R.  one  must  solve  (52)  for  d.  d  =  .0  for  0,  (71)  for  Oj,  0  =  0^  +  ©2 
for  ©2.  (64)  to  (66)  for  h|,  r^,  and  M^,  **  =  r^  +  r2  for  r,  and  finally,  (67)  to  (70)  for  A, 
and  S).  In  general,  this  procedure  will  be  very  laborious  since  (71)  is  a  quartic  in  cos©^  A 
much  simpler  procedure  which  can  be  used  if  one  does  not  need  to  evaluate  A,  s^.  and  f  at 
a  specific  point  (h^,  h2,  R)  but  only  to  plot  a  curve  of  these  functiord  for  variable  (h|,  h^.  R),  is 
to  start  off  by  choosing  (h^,h2.  ©|).  compute  by  (72),  8  by  ©  =  0^  +  ©2.  d  by  d  =  •©,  R  by 
(52).  and  then  continue  as  before.  A  second  such  method,  based  on  the  parameter  is  to  be* 
gin  by  choosing  (h^,  ♦).  compute  Mj,  /Jj,,  and  R  by  (73)  to  (75),  A  bv  A  =  +  M2  -  R,  d  by 

(52).  ©  by  d  =  .©,  and  ®  by  (70). 


2.  Approximations  Based  on  Cvltic  in  d^ 

Assuming  that  h^/«  ,  d^/«  ,  and  h^/d^  are  small,  one  can  replace  £qs.  (52),  (64)  t«.  (68),  and 


(70)  to  (72)  by 


R  =  d 


h!  =  h,  - 


d.*- 


i  “i  2- 


ri^di 


M.  =  d. 
1  1 


<l>  ~  tan 


=  tan 


Ihi  +  h2 

,  ‘*2  11 

0 

2dj  -  3dd^  +  (d^  -2*hj  -Z^hj)  d^  +  Z-h^  d  =  0 

d^-Z.h^  |(d2-2.h/ 

2<»1  M  4d2 


If  hj  >  h2.  £q>  (84)  has  the  solution 
d  _ , »  b  w 


d,  = 


where 


i  =  7-P~»(^) 


P=  C-  (hj+h2)+ 


(7?) 

(78) 

(79) 

(80) 

(81) 


(82) 

(83) 

(84) 

(85) 


(86) 
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^  =  cos 


[Z.eKhg  -h 


in  the  rentatnlng  discussion,  whenever  an  equation  based  on  (84)  is  considered,  it  will  be  as> 
sumed  that  h4  >h,.  When  h,.  <  h,.  the  appropriate  equations  can  be  derived  by  interchanging 
1  and  2. 

These  approximations  are  essentially  the  same  as  those  used  at  the  Radiation  Laboratory, 
and  a  considerable  amount  cf  work  has  been  done  with  them  by  many  people.  (See,  for  example. 
Fishback  in  Kerr^^  and.  also.  Burrows  and  Attwood.^^)  Not  only  have  a  number  of  transforma¬ 
tions  been  found  which  normalise  the  variables  in  such  a  way  that  the  einount  of  computing  is 
considerably  reduced,  but  much  of  this  computing  has  already  been  done  and  is  available  in 
graphical  form.  For  example.  Burrows  and  Attwood  Itave  used  the  transformation 

b  =  (dj  -  d2)/(dj  +  dj)  (8V) 

c  =  (h^-h2)/(hj +  h2)  (88) 

m  »  dVd-(hj  +  hj)  (89) 

to  transform  (84)  into  the  simple  cubic 

c  =  b  +  bm(l  —  b^)  (90) 

and  have  plotted  a  graph  relating  b.  m.  and  c  (see  Fig.  35).  The  curved  contour  on  the  right 
side  of  this  graph  is  determined  by  the  horizon  and  represents  the  intersection  of  c  =  b  4 
bm(l  —  b^)  and  c  =  2mb.  In  order  to  compute  A.  g,  and  3) .  one  now  need  only  to  compute 
m  and  c  by  Eqs.  (88)  and  (89).  obtain  b  from  Fig.  (35).  obtain  d^  and  d2  by  £q.  (87)  and  d  - 
d^  d2.  obtain  A.  g.  and  $  by  (81).  (82).  and  (83).  and  obtain  by  d.  -  a6j^  and  (69).  All  these 
operations  can  be  perfoxmied  quite  rapidly.  Burrows  and  Attwood  have  also  used  the  transforms- 


«  "1  (2.h/'‘ 

and  have  plotted  graphs  of  8(u.  v).  A(u.  v).  and  S  (u.  v).  Fishback.  on  the  other  hand,  has  used 
the  transformation 


Letting 


.h  vl/2 

^  =  (sf) 


(2,hj)’'‘  +  (2.h2)’ 


(2.hj)* 


J(S,T)  =  (1-sf)  (1-S|) 


K(S.T)  = 


fl-S^^4T^(l-S2^)l 


one  can  rewrite  Eqs.  (81).  (82).  and  (83)  as 


2h  h, 

A  s  JiS,  T) 

(94) 

.  fh.  +  h,  1 

..•111  ^  1 

\y^f 

^  —  a«iss  ^  ^ 

f  4S?S,T 

5  *  |l  + - ^ - j 

(96) 

1  S(i-Sf)(14T)i 

Graphs  of  J(S,  T),  K<S.  T).  and  P<S.  T)  are  avaUable  both  in  Kerr*^  and  in  Flshbacki** 

3.  Aj^roxiniations  for  SmaU  b2/h^ 

Another  set  of  ai^roximations  to  A,  dr  0^^  wid  S  which  is  of  use  when  considering  targets 
very  close  to  the  surface  is  the  following*  Asannte  that  ,  d^/*,  hj/dj,  and  hj/hj,  are  small. 
Then 


If..  =»j.^  {»- j^) 

(97) 

*1  ’*  =  r 

(98) 

dg  =  -o.  ♦  (o^.^  +  2.hj)‘/* 

(99) 

(100) 

(101) 

In  the  event  that  hj  «  o*  •,  then  dj  reduces  to  dg  *  hg/o,  and  A  and  9  can  be  approximated 
by  A  =  Zohg  and  5*1. 

4.  Series  Expansion  of  h^  in  r/« 

One  ftirther  method  for  approximating  these  functions  is  based  on  the  expension  of  hj  in  the 
variable  r/*: 


r»n 


n*0 

Assuming  is  small.  repUcing  cooGj  by  H  +  tan^Gj)"*/*  in  (64).  expanding  (4  +  tan^ej)**/^ 
in  a  power  series  inTan^Oj,  and  applying  (76)  to  eliminate  tan©j.  one  U  lead  to  the  expression 


**l**»i"f (hj^) 


(402) 


Replacing  h{  on  both  sides  of  this  expression  by  the  above  series  and  emnpoting  the  first  flve 

i  ix\ 


coefficients  one  obuina  the  approximation 
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kl  .  k 


(hj  ♦  h/  2.V,  ^  hj)* 


(103) 


Th*  closer  one  comes  to  the  horison*  the  more  terms  that  wiU  be  needed  for  s  given  degree  of 
sceursey.  The  corresponding  approximations  to  tan#  and  cos 8^  (using  (64)  and  (68)]  are  given 
by 


cos  8j  »  I  - 


hM 

Z.Vh.  +  h,)^ 


(hf  4hj-h^(h^4h^)|h^V 


T'  — r 

2.*(h,  ♦  hj)* 


tan# 


hj  +  hj 


[  2.(h,  th^  J 


(104) 


(105) 


If  (h^  *  h^)/r  is  small  so  that  (6?)  can  be  replaced  by  >  2h|h^/r,  A  can  be  written 


2h^h 


1“2 


1  - 


hihjr- 


2*ihj  ♦  hj)  *  4.*(hj  ♦  hj^ 


(106) 


Eiq>reasions  for  other  variables  can  be  derived  similarly. 

5.  Approximations  for  Small  # 

Using  Eqs.(7S).  (75),  and  A  «  4  M2  —  R,  and  assuming  that  #  and  ((C^  — C2)  tan#|/(C^  4 

C2  *-  2«  sin#)  are  small,  one  can  approximate  R  and  A  by 

R  *  C|  4  C2  —  2«# 


A  *  (ReVz)  (1  -  (Cj  -  Cj^/R*!  . 

If,  in  addition,  hj/«  and  #  «/h|  are  small,  then 

R  »  (2.hj)*/*  4  (2«h2)*/*  -  2.#  4  (.  V/2)  (C2*h/*/*  4  (g.hj)’*/*] 

A  R  t(2  «h^)*/^  4  (2«h2)*^^  -  nf  |(2>h^)*/^~(2«h2)*^Yj 


(107) 

(108) 

(109) 

(110) 


Equations  (107)  and  (108)  are  accurate  over  a  wide  range  of  interest.  Equations  (109)  and  (110) 
are  of  ttsc  only  in  regions  extremely  close  to  the  iMuiscn. 

D.  fbnetloosh2(b|,r|,R)aiidli|{h|,A,R) 

1.  h2(hj,  R) 

The  exact  equation  for  this  Ibnction  is  given  by 

hj  *  "•  •  ♦  ({•  4  hjl?  4  R*  4  2R{ .  4  hj)  siny,]*'^*  .  (Ill) 

If  hj/« ,  h|/R  and  (R^  4  2«  (h|  4  R  siny^)]/*  ^  are  small,  then 
hj  *  hj  4  R  sinVj  4  R*/2« 
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(112) 


1.  h,(h..A.R) 

\ 

Referrinff  to  ths*  tranflfnrmrition  /Qt\  unrl  rfAfinin#*  <i*  Kw  Q*  -  aIi‘>  -  w  c  -  ct  / 

V  -  - - - - -  ^  .  w  ^  \yS(V  licto  *»»  »  tJ  / 

z  *  z  * 

(1  +  T).  Thus,  N  =  T'  J(S,  T)  -  Ad/2h^  can  be  regarded  as  a  function  of  S'  and  T.  Since  h^, 

A,  and  d  deterr'ine  N  and  S',  the  problem  of  determining  A,  R)  is  equivalent  to  solving 

N(S',  T)  for  T.  A  graphical  solution  of  N(S'.  T)  is  available  in  Fishback'*^ 

Another  method  for  determining  h^lhj,  A.  R)  has  been  aerived  by  McCracken.^^  Substituting 

^2  *  ^  ®n“  into  (106)  and  evaluating  the  coefficients  a  ,  one  obtains  the  approximation 

n=0 


For  a  more  comprehensive  discussion  of  spherical-earth  formulas,  the  reader  is  referred 

44 

to  the  report  by  Durlach,  et  al.  This  report  contains  additional  exact  equations,  additional 
approximations,  detailed  graphical  comparisons  of  the  approximations  with  the  exact  equations, 
and  finally,  a  large  number  of  graphs  computed  by  use  of  the  exact  equations. 
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